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ABSTRACT 

In  this  program  we  are  investigated  the  dynamics  of  cavitation  as  it  occurs  on  or  near  rotating  propulsors.^  While 
there  are  many  different  forms  of  cavitation  on  propulsors,  we  are  concentrating  our  effort  on  understanding  the 
flow  processes  resulting  in  "limited  event  rate"  cavitation  occurring  in  the  vortical  flows  in  the  tip  region  of  ducted 
propulsors.  We  completed  the  construction  of  a  new  recirculating  water  channel.  We  are  examining  the  process 
of  nuche  capture,  cavitation,  and  noise  production  of  vortex  cavitation.  We  have  been  conducting  experirnents  both 
at  our  laboratory  at  the  University  of  Michigan  and  at  the  Naval  Surface  Warfare  Center,  Carderock  Division.  Our 
primary  effort  has  been  examining  the  tip  leakage  flow  of  a  large  ducted  rotor  in  the  36mch  water  tunnel  at 
NSWC-CD.  This  experiment  is  conducted  in  cooperation  with  Dr.  Jessup  of  NSWC. 


RESULTS 

The  main  findings  of  this  project  are  contained  in  the  four  attached  manuscripts.  Three  other  manuscirpts  are  in 
preparation. 
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The  influence  of  developed  cavitation  on  the  flow  of  a  turbulent  shear  layer 

Claudia  O.  Iyer  and  Steven  L.  Ceccio 

Department  of  Mechanical  Engineering,  University  of  Michigan,  Ann  Arbor,  Michigan  49109-2121 
(Received  28  January  2002;  accepted  27  June  2002;  published  28  August  2002) 

Developed  cavitation  in  a  shear  layer  was  studied  experimentally  in  order  to  determine  the  effect 
that  the  growth  and  collapse  of  cavitation  have  on  the  dynamics  of  shear  flows.  Planar  particle 
imaging  velocimetry  (PIV)  was  used  to  measure  the  velocity  field,  the  vorticity,  strain  rates,  and 
Reynolds  stresses  of  the  flow  downstream  of  the  cavitating  and  noncavitating  shear  layer;  the  flow 
pressures  and  void  fraction  were  also  measured.  The  flow  downstream  of  a  cavitating  shear  flow  was 
compared  to  the  noncavitating  shear  flow.  For  cavitating  shear  layers  with  void  fractions  of  up  to 
1.5%,  the  growth  rate  of  the  shear  layer  and  the  mean  flow  downstream  of  the  shear  layer  were 
modified  by  the  growth  and  collapse  of  cavitation  bubbles.  The  cross-stream  velocity  fluctuations 
and  the  Reynolds  stresses  measured  downstream  of  the  cavitating  shear  layer  were  reduced 
compared  to  the  entirely  noncavitating  flow.  This  result  is  inconsistent  with  a  scaling  of  the  shear 
stress  within  the  shear  flow  based  on  the  mean  flow.  The  decrease  in  the  cross-stream  fluctuations 
and  Reynolds  stresses  suggests  that  the  cavitation  within  the  cores  of  strong  streamwise  vortices  has 
decreased  the  coupling  between  the  streamwise  and  cross-stream  velocity  fluctuations.  ©  2002 
American  Institute  of  Physics,  [DOI:  10.1 063/1 . 1501541] 


I.  INTRODUCTION 

Cavitation  has  been  observed  to  occur  in  turbulent  shear 
flows  when  the  local  liquid  pressure  drops  below  the  liquid 
vapor  pressure  and  vapor  bubbles  form  in  the  low-pressure 
cores  of  fluid  vortices.  Such  “vortex  cavitation”  can  occur 
within  strong  liquid  jets,  within  regions  of  flow  separation, 
and  in  the  turbulent  wakes  of  bluff  objects  (see  Brennen^  for 
a  recent  review).  Cavitating  shear  flows  occur  in  the  closure 
region  of  attached  cavities  and  are  associated  with  the  for¬ 
mation  of  strong  vortices  downstream  of  a  partial  cavity 
(Laberteaux  and  Ceccio,^  Laberteaux  et  Gopalan  and 
Katz,"*  and  Laberteaux  and  Ceccio^’^).  The  cavitating  turbu¬ 
lent  shear  flow  associated  with  attached  cavitation  plays  an 
important  role  in  the  formation  of  cloud  cavitation  and  in  the 
processes  that  cause  cavitation  erosion  (Kawanami  et  al.,^ 
Callenaere  et  al}). 

Developed  cavitating  shear  layers  are  complex  multi¬ 
phase  flows.  The  continuous  flow  field  is  turbulent  and  char¬ 
acterized  by  relatively  high  Reynolds  numbers  (typically 
greater  than  10000  based  on  integral  length  scales).  The  dis¬ 
perse  (vapor)  phase  has  a  complex,  nonspherical  geometry 
and  is  often  not  in  thermodynamic  equilibrium  with  the  con¬ 
tinuous  phase.  Volume  fractions  can  be  relatively  large  (often 
much  greater  than  1%).  The  vapor  phase  is  often  preferen¬ 
tially  concentrated  within  the  low-pressure  cores  of  the  stron¬ 
gest  vortices,  the  sites  of  cavitation  inception,  and  the  pres¬ 
ence  of  vapor  in  the  shear  layer  can  alter  both  the  large-scale 
and  small-scale  flow  processes.  But,  to  date,  experimental 
observations  have  not  revealed  which  interactions  between 
the  vapor  and  the  liquid  flow  fields  constitute  the  dominant 
mechanisms  responsible  for  overall  flow  modification.  This 
question  will  be  addressed  by  the  present  study  through  the 
experimental  examination  of  a  cavitating  shear  layer. 

Large-scale  spanwise  vortical  structures  are  an  important 


feature  of  turbulent  shear  layers  (Brown  and  Roshko,^ 
Winant  and  Brown**^).  Counter-rotating  pairs  of  streamwise 
vortices  can  occur  between  the  larger  spanwise  vortices 
(Bernal  et  Lasheras  et  and  Lasheras  and  Choi^^), 
With  a  reduction  in  free-stream  pressure,  the  cores  of  the 
spanwise  and  streamwise  vortices  can  fill  with  gas  and  vapor. 
Cavitating  shear  flows  have  been  investigated  by  a  number 
of  researchers,  and  a  review  of  the  older  literature  is  pro¬ 
vided  by  Arndt.  Katz*^  studied  the  cavitation  inception 
within  shear  layers  formed  on  axisymmetric  bodies  whose 
boundary  layers  underwent  a  laminar  separation  and  subse¬ 
quent  turbulent  reattachment.  Katz  and  O’Hera^^  and 
O’Hem^^  examined  cavitation  inception  and  development  in 
a  plane  shear  layer  for  Reynolds  numbers  up  to  2  X  10^.  Katz 
and  O’Hem*^  showed  that  cavitation  inception  occurs  first 
within  the  streamwise  vortices  of  the  shear  layer.  The  stream- 
wise  vortices  are  stretched  between  the  spanwise  vortices, 
resulting  in  a  reduction  in  the  vortex  diameter  and  core  pres¬ 
sure.  A  further  reduction  in  the  overall  flow  pressure  will 
then  lead  to  cavitation  in  the  streamwise  vortices.  In  a  con¬ 
tinuation  of  this  study,  O’Hem^^  did  not  observe  significant 
changes  in  the  shear  layer  growth  rate  of  a  vortex  roll-up 
process  as  the  cavitation  developed. 

Cavitation  forming  in  the  turbulent  wake  of  bluff  objects 
has  been  investigated  by  many  researchers,  including 
Kermeen  and  Parkin,*^  Young  and  Holl,*^  Amdt,^®  Arndt  and 
George, and  Belahadji  et  Shedding  of  wake  vortices 
can  be  altered  by  the  development  of  cavitation  in  the  vortex 
cores,  and  developed  wake  cavitation  will  typically  increase 
the  shedding  frequency  of  the  Karman  vortices  observed  in 
the  wake  by  up  to  25%.  Belahadji  et  alP  employed  high¬ 
speed  photography  to  investigate  cavitation  in  the  rotational 
structures  of  the  turbulent  wake  behind  a  two-dimensional 
obstacle.  They  concluded  that  incipient  cavitation  in  the 
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wake  acted  as  a  passive  agent  of  flow  visualization,  but  de¬ 
veloped  cavitation  actively  altered  the  dynamics  of  the  tur¬ 
bulent  wake. 

The  presence  of  cavitation  within  a  turbulent  shear  flow 
can  potentially  alter  the  dynamics  of  the  flow  through  both 
local  and  global  mechanisms.  The  presence  of  developed 
cavitation  can  lead  to  overall  flow  changes  through  decreases 
in  mean  flow  density,  increases  in  the  overall  flow  volume,  or 
modification  of  the  mean  pressure  field.  Cavitation  can  also 
modify  the  vortical  flow  locally.  Belahadji  et  alP'  have  sug¬ 
gested  that  cavitation  in  the  cores  of  vortices  will  affect  the 
process  of  vortex  stretching  by  the  decoupling  of  the  vortex 
strain  and  rotation  rate.  Gopalan  and  Katz"^  demonstrated  that 
significant  flow-vapor  interactions  exist  in  the  cavitating 
shear  flow  downstream  of  a  partial  cavity.  And,  Laberteaux 
and  Ceccio^’^  showed  how  the  growth  and  collapse  of  dis¬ 
perse  cavitation  bubbles  led  to  the  production  of  small-scale 
turbulence  in  the  wake  of  partial  cavities. 

In  the  present  study,  we  investigated  how  the  develop¬ 
ment  and  collapse  of  cavitation  within  a  plane  shear  layer 
alters  the  dynamics  of  the  shear  layer.  The  shear  flow  was 
examined  from  cavitation  inception  to  full  development.  Pla¬ 
nar  particle  imaging  velocimetry  (PIV)  was  used  to  examine 
the  flow  downstream  of  the  cavitating  shear  layer,  and  flow 
pressures  and  void  fractions  were  measured.  We  will  show 
that  the  growth  and  collapse  of  cavitation  within  the  shear 
layer  altered  the  mean  and  fluctuating  velocities  of  the  flow 
downstream  of  the  cavitating  region, 

II.  EXPERIMENTAL  SETUP  AND  METHODS 
A.  Flow  facility  and  test  model 

Experiments  were  conducted  in  the  Model  Large  Cavi¬ 
tation  Channel  (MLCC)  in  the  Cavitation  and  Multiphase 
Flow  Laboratory  described  by  George  and  Ceccio.^  The 
MLCC  is  a  recirculating  water  tunnel  with  a  capacity  of  ap¬ 
proximately  2.8  m^.  Most  of  the  lower  leg  has  circular  cross 
sections  upstream  and  downstream  of  the  impeller;  the  upper 
leg,  including  the  contraction,  test  section,  and  diffuser,  has  a 
square  cross  section.  The  flat  top  of  the  MLCC  allows  low 
pressures  in  the  test  section  without  subvapor  pressures  and 
cavitation  occurring  in  adjacent  parts  of  the  tunnel.  As  a 
result,  the  6:1  contraction  upstream  of  the  test  section  is  non- 
symmetric,  and  internal  flow  straighteners  are  angled  upward 
to  make  the  velocity  profile  uniform  within  the  test  section. 
The  rectangular  test  section  has  interior  dimensions  of  21.6 
X21.6X  1 18  cm,  with  a  flow  area  of  446  cm^  due  to  fillets  in 
the  comers.  The  four  sides  of  the  test  section  hold  removable 
windows  made  of  acrylic  or  aluminum.  Before  experiments, 
water  is  deaerated  by  circulation  between  the  MLCC  and  a 
2.27  m^  deaeration  tank;  the  oxygen  content  of  the  water  can 
be  reduced  to  about  40%  saturation  in  three  hours.  Flow  in 
the  MLCC  is  generated  by  an  impeller  in  the  bottom  leg.  The 
MLCC  can  operate  at  constant  test  section  velocities  of  up  to 
15  m/s.  The  overall  air  content  of  the  MLCC  can  be  modified 
using  a  deaeration  system.  However,  the  nuclei  content  of  the 
flow  within  the  tunnel  cannot  be  directly  controlled.  The 


oxygen  content  of  the  water  was  measured  with  a  dissolved 
oxygen  meter  (Orion  Model  810).  Values  of  the  oxygen  con¬ 
tent  are  reported  at  atmospheric  pressure. 

To  create  a  shear  layer,  a  sharp-edged  plate  spanning  the 
full  21.6  cm  width  of  the  test  section  was  placed  vertically 
into  the  flow  in  a  manner  similar  to  that  of  Katz  and 
O’Hem*^  and  O’Hem.^^  The  plate  was  90  mm  from  the  en¬ 
trance  to  the  test  section,  and  mounted  on  a  false  floor  to 
allow  the  boundary  layer  flow  on  the  test  section  walls  to 
pass  beneath  through  a  gap  12.5  mm  tall.  The  false  floor 
consists  of  a  58  mm  tall  and  710  mm  long  wedge.  The  height 
of  the  plate,  h,  is  equal  to  31  mm.  The  top  of  the  plate 
consists  of  a  chamfered  edge  with  a  broken  apex.  The  shear 
layer  produced  at  the  apex  of  the  wedge  flows  into  a  diffuser 
with  an  angle  of  9°.  The  pressure  recovery  in  this  diffuser  is 
sufficient  to  cause  desinence  of  the  cavitating  shear  layer.  A 
schematic  diagram  of  the  plate  within  the  test  section  is 
shown  in  Fig,  1, 

A  differential  pressure  transducer  (Mid-West  Instrument 
Model  120)  was  used  to  measure  the  pressure  between  the 
entrance  to  the  contraction  and  the  wall  above  the  apex  of  the 
plate.  This  pressure  difference  was  used  to  determine  the 
free-stream  velocity,  f/oot  after  calibration  with  a  single¬ 
component  laser  Doppler  anemometry  (LDA)  system  (Dan- 
tec  DISA  TVpe  55L90a)  with  the  probe  volume  located  mid¬ 
way  between  the  apex  of  the  plate  and  the  top  window  of  the 
test  section.  This  tap  was  also  used  to  measure  the  static 
pressure  at  the  inlet  of  the  test  section  with  a  Test  Gauge  by 
Ashcroft.  Pressure  was  measured  at  this  tap  to  determine  the 
cavitation  number,  o*=(P/”Py)/5-pf/L  The  temperature  of 
the  flow  was  measured  with  a  thermocouple  and  used  to 
determine  the  liquid  density,  p,  liquid  kinematic  viscosity,  v, 
and  vapor  pressure,  .  A  second  pressure  tap  was  located 
downstream  of  the  plate,  and  the  differential  pressure  be¬ 
tween  the  upstream  and  downstream  taps  was  measured  with 
a  differential  pressure  transducer  (Omega  PX26-005DV). 
Moreover,  dynamic  pressure  transducers  were  placed  near 
these  two  pressure  taps  (PCB  Piezotronics  Model  102A05) 
to  measure  the  fluctuating  pressure  component.  Spectra  of 
the  fluctuating  pressures  were  measured  with  a  Stanford  Re¬ 
search  Systems  SR770  FFT  network  analyzer.  A  bandwidth 
of  1  Hz  was  used  over  a  range  of  1024  kHz,  and  1000  indi¬ 
vidual  spectra  were  averaged.  Side  view  and  top  view  images 
of  the  shear  layer  were  recorded  with  a  Nikon  35  mm  camera 
using  flash  photography  and  with  a  high-speed  digital  imag¬ 
ing  system  (Kodak  EktaPro  HS  Motion  Analyzer  Model 
4540mx). 

B,  Void  fraction  measurement 

The  average  void  fraction  of  the  shear  layer  was  mea¬ 
sured  using  two  different  techniques.  Light  scattering  was 
used  to  detect  the  passage  of  bubbles  in  a  small  probe  vol¬ 
ume.  The  LDA  system  was  used  to  detect  the  passage  of 
bubbles  within  its  probe  volume.  As  bubbles  passed  through 
the  probe  volume  (with  dimensions  of  1 ,0  mm  for  the  probe 
volume  length  and  0.12  mm  for  the  beam  waist  diameter),  a 
strong  flash  of  laser  light  was  detected  by  the  photomultiplier 
tube.  This  signal  saturated  the  tube,  creating  an  electrical 
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pulse  with  time  width  much  longer  that  the  actual  bubble 
passage.  Thus,  the  light  scattering  system  could  detect  the 
passage  of  a  bubble,  but  not  the  residence  time  within  the 
probe  volume.  Instead,  the  residence  time  was  measured 
with  a  high-speed  digital  imaging  system  (Kodak  EktaPro 
HS  Motion  Analyzer  Model  4540mx).  The  duration  of  the 
light  pulses  that  resulted  from  the  passage  of  bubbles  within 
the  probe  volume  was  recorded  and  was  typically  250  yus. 
Bubble  event  rates  were  measured  for  time  intervals  of  60  s 
and  then  averaged.  With  the  bubble  event  rate,  /q  ,  and  resi¬ 
dence  time,  tgj  the  local  void  fraction,  a,  can  be  calculated 
as 

ot-fsh' 

The  above  equation  provides  a  rough  estimate  of  the 
void  fraction.  Second,  a  laser  light  sheet  was  used  to  illumi¬ 
nate  a  plane  of  the  cavitating  shear  layer  and  digital  images 
of  the  shear  layer  were  recorded.  (More  details  of  the  light 
sheet  and  camera  setup  are  provided  below.)  The  image  of 
the  shear  layer  contained  images  of  the  bubbles  within  the 
shear  layer  with  some  in  focus  and  others  out  of  focus.  The 
image  was  thresholded  and  scaled  to  produce  a  binarized 
image  that  recovered  only  the  in-focus  bubbles.  A  pixel  value 
of  “0”  was  liquid  and  a  value  of  “1”  was  gas.  Five-hundred 
images  were  averaged  to  yield  the  overall  average  pixel 
value  for  the  image  and  this  was  considered  the  void  frac¬ 
tion.  The  void  fraction  values  computed  this  way  were  con¬ 
sistent  with  the  values  obtained  from  light  scattering,  al¬ 
though  the  uncertainty  of  both  measurements  is  relatively 
large. 

C.  Particle  imaging  velocimetry 

The  noncavitating  and  cavitating  flows  were  examined 
with  planar  particle  imaging  velocimetry  (PIV).  PIV  mea¬ 


surements  were  performed  in  region  A,  upstream,  for  non¬ 
cavitating  flows  and  in  region  B,  downstream,  for  both  non¬ 
cavitating  and  cavitating  flows,  as  shown  in  Fig.  1.  Two 
frequency-doubled  Quanta-Ray  Spectra-Physics  GCR-130 
Nd:YAG  lasers  were  used  to  create  a  double -pulsed  light 
sheet  in  the  test  section.  These  lasers  emit  light  pulses  with 
10  ns  duration  at  532  nm  wavelength  with  a  repetitive  rate  of 
30  Hz  and  30  mJ/pulse,  The  light  sheet  was  approximately  1 
mm  thick  and  100  mm  wide  and  was  placed  parallel  to  the 
mean  flow  direction.  Single-frame,  double-pulsed  images 
were  recorded  at  90°  to  the  light  sheet  plane  with  a  Silicon 
Mountain  Design  SMD  4M4  digital  camera  that  frames  at  4 
Hz  and  has  a  resolution  of  2048  X  2048  pixels.  The  camera 
was  mounted  on  an  X-T  traverse  for  easy  translation  along 
the  test  section.  The  light  sheet  making  optics  was  also 
placed  on  a  traverse  that  could  be  translated  from  region  A  to 
region  B  of  the  test  section.  The  camera  has  an  image- 
shifting  feature  whereby  the  image  created  by  the  first  laser 
pulse  is  shifted  vertically  by  a  finite  number  of  lines  before 
the  second  image  is  recorded.  This  image  shifting  permits  the 
introduction  of  a  velocity  bias  in  the  double-pulsed  image, 
thus  eliminating  ambiguity  in  the  sign  of  the  velocity.  Images 
were  acquired  at  3  Hz,  and  the  time  between  light  pulses,  Ar, 
was  measured  within  an  accuracy  of  50  ns  with  an  HP  54620 
A  logic  analyzer.  The  time  separation,  Ar,  determines  the 
maximum  and  minimum  velocities  that  can  be  measured. 
The  time  separation,  Ar,  used  in  region  A,  upstream,  was 
21.6  yus  and  in  region  B,  downstream,  53,6  yots.  PIV  measure¬ 
ments  were  made  with  a  field  of  view  of  5,5  X  5.5  cm.  Spatial 
calibration  of  the  image  was  performed  by  registering  the 
image  on  a  fixed  object  after  a  known  motion  on  the  traverse, 
and  the  traverse  motion  was  measured  to  an  accuracy  of  50 
yum  over  a  span  of  55  000  yum.  The  flow  was  seeded  with 
fluorescent  latex  particles  80  fjm  in  diameter  with  a  specific 
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gravity  between  0.9  and  1.1.  A  filter  was  placed  on  the  cam- 
era  to  block  the  green  light  of  the  lasers  (\  =  532nm)  and 
only  allow  the  fluorescent  return  of  the  particles  (\ 
=  580  nm).  In  this  way,  images  of  vapor  bubbles  could  be 
excluded  from  the  PIV  image. 

The  double-pulsed  images  were  interrogated  using  auto¬ 
correlation  with  the  software  package  visiFLOW  by  AEA 
Technology.  An  interrogation  region  of  128  X 128  pixels  was 
used  with  an  overlap  of  50%,  resulting  in  a  resolution  of  1.7 
mm  with  a  field  of  view  of  5.5 X  5.5  cm.  There  were  typi¬ 
cally  8-15  particle  pairs  within  each  interrogation  region. 
With  a  time  separation  of  21.6  fis  in  region  A  and  53.6  pis  in 
region  B,  a  velocity  of  5  m/s  would  correspond  to  a  particle 
displacement  of  four  pixels  and  ten  pixels,  respectively.  The 
number  of  bad  vectors  in  each  image  was  typically  5%  or 
less  and  interpolation  was  used  to  replace  them.  Each  pro¬ 
cessed  image  yielded  32X32  instantaneous  local  velocity 
vectors  in  the  plane  of  interest.  The  associated  vorticity  and 
strain-rate  distributions  were  computed  from  the  two- 
dimensional  velocity  field  using  centered  finite  differences: 


dv  du 
dx  dy 


2^x 


2Ay 


SxyiiJ) 


dv  du  v{i+\J)  —  v{i  —  \,j) 
'dx^'dy~  2Kx 


2Ay 


where  u  and  v  are  the  horizontal  and  vertical  velocity  com¬ 
ponents,  i  and  j  are  the  horizontal  and  vertical  indices,  and 
Ax  and  Ay  are  the  horizontal  and  vertical  grid  spacing.  Up 
to  1000  images  were  averaged  to  find  the  average  and  fluc¬ 
tuating  component  of  these  quantities. 


D.  Experimental  uncertainty 

The  estimated  uncertainty  in  the  measured  mean  pres¬ 
sures  and  velocity  is  ±2%.  The  temperature  was  accurate  to 
within  ±0.2  °C.  This  leads  to  an  uncertainty  in  the  cavitation 
number  of  ±6%.  The  uncertainty  of  the  PIV  velocity  is  re¬ 
lated  to  the  uncertainty  in  the  spatial  and  temporal  measure¬ 
ments  for  the  system.  The  uncertainty  of  the  spatial  measure¬ 
ment  is  related  to  the  magnification  of  the  image  area,  the 
particle  displacement,  and  the  size  of  the  PIV  autocorrelation 
interrogation  area.  For  a  speed  of  5  m/s,  the  average  particle 
displacement  was  10  pixels.  The  displacements  are  known  to 
subpixel  resolution  of  ±0.2  pixels.  This  yields  an  uncertainty 
in  the  displacements  of  around  ±2%.  The  uncertainty  is 
larger  for  the  cross-stream  velocities  since  they  are  on  the 
order  of  0.5  m/s  or  less.  The  uncertainty  in  the  spatial  cali¬ 
bration  is  0.1%  (based  on  the  traverse  accuracy  of  50  fxm 
over  a  span  of  55  000  /.tm).  The  time  separation  between 
successive  laser  light  pulses  is  measured  very  accurately 
within  50  ns,  which  makes  the  uncertainty  of  the  measured 
time  interval  less  than  ±0.1%.  Thus,  we  estimate  that  the 
uncertainty  in  the  measured  velocity  is  approximately 


1  %-2%  for  the  streamwise  («  component)  of  the  velocities, 
and  10%-20%  for  the  cross-stream  {v  component)  of  the 
velocities.  The  uncertainty  in  the  measurement  of  the  aver¬ 
age  void  fraction  is  about  25%  of  the  reported  average  value. 

The  smallest  resolvable  velocity  fluctuation,  Aw/Ufj,  is 
±0.02  based  on  the  minimum  resolution  of  the  PIV  process¬ 
ing  of  ±0.2  pixels.  Here  N=  1000  images  were  averaged  to 
achieve  .a  higher  level  of  precision,  and  the  resulting  uncer¬ 
tainty  is  approximated  by  Au/Uf^'(\/ yfN) ~  ± 0.00 1 .  Differ¬ 
ences  of  less  than  ±0.002  are  within  the  uncertainty  of  the 
measurement.  If  the  integral  length  scale  of  the  flow  is 
'^2<5=L/,  the  spatial  resolution  of  the  PIV  images  is 
~0.03*L/.  This  compares  with  an  estimate  of  the  Taylor 
length  scale  L),<-‘L,^(v/k''%)--^L,^{vl0.2- Uf,-2S) 
^  0.003 -L/,  where  the  turbulent  kinetic  energy,  k,  is  ap¬ 
proximated  with  ^/I^v  ^0.2t/f5.  Consequently,  these 
measurements  do  not  spatially  resolve  the  smallest  turbulent 
length  scales  of  the  flow  and  the  implication  for  not  resolv¬ 
ing  all  scales  is  that  we  will  likely  underestimate  the  local 
turbulent  intensity. 

The  smallest  resolvable  Reynolds  stress,  Au^u/U\,  is 
±0.04,  based  on  the  minimum  resolution  of  the  PIV  process¬ 
ing  of  ±0.2  pixels.  The  V=1000  images  were  averaged  to 
achieve  a  higher  level  of  precision,  and  the  resulting  uncer¬ 
tainty  is  approximated  by  A(^uAu/U%)  =  ^uAu/U%’(\/ 
VV)~ ±0.001.  Differences  of  less  than  ±0.002  are  within 
the  uncertainty  of  the  measurement. 

III.  RESULTS  AND  DISCUSSION 

A.  The  non-cavitating  flow  of  the  shear  layer  in 
region  A 

PIV  was  used  to  examine  the  noncavitating  flow  of  the 
shear  layer  near  the  origin  of  the  shear  layer  in  region  A.  The 
free-stream  velocity  measured  above  the  shear  layer  was 
f/M=9m/s,  giving  a  Reynolds  number  based  on  the  plate 
height,  /z  =  31  mm,  of  Re =2hU^/v— 5.56X10^  (O’Hem*^ 
reported  data  for  Re-2hU„lv=9.57x  10^).  The  free-stream 
pressure  was  sufficiently  high  to  prevent  cavitation  (a 
=  2.3).  Figure  2  shows  (a)  the  mean  in-plane  velocity,  u 
“^Mean,A»  W  the  mean  vorticity,  (c)  the  mean  strain 
rate,  Sj^y ,  and  (d)  the  mean  Reynolds  stress,  u'v' 

Here  and  is  the  average  velocity  deter¬ 

mined  by  averaging  all  32X32  vectors  of  the  mean  velocity 
field  of  region  A.  700  individual  vector  realizations  were 
averaged.  As  expected,  high  levels  of  vorticity,  strain  rate, 
and  Reynolds  stresses  are  observed  in  the  middle  of  the  shear 
layer. 

The  data  from  Fig.  2  was  used  to  determine  the  nondi- 
mensional  mean-axial  velocity  profiles  (u 
-t/Bonom,A)/(t/Top,A-t^Bottom,A)-  Thcse  BTC  plotted  Bgainst 
the  nondimensional  vertical  distance,  y/<5,  at  three  cross- 
sections  of  region  A;  x/h  =  5.53,  x/h  =  5.96,  and  x/h  =  6.4,  in 
Fig.  3.  The  y  axis  was  first  shifted  by  <5 such  that  the  origin  of 
the  y  axis  corresponds  to  the  point  where  (u 

-t/Bottora,A)/(t/Top,A-t/Botton.,A)  =  0.5  and  then  the  >>  values 

were  nondimensionalized  by  S=  26.2  mm,  corresponding  to 
approximately  half  of  the  local  shear  layer  thickness.  Here 
^Top.A^^.O  m/s  and  UBoitom,A~0.5  m/s  are  the  average  ve- 
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FIG.  2.  Average  noncavitating  flow  in  region  A  presenting  (a)  mean  velocity,  u- Umcwa  (^s),  (b)  vorticity,  (1/s),  (c)  strain  rate,  S^y  (I/s);  and  (d)  the 
Reynolds  stresses  -m'o  Vt/Mean.A*  Th®  from  left  to  right.  Here  1000  individual  images  were  averaged. 


locities  of  the  high-  and  low-speed  free-stream  flows  mea¬ 
sured  at  the  top  and  bottom  boundary  of  the  PIV  images  in 
region  A,  Three  averaged  velocity  profiles  across  the  shear 
layer  are  shown  for  three  distances  measured  from  the  vertex 
of  the  plate.  In  this  figure  the  downstream  distance  x  was 
nondimensionalized  by  the  plate  height  h,  (This  length  scale 
was  similarly  chosen  by  O’Hem*^  and  it  is  on  the  order  of 
the  vortex  spacing  in  region  A.) 

The  velocity  profiles  indicate  that  the  shear  layer  is  close 
to  achieving  a  similarity  condition.  Using  the  average  shear 
layer  thickness,  2S,  in  region  A  and  the  velocity  difference, 
(^Top.A““f^Bottom.A)»  l^e  average  Reynolds  number  is  calcu- 
lated  as  Re5=2^l7xop,A-  i^Bottom,A)/»'=3.35x  10*.  Wygnan- 
ski  and  Fiedler^'*  used  hot-wire  probes  to  examine  a  two- 
dimensional,  incompressible  shear  layer  of  air  forming 
downstream  of  a  backward  facing  step.  This  flow  geometry 
is  analogous  to  that  of  the  present  study.  Their  flow  had  a 
local  Reynolds  number  on  the  order  of  Re/^IO"^  after  the 
shear  layer  achieved  similarity.  Figure  4  presents  the  distri- 
bution  of  the  nondimensional  velocity  fluctuation  profiles, 

U  U  /(  C/xop.A““  ^BoUom.A)  2^nd  U  U  /( U-pop. A  ~  ^Bottom, a)  » 


FIG.  3.  Mean  axial  velocity  profiles  in  region  A  derived  from  the  data  in 
Fig.  2.  Here  («- f/Botiom.A)^(t^Top.A- ^/Boitom,A)  is  plotted  across  the  shear 
layer  for  three  downstream  positions,  x/h.  The  symbols  O,  •,  and  A  indi¬ 
cate  x//i=  5.53,  5.97,  and  6.40,  respectively. 
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FIG.  4.  Velocity  fluctuation  profiles  (a)  (b) 

u't;V(£/Top  A-f/Boiiom.A)^  “  region  A  derived  from  the  data  in  Fig.  2  for 
three  downstream  positions,  xlh.  The  symbols  O,  #,  and  A  indicate  xlh 
=  5.53,  5.97,  and  6.40,  respectively. 

three  cross-sections.  The  maximum  velocity  fluctuation  am¬ 
plitude  is  4.4%,  which  qualitatively  agrees  with  a  maximum 
value  of  18%  reported  by  Wygnanski  and  Fiedler^'^  for  the 
nondimensional  root-mean-square  velocity  fluctuations.  It 
should  also  be  noted  that  the  turbulence  level  does  not  tend 
to  zero  outside  of  the  shear  layer,  and  this  suggests  that  the 
free -stream  turbulence  levels  are  relatively  high  (on  the  order 
of  2%  of  the  free-stream),  even  with  the  flow  conditioning 
and  6:1  contraction.  This  indicates  that  the  shear  layer  in 
region  A  is  reasonably  two-dimensional  and  has  becomes 
self-similar  before  reaching  the  diffuser  in  region  B. 

B.  Examination  of  the  cavitating  shear  layer 

While  maintaining  the  free-stream  velocity  at  10  m/s,  the 
free-stream  pressure  was  reduced,  and  cavitation  was  ob¬ 
served  in  region  A.  The  inception  cavitation  number  was 
approximately  £r=0.65  for  a  dissolved  oxygen  content  of  3 


ppm  and  Gr=0.57  for  a  dissolved  oxygen  content  of  1.7  ppm. 
O’Hem^^  reported  a  cavitation  inception  index  of  cr=  1.1  for 
a  dissolved  air  content  of  3.6  ppm  and  cr=  1.4  for  an  air 
content  of  7-12  ppm.  The  extent  of  cavitation  developed  at  3 
ppm  is  similar  to  that  reported  by  O’Hem.*^  (See  Yu  and 
Ceccio^^  for  a  discussion  of  the  relationship  between  dis¬ 
solved  oxygen  content  and  dissolved  air  content.)  In  both  the 
present  experiment  and  those  of  O’ Hem,  the  free-stream 
nuclei  content  was  not  controlled  or  measured.  These  differ¬ 
ences  in  inception  indices  between  the  work  of  O’ Hem  and 
the  present  work  can  be  partially  explained  by  the  free- 
stream  nuclei  content.  However,  our  method  used  to  call  in¬ 
ception  in  the  present  work  is  “visual  calls’’  and  it  leads  to  a 
delayed  cavitation  call  compared  to  the  method  of  O’ Hem. 

Figure  5  presents  side-view  photographs,  taken  under 
stroboscopic  light,  of  different  levels  of  cavitation  in  the 
shear  layer  (with  an  oxygen  content  of  2.7  ppm).  Figure  5(a) 
shows  that  inception  cavitation  occurs  in  the  form  of  long, 
thin  streamwise  cavities  oriented  at  about  45°  to  the  mean 
flow,  visualizing  the  low-pressure  cores  of  streamwise  vorti¬ 
ces  of  the  shear  layer  that  form  between  adjacent  span  wise 
vortices.  At  lower  cavitation  numbers,  the  entire  turbulent 
stmcture  of  the  shear  layer  can  be  visualized  by  cavitation 
within  both  streamwise  and  spanwise  vortices,  as  seen  in 
Figs.  5(b) -5(e).  At  the  lowest  cavitation  number  [Fig.  5(e)] 
the  recirculating  region  beneath  the  cavity  becomes  filled 
with  a  bubbly  mixture.  These  observations  of  incipient  and 
developed  shear  layer  cavitation  are  consistent  with  those  of 
Katz  and  O’Hem^^  and  O’Hem.^^  Figure  6  presents  a  time 
series  of  the  developed  cavitating  shear  layer  formed  behind 
the  sharp-edged  plate  (cr=0.52,  I7oo  =  9.0m/s,  and  cq 
=  2.7  ppm,  where  cq  is  the  oxygen  content).  Such  high¬ 
speed  video  images  were  analyzed  to  determine  the  shear 
layer  growth  rate  and  vortex  spacing. 

O’ Hera  observed  that  the  growth  of  the  cavitating 
shear  layer  was  little  changed  from  that  of  the  noncavitating 
shear  layer,  and  a  similar  result  was  found  in  the  present 
study.  Figure  7  shows  a  plot  of  the  dimensionless  shear  layer 
thickness,  versus  the  dimensionless  downstream  dis¬ 

tance,  x/h.  The  growth  rate,  had  a  slope  of  0.25,  and 
this  compares  well  to  the  value  of  0.23  found  by  O’Hem^^ 
for  a  cavitating  shear  layer,  and  to  the  value  of  0.22  found  by 
Wygnanski  and  Fiedler. The  spacing  of  the  spanwise  cavi¬ 
tating  vortices  increases  with  an  increasing  downstream  dis¬ 
tance.  Figure  8  shows  the  average  spanwise  vortex  spacing, 
/,  as  a  function  of  the  downstream  distance,  x,  for  two  cavi¬ 
tation  numbers.  The  spanwise  vortex  spacing  constant  is  de¬ 
fined  as  //x,  and  this  was  found  to  be  0.47  for  cr=0.45  and 
0.51  for  a'=0.57.  A  value  of  0.50  was  found  by  O’Hem^^  for 
a  shear  flow  with  developed  cavitation.  Bemal  eta/.,** 
Heman  and  Jimenez, and  Jimenez^^  found  values  of  0.56 
for  noncavitating  turbulent  plane  shear  layer  with  Re^ 
>10^  The  values  found  for  the  cavitating  shear  layer  are 
equivalent  to  those  measured  for  the  noncavitating  layer, 
within  the  uncertainty  of  the  measurement. 

Figure  9  presents  the  average  spacing  and  diameter  of 
the  spanwise  vortices,  and  the  length  of  the  streamwise  vor¬ 
tices  as  a  function  of  cavitation  number.  These  quantities 
were  averaged  over  the  entire  field  of  view  (about  10  cm  in 


Downloaded  05  Aug  2003  to  141.212.134.149.  Redistribution  subject  to  AlP  license  or  copyright,  see  http://ojps.aip.org/phf/phfcr.jsp 


3420  Phys.  Fluids,  Vol.  14,  No.  10,  October  2002 


C.  O.  Iyer  and  S.  L  Ceccio 


FIG.  5.  (a)  Photograph  of  incipient  cavitation  for  (r=0.65.  f/»  =  9.0  m/s.  (b)  Photograph  of  incipient  cavitation  for  C/*  =  9.0  m/s.  (c)  Photograph  of 

incipient  cavitation  for  cr^O.SO,  f/oc  =  9.0  m/s.  (d)  Photograph  of  incipient  cavitation  for  (r=0.41,  f/„  =  9.0m/s.  (e)  Photograph  of  incipient  cavitation  for 
o-=0.33,  U^^9.0m/s. 


the  streamwise  direction  and  15  cm  in  the  spanwise  direc- 
tion)  for  the  growing  shear  layer.  The  lengths  in  this  figure 
were  nondimensionalized  by  2^,  where  ^=26.2  mm  is  half 
the  shear  layer  thickness  measured  in  region  A.  Only  small 
changes  are  seen  with  a  variation  in  the  cavitation  number.  In 
Fig,  9,  the  average  spacing  of  the  spanwise  vortices,  1,  de¬ 
creases  with  decreasing  cavitation  number  as  the  number  of 
vortices  observed  increases.  In  the  same  figure,  a  slight  in¬ 
crease  in  the  diameter  of  the  spanwise  vortices  and  the  length 
of  the  streamwise  vortices  is  observed,  although  these 
changes  are  well  within  the  range  of  uncertainty.  Note  that 
there  is  a  slight  increase  in  the  diameter  and  the  length  of 
these  vortices  for  cr=0.50  as  the  cavitation  number  is  low¬ 
ered  to  0.45. 


If  a  noncavitating  streamwise  vortex  is  simply  modeled 
as  a  Rankine  vortex  with  core  radius,  a,  the  pressure  within 
the  core,  ,  is  given  by 


where  is  the  strength  of  the  vortex.  At  the  point  of  cavi¬ 
tation  inception  in  the  vortex,  the  core  pressure  reaches  the 
vapor  pressure,  P^=Py,  assuming  that  the  water  is  suffi¬ 
ciently  populated  with  cavitatable  nuclei  (i.e.,  the  water  is 
“weak’')*  In  this  case,  the  vortex  strength  would  be  given  by 

Ts=27raU„\f^i, 
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FIG.  6.  A  time  series  of  images  of  shear-flow  cavitation  for  cr=0.52  and  U^-9.0  m/s.  The  flow  is  from  left  to  right  in  each  image.  The  series  progresses  from 
top  to  bottom  and  left  to  right.  The  time  interval  between  each  frame  is  MI  ms. 


where  cr,  is  the  cavitation  number  at  inception.  O' Hem 
examined  cinemagraphic  records  of  stream  wise  vortex  cavi¬ 
tation  inception,  and  the  size  of  the  initial  bubble  in  the  vor¬ 
tex  at  inception  was  determined.  It  was  then  assumed  that  the 
bubble  radius  was  approximately  equal  to  the  original  core 
radius,  and  this  measurement  was  used  to  determine  F  j .  The 
spanwise  vortex  strength,  F p ,  was  estimated  to  be  Fp 
where  At/  is  the  velocity  difference  across  the 
layer  and  1  is  the  local  spanwise  vortex  spacing.  O’ Hem 
concluded  that  the  ratio  Fp/F^  was  between  0.01  and  0.10 
near  the  inception  point  of  the  streamwise  vortices.  This  ratio 
was  smaller  than  the  value  suggested  by  Jimenez  et 


(who  reported  F5/Fp«=^0.6±0.8  for  a  single-phase  plane 
shear  layer).  Recent  investigations  into  developed  tip-vortex 
cavitation  suggest  that  incipient  vortex  cavitation  does  not 
fill  the  original  viscous  core  (Brian9on-Maijollet  and 
Merle^^).  As  the  vortex  fills  with  vapor,  the  diameter  of  the 
vapor  tube  may  be  several  times  the  size  of  the  original  core 
size,  but  this  is  principally  due  to  diffusion  of  gas  into  the 
vortex  core.  Arndt  and  Keller^®  analyzed  a  Rankine  vortex 
before  and  after  the  advent  of  developed  cavitation.  If  angu¬ 
lar  momentum  is  conserved,  the  vapor  core  diameter  will  be 
al^/2.  Examining  the  data  at  low  air  content  (o-j—l.l),  we 
will  assume  that  the  core  size  of  noncavitating  the  vortex 
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FIG.  7.  Growth  of  shear  layer  thickness,  as  a  function  of  downstream 
distance,  xfh^  for  cr— 0.50  and  t/*  =  9m/s  (b).  The  solid  line  is  the  linear 
regression  =  0.25(jc/A  +  0.02) .  These  data  are  compared  to  that  of 

O’Hem  (1990),  who  found  a  relationship  of  =  0.225(jf//i  +  0.05) 
shown  with  the  dashed  line. 

was  1.4  times  the  vapor  core  diameter.  The  average  value  of 
Tj  is  then  about  0.14  m^/s,  and  T/)  is  approximately  0.7  m^/s. 
Then,  r5/r,.'«0.2.  As  we  have  assumed  a  larger  initial  core 
radius,  this  value  is  somewhat  larger  than  that  of  O’Hem. 
However,  it  is  still  smaller  than  that  expected  from  examina¬ 
tion  of  noncavitating  mixing  layers. 

C.  Void  fraction  of  the  cavitating  shear  layer 

Figure  10  presents  a  plot  of  the  local  void  fraction  mea¬ 
sured  in  region  A.  Values  reported  are  for  two  oxygen  con¬ 
tents,  and  for  both  the  light  scattering  and  digital-imaging 
measurement  methods.  As  the  cavitation  number  decreases, 
the  void  fraction  of  the  shear  layer  increases,  as  expected. 
However,  the  light  scattering  method  consistently  measured 
lower  void  fractions  compared  to  the  digital  images  measure¬ 
ment  at  higher  cavitation  numbers.  The  calculation  of  the 
void  fraction  from  the  bubble  event  rate  employed  a  fixed 
bubble  residence  time,  of  250  yus,  in  the  measurement  vol¬ 


FIG.  8.  Growth  of  the  spanwise  vortex  spacing.  Vh,  as  a  function  of  down¬ 
stream  distance,  xfh,  for  7=  0.57  (•)  and  cr=0.45  (■)  and  11^  =  9  mis. 
The  solid  lines  are  linear  regressions  l//j -0.51  (x//2  + 0.66)  and  Mh 
=  0.47(.t//i  +0.85).  These  data  are  compared  to  that  of  O’Hem  (1990),  who 
found  a  relationship  of  I//1  =  0.50(a-//i +0.20)  shown  with  a  dashed  line. 


FIG.  9.  The  average  spacing  (■),  length  (A),  and  diameter  (•)  of  the 
streamwise  vortices  as  a  function  of  cavitation  number  observed  in  region  A. 
The  lengths  are  nondimensionalized  by  IS. 


ume.  However,  the  actual  residence  time  will  likely  increase 
with  decreasing  cavitation  number  as  the  average  length  of 
the  streamwise  bubbles  increases.  Moreover,  the  light  scat¬ 
tering  method  measured  the  void  fraction  at  a  small  volume 
within  the  center  of  the  shear  layer,  while  the  light-sheet 
method  averaged  the  void  fraction  over  the  entire  region  A. 
Yet,  even  with  these  differences,  the  void  fractions  are  quali¬ 
tatively  similar  with  maximum  average  void  fractions  of 
1%-1.5%  at  the  lowest  cavitation  numbers.  Note  that  the 
average  void  fraction  is  reduced  as  the  dissolved  oxygen 
content  is  reduced  for  a  given  cavitation  number.  A  reduction 
in  dissolved  gas  content  is  often  associated  with  a  reduction 
in  the  free  gas  content  (i.e.,  the  number  of  small  free-stream 
bubbles),  and  with  a  reduction  in  the  number  of  free-stream 
nuclei  the  extent  of  the  cavitation  will  be  reduced. 

D.  Mean  pressure 

The  mean  pressure  difference  measured  between  the  sur¬ 
face  taps  T1  and  T2,  shown  in  Fig.  1,  was  5.8±0.2  kPa.  This 


Cavitation  Number 


FIG.  10.  Void  fraction  in  region  A  as  a  function  of  cavitation  number.  Two 
methods  were  used  to  determine  the  void  fraction.  Light  scattering  was 
employed  at  a  single  point  in  the  center  of  region  A  for  3  ppm  (•)  and  1 
ppm  (A)  dissolved  oxygen  content.  The  void  fraction  was  inferred  from 
planar  digital  images  in  region  A  for  3  ppm  dissolved  oxygen  content  (□). 
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FIG.  11.  Average  flow  in  region  B  presenting  the  nondimensional  mean  velocity,  («- t/Mean.B)/^^Mean,B  the  nondimensional  vorticity  contours, 
0)j2SfUMem,B'  Thc  flow  is  from  left  to  right.  Here  1000  individual  images  were  averaged. 


value  was  not  influenced  by  the  extent  of  cavitation  within 
the  test  section  over  the  range  of  cavitation  numbers  exam¬ 
ined.  This  suggests  that  the  mean  flow  through  the  test  sec¬ 
tion  was  not  substantially  altered  by  the  development  of 
cavitation  in  the  shear  layer. 

E.  Comparison  between  the  noncavltating  and  the 
cavitating  mean  flows  in  region  B 

As  the  shear  layer  enters  the  diffuser,  the  vapor  pockets 
in  the  cavitating  vortices  collapse  leaving  only  small  gas 
bubbles  in  the  flow.  PIV  measurements  were  performed  in 
this  downstream  region  of  the  cavitating  flow  using  only 
seed  particles  as  flow  tracers.  Measurements  were  performed 
downstream  of  the  noncavitating  shear  flow  (or =2.3)  and 
downstream  of  three  cavitating  shear  flows,  with  cavitation 
numbers  of  0.55,  0.50,  and  0.45.  As  shown  in  Fig.  10,  the 
void  fraction  of  the  cavitating  shear  layer  increases  sharply 
between  or=0.50  and  0.45.  By  conducting  measurements  in 
region  B,  the  lasting  influences  of  the  growth  and  collapse  of 
the  cavitating  vortices  on  the  turbulent  shear  flow  are  exam¬ 
ined. 

One-thousand  PIV  images  were  acquired  in  region  B  and 
averaged  for  each  cavitation  number.  Figure  1 1  presents  the 


dimensionless  mean  in-plane  velocity,  (u 
“  ^^Mcan.B)/^Mean.B »  superimposed  on  the  dimensionless  out- 
of-plane  vorticity,  Wj2^/t/iviean.B  *  Th^  mean  velocity, 
f^Mean,B »  has  been  subtracted  from  each  image.  Some  differ¬ 
ences  are  observed  in  the  mean  velocity  field  of  the  noncavi¬ 
tating  flow  compared  to  the  cavitating  shear  flow.  Figure  12 
shows  the  average  velocity  profiles,  ufUf^,  for  three  cross- 
stream  sections  of  region  B:  x/h  =  20.54  in  the  left  half  (up¬ 
stream)  of  the  image,  x/h  =  20,96  in  the  middle  of  the  image, 
and  x/h-2\.3S  in  the  right  half  (downstream)  of  the  image. 
The  average  free-stream  velocity,  U fs ,  was  computed  by  av¬ 
eraging  25  velocity  vectors  from  the  top  left  comer  of  the 
average  velocity  field. 

At  the  lowest  cavitation  numbers,  the  gradient  of  the 
velocity  across  region  B  is  increased.  The  influence  of  the 
cavitation  decreases  with  increasing  downstream  distance. 
Convection  of  the  flow  across  region  B  represents  nearly  one 
turnover  of  the  spanwise  eddies.  Figures  13  and  14  present 
profiles  of  the  average  vorticity  and  strain  rate  for  the  same 
three  cross-streams.  These  data  are  normalized  by  2S,  the 
average  shear  layer  thickness  of  region  A,  since  for  the  pur¬ 
pose  of  scaling  2S  is  the  order  of  magnitude  of  the  shear 
layer  thickness  in  region  B.  With  Aw/C/f5~±0.02,  and  with 
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FIG.  13.  Average  vorticity  profiles,  as  a  function  of  cavitation 

number  for  three  cross-stream  sections  of  region  B:  (a)  xlh  -  20.54,  (b) 
jc/A  =  20.96,  (c)  =  21.38.  O  is  the  noncavitating  case,  A  is  cr= 0.55,  •  is 

<7=0.50,  and  A  is  o-=0.45. 


FIG.  12.  Average  velocity  profiles,  w/t/fj,  as  a  function  of  cavitation  num¬ 
ber  for  three  cross-stream  sections  of  region  B  derived  from  the  data  in  Fig. 
13.  (a)  xlh^lOM,  (b)  xlh^2Q.96,  (c)  xfh^2\3%.  The  symbols  O  is  for 
the  noncavitating  case,  A  is  <7=0.55,  •  is  cr=0.50,  and  A  is  <7=0.45.  Here 
1000  individual  PIV  images  were  averaged. 


an  uncertainty  in  the  position  of  A.x/2^~ 0.001,  the  uncer¬ 
tainty  in  the  nondimensionalized  vorticity  and  strain  rate  is 
approximately  ±0.10.  After  averaging,  the  uncertainty  in  the 
vorticity  and  strain  profiles  is  ±0.003.  Differences  of  less 
than  ±0.006  are  within  the  uncertainty  of  the  measurement. 
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FIG.  14.  Average  strain  rate  profiles,  as  a  function  of  cavitation 

number  for  three  cross-stream  sections  of  region  B:  (a)  =  20.54,  (b) 

x//i  =  20.96,  (c)j://i  =  21.38.  O  is  the  noncavitating  case,  A  is  (r=0.55,  ■  is 
cr=0.50,  and  A  is  cr=0.45. 


The  maximum  absolute  vorticity  levels  are  higher  in  the 
cavitating  flows  than  in  the  noncavitating  flow,  especially  for 
the  cavitation  numbers  cr=  0.50  and  cr=0,45.  Similar  results 
for  the  strain  rate  profiles  are  shown  in  Fig.  14.  The  higher 
vorticity  and  strain  rates  are  consistent  with  the  decrease  in 


the  shear  layer  thickness  at  lower  cavitation  numbers.  The 
case  of  (7=0.50  is  again  anomalous  with  the  highest  maxi¬ 
mum  vorticity  and  strain  rate  magnitudes.  A  possible  expla¬ 
nation  for  this  may  be  the  significant  increase  in  the  void 
fraction  beneath  the  shear  layer  at  the  lowest  cavitation  num¬ 
bers.  This  is  observed  in  Figs.  5(d)  and  5(e).  Increasing  con¬ 
centration  of  bubbles  beneath  the  shear  layer  reduces  the 
mean  density  of  this  recirculating  flow.  This,  in  turn,  would 
reduce  the  shear  stresses  across  the  layer  upstream  of  region 
B. 

F.  Comparison  between  the  turbulent  fluctuations  of 
the  noncavitating  and  the  cavitating  flows  in 
region  B 

Figures  15,  16,  and  17  present  profiles  of  the  mean  Rey¬ 
nolds  stresses  u*u'IU\,  v^v*fU\,  and  —u*v^iu\  at  the 
three  cross-stream  locations.  The  average  velocity  at  each 
vector  location  was  determined,  and  this  local  average  veloc¬ 
ity  was  subtracted  from  each  vector  to  determine  the  fluctua¬ 
tion  velocity  components  u*  and  v'  for  each  image.  This 
data  was  then  used  to  compute  the  average  Reynolds 
stresses.  Figure  18  shows  the  convergence  plots  for  the  Rey¬ 
nolds  stresses  at  a  point  in  the  middle  of  the  domain.  It 
indicates  that  1000  images  are  sufficient  to  achieve  a  statis¬ 
tically  steady  average.  _ 

Consider  the  trends  in  u'u'fU\  in  Fig.  15.  First,  it 
should  be  noted  that  the  location  of  maximum  fluctuations  is 
not  in  the  center  of  the  shear  layer.  By  the  time  the  flow  is  in 
region  B,  the  shear  layer  is  in  the  diffuser  and  is  no  longer 
similar  to  the  planar  mixing  layer  of  region  A.  The  u'u'/Uj^ 
profiles  are  not  significantly  changed  with  decreasing  cavita¬ 
tion  number,  and  the  magnitude  of  the  fluctuations  increases 
with  the  extent  of  cavitation  upstream,  Laberteaux  and 
Ceccio^  showed  that  the  collapse  of  cavitation  bubbles  leads 
to  the  production  of  turbulence  within  the  flow,  and  this  in¬ 
crease  in  streamwise  turbulence  may  result  from  this  mecha¬ 
nism.  Again,  the  case  of  C7=0.50  is  anomalous,  as  the  level 
of  turbulence  observed  at  the  lower  boundary  of  the  shear 
layer  is  reduced  for  just  this  case. 

For  the  measurements  of  v'v'/U\  (Fig,  16),  the  pres¬ 
ence  of  cavitation  has  reduced  the  extent  of  the  maximum 
cross-stream  fluctuations  compared  with  the  noncavitating 
flow  by  as  much  as  30%.  These  results  suggest  that  the  pres¬ 
ence  of  cavitation  either  damps  out  the  cross-stream  fluctua¬ 
tions  or  reduces  their  production.  Similar  trends  are  found  in 
the  measurements  of  -u'v'/u\  (Fig.  17),  where  an  increase 
in  the  extent  of  cavitation  reduced  the  correlation  between 
the  streamwise  and  cross-stream  velocities.  The  streamwise 
stretching  of  a  vortex  will  lead  to  an  increase  in  the  rotation 
rate  of  the  vortex;  thus,  streamwise  velocity  fluctuations  are 
coupled  to  cross-stream  fluctuations.  Belahadji  et  al}^  sug¬ 
gested  that  the  presence  of  vapor  in  the  cores  of  such  vortices 
modifies  the  vortex-stretching  process  by  decoupling  the  rate 
of  vortex  straining  from  the  rotation  rate.  As  a  noncavitating 
vortex  is  stretched,  the  conservation  of  angular  momentum 
implies  that  the  rotation  rate  will  increase  and  the  pressure  in 
the  viscous  core  will  decrease.  However,  the  core  of  a  cavi¬ 
tating  vortex  can  maintain  a  constant  pressure,  and  stretching 
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FIG.  15.  Mean  Reynolds  stress  profiles,  u' u’ HJ\,  as  a  function  of  cavita¬ 
tion  number  for  three  cross-stream  sections  of  region  B:  (a)  jr//i==  20.54,  (b) 
xlh  =  20.96,  (c)  xlh  -  21.38.  O  is  the  noncavitating  case,  A  is  or=  0.55,  •  is 
cr=0.50,  and  A  is  o-=0.45. 


FIG.  16.  Mean  Reynolds  stress  profiles,  u'u'/f/J,,  as  a  function  of  cavita¬ 
tion  number  for  three  cross-stream  sections  of  region  B:  (a)  jr//i- 20.54,  (b) 
;c//i  =  20.96,  (c)  xlh  =  2X3%.  O  is  the  non-cavitating  case,  A  is  cr=0.55,  • 
is  O' -0.50,  and  A  is  cr=0.45. 


will  result  only  in  the  production  of  more  core  vapor  with 
little  change  in  the  vortex  core  diameter.  Such  a  mechanism 
could  be  responsible  for  a  decrease  in  the  production  of 
cross-stream  fluctuations  and  a  reduction  in  the  stress 
“w'uVf/l,  as  seen  in  the  data.  However,  these  distributions 


also  exhibit  some  anomalies.  Once  more,  the  case  of  <t 
=  0.50  shows  the  largest  reductions  in  average  Reynolds 
stresses  compared  to  the  noncavitating  case.  And,  the  case  of 
(7=0.45  shows  increasing  values  of  /u\  near  the  top 
edge  of  the  shear  layer. 
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FIG.  17.  Mean  Reynolds  stress  profiles,  -u'v' IU\y  as  a  function  of  cavi¬ 
tation  number  for  three  cross-stream  sections  of  region  B:  (a)  =  20.54, 

(b)  20.96,  (c)  ;c//i  =  21.38.  O  is  the  non-cavitating  case,  A  is  cr 

=  0.55,  •  is  cr=0.50,  and  ▲  is  o-=0.45. 


The  PIV  images  were  analyzed  to  reveal  the  probability 
distributions  of  the  instantaneous  vorticity  and  strain  rates 
measured  within  region  B.  Figures  19  and  20  show  the  prob¬ 
ability  distribution  for  the  values  of  vorticity  and  strain  rate 
taken  from  1000  images.  Differences  are  only  observed  at 


FIG.  18.  Convergence  plot  for  the  mean  Reynolds  stresses  at  a  point  in  the 
middle  of  region  B:  u'u' IU\  (solid  line);  v*v’IU\  (short-dashed  line);  and 
-u'v'lU\  (long-dashed  line). 


the  extreme  positive  and  negative  values,  and  the  trends  with 
cavitation  number  observed  in  the  profiles  can  be  detected  in 
variations  of  the  histograms.  For  the  case  of  cr=0.50,  there  is 
a  slight  shift  toward  a  higher  probability  of  strong  negative 
vorticity  and  strong  positive  strain  rates,  and  these  trends 
reverse  for  the  condition  of  or =0.45. 


FIG.  19.  Probability  density  function  of  the  normalized  vorticity,  Wj25/£/fs , 
measured  in  region  B  for  varying  cavitation  numbers.  The  measurements 
from  1000  PIV  images  are  presented.  O  is  the  noncavitating  case,  A  is  o* 
=  0.55,  •  is  cr=0.50,  and  A  is  (7=0.45. 
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FIG.  20.  Probability  density  function  of  the  normalized  strain  rate, 
SgySIU^i,  measured  in  region  B  for  varying  cavitation  numbers.  The  mea¬ 
surements  from  1000  PIV  images  are  presented.  O  is  the  noncavitating  case, 
A  is  (7—0.55,  •  is  (7=0.50,  and  A  is  (7=0.45. 


Q.  Influence  of  cavitation  on  the  shear  stress  across 
the  shear  layer 

The  data  presented  above  indicate  that  the  shear  layer 
was  modified  by  the  cavitation,  but  only  subtly.  As  we  are 
examining  the  flow  after  the  collapse  of  the  cavitation 
bubbles,  interpretation  of  the  results  can  be  problematic,  but 
it  appears  that  the  cavitation  resulted  in  an  increase  of  the 
mean  shear  rate  across  the  layer  and  a  decrease  in  the  Rey¬ 
nolds  stresses  within  the  layer.  While  the  flow  examined  here 
is  not  a  classical  mixing  layer,  the  scaling  of  a  variable  den¬ 
sity  mixing  layer  may  be  useful  in  examining  this  flow  (see 
Hermanson  and  Dimotakis^*  for  a  complete  derivation).  A 
similarity  variable  can  be  defined  as 

r)=ylS{x), 

where  x  is  the  streamwise  direction  and  y  is  the  cross-stream 
direction.  The  velocity  distribution  across  the  shear  layer 
takes  the  form 


U{  v)  =  sU 


dfiv) 

dr] 


and 


V{  77)  =  5f/i 


V' 


df{ri) 

dr] 


where  U  is  the  mean  streamwise  velocity,  pV=  pV~^ p'v' j 
and  U 1  is  the  velocity  of  the  high-speed  stream  on  the  mix¬ 


ing  layer.  The  density  ratio  s(r])  =  p/ps  compares  the  free- 
stream  density  to  the  local  shear  layer  density.  Substitution 
into  the  streamwise  momentum  equation  yields 

u]  dS  J  ds  J^f\_  1  dr  7)  dP  dS 

'KIP  dx^\drj^^^T^j  "pKlP  drj  JKIF  di]  dx’ 

where  T=-pu'v'  is  the  Reynolds  shear  stress  plus  higher- 
order  terms  and  AU  is  the  velocity  difference  across  the 
shear  layer.  Here,  it  is  assumed  that  streamwise  velocity  fluc¬ 
tuation  correlations  are  small  compared  to  and  that  the 
Reynolds  stresses  are  much  larger  than  the  viscous  stresses. 
The  pressure  gradient  term  on  the  right-hand  side  is  small 
compared  to  the  shear-stress  term.  The  second  term  on  the 
left  results  from  density  gradients  in  the  shear  layer.  The 
effect  of  buoyancy  has  been  neglected  here  as  the  flow  is 
momentum  dominated.  The  Richardson  number  is 


n- 

R>’“-77r-' 


where  g  is  the  gravitational  constant. 

While  the  advent  of  cavitation  alters  the  local  density  of 
the  shear  layer,  the  maximum  reduction  in  the  mean  density 
within  the  shear  layer  is  less  than  1%,  based  on  the  void 
fraction  measurements.  This  is  in  contrast  to  reacting  shear 
layers,  where  the  density  reduction  variation  can  be  as  high 
as  60%.  We  will  assume  here  that  the  gradient  of  density 
across  the  shear  layer  is  small  and  can  be  neglected.  The 
streamwise  momentum  equation  becomes 

^  u]  dS  d^f  ^  1  dr 

~^'KlFlxl^^‘‘'pKlF  dr)' 

This  equation  can  be  integrated  numerically  for  a  given  mea¬ 
sured  velocity  distribution  with  the  constraint  that  the  shear 
stresses  vanish  in  the  two  free-streams.  A  constraint  based  on 
mass  conservation  can  be  applied  to  achieve  this  result  (Her¬ 
manson  and  Dimotakis^^). 

Figure  21  shows  the  velocity  profiles  measured  in  region 
B  [from  Fig.  12(a)]  scaled  with  the  measured  velocity  differ¬ 
ence  across  the  layer.  The  profiles  are  nearly  identical,  and 
they  are  approximated  by  the  function  tanh[2(y/(5-l)],  and 
this  is  approximately  the  form  of 

dfi  v)  AU 

The  function  /(rj)  can  be  found  by  integrating  df{r])ldr]. 
The  maximum  value  of  the  shear  stress  is  found  on  the  di¬ 
viding  streamline,  where  V—rf^U,  Let  us  assume  that  the 
profiles  across  the  layer  are  symmetric  in  {rj-  7/0)=  Both 
the  shear-stress  distribution  r(  fj)  and  d^f{  f])ld  ff'  are  even 
functions  of  f].  Consequently, /(?))  is  an  odd  function  of  f]. 
An  approximate  expression  for  f{f])  is 

AU 

/( 

Also,  the  derivative  of  df{f])fdf]  yields 


Downloaded  05  Aug  2003  to  141.212.134.149.  Redistribution  subject  to  AlP  license  or  copyright,  see  http://ojps.aip.org/phf/phfcr.lsp 


Phys.  Fluids,  Vol.  14.  No.  10,  October  2002 


The  influence  of  developed  cavitation  3429 


FIG.  21.  Scaled  velocity  profiles,  ^  ^  funciioa  of  cavita¬ 

tion  number  for  three  cross-stream  sections  of  region  B  for  xlh  ~  20.54 
derived  from  the  data  in  Fig.  14(a):  O  is  the  noncavitating  case,  A  is  tr 
=  0.55,  •  is  cr=0.50,  and  ▲  is  o-=0.45.  Also  shown  is  the  velocity  profile 
tanh[2(y/<5-l)]. 


d  rf" 


AU 

—  sech^(2^). 


These  expressions  can  be  integrated  to  determine  an  expres¬ 
sion  for  the  shear  stress  across  the  shear  layer: 


d^f  AU 

—  [ 2  7/  tanh( 2  ^)  -  In  cosh(2  v)-C], 


manner  similar  to  a  rearward  facing  step,  where  AU^Uf^, 
With  this  assumption,  (- w' u'/C/fs) Max 0-02,  which  is  close 
to  the  values  measured  here. 

The  above  relationship  can  be  used  to  determine  trends 
in  the  shear  stress  with  changes  in  the  extent  of  cavitation 
within  the  shear  layer.  As  the  cavitation  number  decreases, 
the  mean  density  in  the  cavitating  shear  layer  decreases  by 
less  than  2%,  and  no  significant  variation  was  observed  in 
the  growth  rate,  although  there  is  considerable  uncertainty  in 
this  measurement.  The  mean  profiles  in  the  region  just  down¬ 
stream  of  the  cavitation  suggest  that  AU  in  region  B  in¬ 
creases  by  about  10%  at  the  lowest  cavitation  number.  A 
reduction  in  mean  density  would  suggest  a  reduction  in  the 
shear  stresses  in  the  cavitating  layer,  while  an  increase  in  the 
strain  rate  within  the  layer  would  suggest  an  increase.  An 
examination  of  the  Reynolds  stresses  in  Fig.  17  indicates  that 
the  advent  of  cavitation  reduced  the  stresses  by  about  a  factor 
of  2,  and  this  is  accompanied  by  a  reduction  of  cross- stream 
fluctuations  shown  in  Fig.  16.  The  changes  in  mean  density 
or  growth  rate  do  not  account  for  such  a  difference,  and  the 
changes  in  the  strain  rate  across  the  layer  suggest  an  increase 
in  shear  stress  with  cavitation. 

Reduction  in  the  Reynolds  stresses  is  most  pronounced 
for  the  case  of  o'=0.50.  This  is  the  condition  when  the  cavi¬ 
tation  has  developed  the  most  without  the  formation  of  a 
large  bubbly  mixture  beneath  the  cavity.  For  the  case  of  cr 
=  0.45,  the  flow  beneath  the  shear  layer  is  full  of  large  vapor 
bubbles  that  are  trapped  and  recirculated  [see  Figs.  5(d)  and 
5(e),  for  example].  It  is  possible  that  the  extent  of  the  shear 
across  the  shear  layer  is  reduced  by  the  presence  of  this 
vapor.  This  could  explain  why  the  trend  reverses  with  a  fur¬ 
ther  decrease  in  cavitation  number  below  0.50. 


where  C  is  a  constant  chosen  such  that  the  shear  stress  tend 
to  zero  for  large  magnitudes  of  In  the  limit  as  ?)— >co, 
2  7}  tanh(2^)-ln  cosh(2^)— >ln(2)= 0.693 The  shear 
stress  distribution  is  then  given  by 

\  dS  ^  T 

-  4  [2  ^  tanh(2  -  In  cosh(2  v) . 

The  maximum  shear  stress  at  ^=0  is  then 

1  dS  Tniax 

^4  dx  pKiP' 

This  scaling  indicates  that  the  magnitude  of  the  maximum 
shear  stress  decreases  linearly  with  a  decreasing  shear  layer 
growth  rate  and  decreasing  mean  flow  density. 

Plane  mixing  layers  have  growth  rates  of  dS/dx^0A6y 
which  suggests  that  Tj^^/pAU^^ 0,027.  This  estimate  is 
nearly  twice  the  peak  measured  Reynolds  stresses  r 
=  -  u^v'/AU^^0A3  (Spencer  and  Jones^^).  With  dSfdx 
•=^0.25,  pAU^^O.OA.  This  suggests  that  the  maximum 
value  of  the  Reynolds  stress  (-w'y 7t7|)Max‘^0.02* (Aty/ 
(7fj.)^,  If  AU  is  the  measured  velocity  difference  across  the 
layer,  (~M'i;7t/fs)Max*^ 0.002,  and  this  is  about  an  order  of 
magnitude  smaller  than  the  measured  value  of  between  0.02 
and  0.03.  However,  in  the  present  experiment,  the  high-speed 
flow  essentially  drives  a  recirculating  low  speed  flow  in  the 


IV.  CONCLUSIONS 

Visual  observations  of  the  shear  layer  suggest  that  the 
overall  formation,  growth,  and  convection  of  the  primary  and 
secondary  vortical  structures  in  the  shear  layer  are  not  sig¬ 
nificantly  affected  by  the  presence  of  cavitation.  The  average 
flow  field  downstream  of  the  cavitating  shear  layer  is  not 
significantly  altered  both  in  terms  of  average  velocities  and 
mean  pressure  drop  across  the  test  section. 

The  largest  differences  between  the  cavitation  and  the 
noncavitating  flows  were  in  the  turbulent  fluctuations  exam¬ 
ined  downstream  of  the  cavitating  shear  layer.  The  stream- 
wise  fluctuations  increased  by  about  15%  compared  to  the 
noncavitating  case  in  the  center  of  the  shear  flow.  However, 
the  maximum  cross-stream  fluctuations  and  Reynolds 
stresses  decreased  by  about  30%.  As  bubbles  collapse,  an 
increase  in  turbulence  levels  is  expected  (Laberteaux  and 
Ceccio^’^).  The  decrease  in  the  cross-stream  fluctuations  and 
Reynolds  stresses  suggest  that  the  cavitation  within  the  cores 
of  streamwise  vortices  has  decreased  the  coupling  between 
the  streamwise  and  cross-stream  velocity  fluctuations.  This  is 
consistent  with  the  hypothesis  of  Belahadji  et  al,^^  whereby 
cavitation  in  the  cores  of  the  streamwise  vortices  decoupled 
the  stretching  and  rotation  rate  of  these  flow  structures.  The 
reduction  in  Reynolds  stresses  downstream  of  the  shear  layer 
is  inconsistent  with  the  scaling  of  the  shear  stress  within  the 
shear  layer  based  on  the  shear  layer  growth  rate,  mean  den- 
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sity,  and  the  mean  flow  profiles.  As  these  measurements  have 
limited  spatial  resolution,  it  may  be  necessary  to  examine 
turbulent  fluctuation  at  smaller  scales  to  resolve  this  issue. 

The  production  and  distribution  of  turbulent  fluctuations 
seem  to  have  been  significantly  altered  by  the  presence  of 
developed  cavitation.  It  is  problematic  to  interpret  the  data 
from  this  experiment  in  the  context  of  a  true  mixing  layer. 
Moreover,  we  are  left  to  infer  the  flow  processes  within  the 
cavitating  layer  from  observations  of  the  noncavitating  flow 
downstream.  It  may  be  necessary  to  examine  the  turbulent 
cavitating  flow  directly  to  resolve  the  questions  regarding  the 
shear  stress  distribution  and  Reynolds  stresses.  However,  the 
mean  shear  and  Reynolds  stress  distribution  in  turbulent 
flows  can  be  substantially  altered  when  the  rheological  be- 
havior  of  the  flow  is  altered.  An  example  is  the  addition  of 
high  surfactants  to  a  turbulent  shear  flow.  Warholic  et  al?^ 
demonstrated  that  the  addition  of  relatively  small  quantities 
of  surfactants  resulted  in  the  substantial  suppression  of  cross¬ 
stream  velocity  fluctuations  and  Reynolds  stresses  in  a  tur¬ 
bulent  channel  flow.  This  resulted  in  a  substantial  reduction 
in  the  shear  stress  at  the  wall  of  the  channel.  Similarly,  re¬ 
searchers  have  noted  that  the  addition  of  microbubbles  will 
lead  to  the  reduction  of  wall  shear  stresses,  and  dynamic 
shear  stress  measurements  at  the  wall  suggest  that  the 
bubbles  are  also  altering  the  dynamics  of  the  turbulent  flow 
(Merkle  and  Deutsch^"*).  It  is  therefore  possible  that  the  pres¬ 
ence  of  the  cavitation  within  the  shear  layer  modifying  the 
distribution  of  mean  shear  stress  across  the  layer  through  a 
change  in  the  effective  rheology  of  the  flow. 

Cavitating  wakes  have  been  significantly  modified  by 
developed  cavitation,  especially  with  respect  to  the  rate  of 
vortex  shedding.  Such  effects  were  not  observed  here.  The 
void  fraction  of  the  cavitating  shear  layer  was,  at  most,  2%. 
It  is  unclear  if  the  void  fraction  in  the  separated  region  be¬ 
hind  the  bluff  bodies  is  significantly  higher,  but  this  may  be 
the  case.  Moreover,  the  dynamics  of  the  shedding  process 
behind  bluff  object  is  altered  by  changes  in  the  pressure  field 
near  the  bluff  body.  The  presence  of  developed  cavitation 
can  alter  the  overall  pressure  field,  especially  if  the  cavitating 
object  is  confined  in  a  flow  channel.  In  the  present  study,  the 
presence  of  the  cavitation  did  not  significantly  change  the 
average  pressure  field. 

Last,  these  results  suggest  that  a  single-phase  simulation 
of  shear  flow  may  be  used  to  determine  the  extent  of  the 
cavitation  within  the  shear  flow.  A  method  is  needed  that  can 
successfully  simulate  the  strongest  vortices  in  the  flow  over  a 
range  of  scales.  While  such  single-phase  simulations  are,  at 
present,  difficult  to  conduct  at  high  Reynolds  numbers,  meth¬ 
ods  such  as  large  eddy  simulations  may  eventually  be  used  to 
reveal  the  largest  vortical  structure  in  a  flow  field.  Once  iden¬ 
tified,  the  extent  of  the  cavitation  can  be  estimated.  This  is 
particularly  true  of  incipient  cavitation  where  the  influence 
on  the  underlying  shear  flow  may  be  neglected. 
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Abstract 

The  tip-leakage  vortex  occurring  on  a  ducted  rotor  was  examined  using  both  three  component  Laser  Doppler 
Velocimetry  (LDV)  and  planar  Particle  Imaging  Velocimetry  (PIV).  The  vortex  strength  and  core  size  were 
examined  for  different  vortex  cross  sections  downstream  of  the  blade  trailing  edge.  The  variability  of  these 
quantities  are  observed  with  PIV  and  the  average  quantities  are  compared  between  LDV  and  PIV.  Developed 
cavitation  is  also  examined  for  the  leakage  vortex.  The  implication  of  vortex  variability  on  cavitation 
inception  is  discussed. 


1  Introduction 


Tip- vortex  cavitation  is  a  source  of  noise,  propeller-induced  hull  vibration  and  erosion  of  downstream  rudders. 
With  a  ducted  propeller,  the  tip  clearance  vortex  can  cause  damage  to  the  duct  surface  as  well.  The  pressure 
difference  across  the  two  sides  of  the  blade  along  with  the  presence  of  a  gap  between  the  duct  and  propeller 
blade  result  in  a  tip  clearance  or  leakage  flow.  The  leakage  flow  interacts  with  the  through-flow  across  the 
suction  side  of  the  blade  to  form  a  sheet  of  vorticity  that  rolls  up  into  a  tip  clearance  vortex.  Tip)- vortex 
cavitation  occurs  when  the  pressure  in  the  core  of  this  vortex  drops  below  vapor  pressure  and  small  bubbles 
within  the  vortex  core  begin  to  explosively  grow.  With  further  decreases  in  pressure,  the  core  of  the  vortex 


can  fill  with  vapor. 

Researchers  have  examined  the  flow  of  tip  and  leakage  vortices  in  order  to  determine  when  cavitation  will 
occur  and  to  develop  predictive  scaling  laws.  Recent  studies  of  tip- vortex  cavitation  are  reviewed  by  Fruman 
et  a/.(1992),  Fruman  et  a/.(1995),  and  Maines  and  Arndt  (1997).  Research  on  leakage-vortex  cavitation  is 
more  limited,  and  a  recent  review  is  given  by  Boulon  et  al.  (1999).  Of  particular  note  is  the  study  of  Farrell 
and  Billet  (1994)  where  the  scaling  of  leakage  vortex  inception  was  presented. 

Scaling  of  vortex  cavitation  inception  often  follows  the  idea  of  McCormick  (1962)  whereby  the  pressure 
within  the  core  of  the  vortex  can  be  related  to  the  viscous  flow  over  the  lifting  surface  and  the  lift  that  the 
surface  produces.  If  axisymmetric  vortex  flow  is  assumed,  the  local  depression  of  pressure  at  the  center  of 
the  vortex  is  proportional  to  the  circulation  of  the  vortex  and  the  size  of  the  viscous  vortex  core.  McCormick 
assumed  that  the  core  size  is  proportional  to  the  mean  boundary  layer  thickness  on  the  lifting  surface,  and 
that  the  vortex  circulation  was  proportional  to  the  lift  coefficient,  Cl>  Reynolds  number  dependence  arises 
from  the  assumption  that  the  boundary  layer  thickness  grows  as  Lastly,  it  is  assumed  that  the  vortex 

will  cavitate  when  the  core  pressure  falls  below  vapor  pressure.  When  combined,  these  relations  lead  to  the 
scaling: 

a  «  ClRe^^^.  (1) 


This  scaling  has  been  successfully  used  to  determining  the  point  of  cavitation  desinence  across  Reynolds 
numbers.  Scaling  of  cavitation  inception  must  also  include  the  effect  of  both  the  dissolved  and  free  gas  content 
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in  the  flow.  Farrell  and  Billet  (1994)  discuss  how  McCormick  scaling  should  be  modified  for  application  to 
tip-leakage  flows,  especially  with  regard  to  the  choice  of  scaling  for  the  core  size  with  varying  tip  clearance. 

While  the  above  scaling  has  proven  useful  and  predictive,  there  are  features  of  tip- vortex  flows  that  are 
not  directly  captured.  These  include  the  dynamics  of  the  vortex  during  roll-up  near  the  point  of  vortex 
formation,  the  effect  of  secondary  vorticity  in  the  roll-up  region,  the  influence  of  axial  flow  within  the  vortex 
core,  and  the  effect  of  turbulence  near  and  within  the  core  (Pauchet,  1997).  For  tip-leakage  flows,  there  is 
the  added  influence  flow  near  the  shroud. 

We  wish  to  examine  the  process  of  tip-leakage  vortex  cavitation  inception  and  its  scaling  as  it  occurs 
on  a  ducted  rotor.  Our  overall  goal  will  be  to  determine  what  physical  processes  are  at  work  within  the 
flow  field  to  lead  to  the  inception  of  discreet  cavitation  events.  As  in  previous  studies  of  vortex  cavitation 
inception,  we  will  examine  the  flow  field  using  three-component  Laser  Doppler  Velocimetry  (LDV).  LDV  pro¬ 
vides  time-averaged  measurements  of  the  three-dimensional  flow-field  with  relatively  high  spatial  resolution. 
However,  it  is  well  known  that  global  unsteadiness  of  the  vortex  flow  field  (“vortex  wandering”)  can  lead 
to  erroneous  interpretation  of  average  flow  fields,  especially  with  regard  to  the  estimation  of  turbulent  flow 
quantities  (Straka  and  Farrell,  1992).  Consequently,  we  have  combined  the  LDV  measurements  with  inst^- 
taneous  velocity  field  measurements  acquired  with  planar  Particle  Imaging  Velocimetry  (PIV)  to  examine 
the  variability  in  the  vortical  flow.  Tip-vortex  cavitation  and  near-tip  velocity  distributions  were  measured 
on  a  three-bladed  propeller  designated  “Propeller  5206”.  The  tests  were  performed  in  the  NSWCCD  36-inch 
water  tunnel.  Presented  here  are  some  preliminary  results  of  this  study  and  an  accompanying  discussion  of 
their  implications. 


2  Experimental  Setup 


Water  Tunnel  All  measurements  were  made  in  the  David  Taylor  36-inch  Variable  Pressure  Water  Tunnel. 
The  tunnel  is  a  recirculating  design  with  interchangeable  test  sections.  (See  Brownell,  1962,  for  a  detailed 
description  of  the  water  tunnel).  The  36-inch  diameter,  open  jet  test  section  was  used  for  these  tests.  For 
these  test,  the  propeller  was  driven  using  the  upstream  d3mamometer.  Inflow  to  the  propeller  was  uniform 
except  for  the  wakes  from  three  upstream  shaft  support  struts. 

Propeller  Model  Propeller  5206  is  a  three  bladed  rotor  designed  to  operate  in  a  cylindrical  duct.  The  duct 
in  this  case  is  a  cylindrical  extension  of  the  NSWCCD  36-inch  water  tunnel  flow  nozzle.  This  configuration 
produces  an  inner  duct  diameter  of  0.8636  m  (34  inches),  the  largest  propeller  operated  in  the  36-inch  water 
tunnel.  The  rotor  can  be  operated  at  speeds  up  to  500  rpm  resulting  in  blade  Reynolds  Numbers  on  the 
order  of  1/4  of  full  scale. 


Figure  1:  Propeller  5206  in  the  36-inch  water  tunnel  at  the  Naval  Surface  Warfare  Center,  Carderock  Division 
and  a  view  of  the  propeller  along  the  shaft  axis 

The  rotor  has  a  diameter  of  0.8504  m  (33.475  inches),  maintaining  a  0.015ilp  (Radius  Rp  =  0.4252  m) 
tip  clearance  resulting  in  a  6.67  mm  (0.2625  in)  tip  gap.  The  blade  chord-length  is  a  constant  0.381  in  (15 
in)  across  the  span.  Blade  thickness  approaches  10%  of  chord  at  the  tip.  The  tip  geometry  is  simple  with  a 
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constant  3.175  mm  (1/8  in)  radius  wrapped  about  the  tip.  The  ratio  of  the  tip  clearance  to  the  maximum 
thickness  at  the  tip  is  0.17,  which  is  near  the  optimum  value  of  0.2  reported  by  Farrell  and  Billet  (1994). 
The  propeller  was  made  using  a  numerically  controlled  milling  process  and  manufactured  to  tolerances  of 
approximately  0.1  mm  (0.004  in)  from  a  monobloc  6061T6  aluminum  forging.  The  blade  section  is  a  NACA 
66,  DTMB  (David  Taylor  Model  Basin)  modified  thickness  form,  with  a  =  0.8  meanline  camber,  with  trailing 
edges  thickened  to  incorporate  a  typical  anti-singing  trailing  edge  bevel.  The  blade  pitch  was  reduced  at  the 
hub  and  tip  to  minimize  root  and  tip  cavitation.  Table  1  summarizes  the  blade  geometry.  Figure  2  shows 
the  test  setup  with  the  propeller,  duct,  LDV  and  PIV  systems  as  well  as  the  trajectory  of  the  developed 
cavitating  vortex. 


Table  1.  Propeller  5206  Blade  geometry  description. 


rIRp 

C/D 

P/D 

It/D 

T/C 

F/C 

Os,  Degrees 

0.416 

0.446 

0.92 

-0.0063 

0.172 

-0.065 

-0.4 

0.5 

0.446 

1.135 

0.0063 

0.148 

0.0 

4.2 

0.6 

0.446 

1.22 

0.0145 

0.119 

0.037 

10.0 

0.7 

0.446 

1.175 

0.0227 

0.102 

0.047 

15.90 

0.8 

0.446 

1.087 

0.0309 

0.100 

0.047 

21.5 

0.9 

0.446 

0.995 

0.0391 

0.100 

0.038 

26.3 

1.0 

0.446 

0.905 

0.0473 

0.100 

0.02 

30.0 

Rp,  Tip  radius 
C,  Chordlength 
P,  Pitch 

Jr,  total  rake,  axial  displacement  of  blade  sections  relative  to  propeller  centerline 
r.  Blade  section  thickness 
P,  Blade  section  camber 

9s,  Blade  section  skew,  angular  displacement  of  blade  sections  relative  to  blade  line  perpendicular  to  the  hub 


Figure  2:  A  drawing  of  the  propeller  with  the  duct,  LDV  system,  PIV  system,  and  vortex  trajectory.  Note 
the  LDV  laser  beams  and  the  PIV  laser  sheet  and  camera.  The  camera  is  looking  up  perpendicular  to  the 
laser  sheet.  The  camera  mount  is  on  the  same  traverse  as  the  LDV  lasers. 

Test  Conditions  Operation  of  the  propeller  poses  interesting  challenges  due  to  its  large  size  and  its  complete 
blockage  of  the  test  section  flow.  Under  initial  operation  of  the  propeller,  attempts  to  run  at  its  design  advance 
coefficient  of  J  =  0.75  required  completely  shutting  down  the  water  tunnel  impeller.  At  around  500  rpm, 
the  rotor  mass  flow  was  such  that  the  impeller  would  begin  to  free  wheel,  reducing  the  flow  restriction  of 
the  stationary  impeller  blades  and  thus  dropping  the  measured  rotor  thrust  and  torque.  A  unique  feature 
of  the  36-inch  Water  Tunnel  is  the  controllable  pitch  impeller  blade,  which  permits  impeller  operation  over 
various  speeds  and  blade  pitch.  It  was  found  that  by  reducing  the  pitch  of  the  impeller  blades  to  around  zero, 
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the  tunnel  flow  could  be  controlled,  enabling  the  rotor  to  operate  over  a  relatively  large  range  of  advance 
coefficients,  0.65  <  J  <  1.0.  The  standard  test  conditions  at  500  rpm,  the  highest  speed  tested,  was  as 
follows  in  the  table  below.  The  inlet  pressure  was  maintained  high  enough  (45  psia  for  500  rpm)  during  the 
LDV  and  PIV  measurements  to  suppress  cavitation. 


Table  2.  Measured  flow  conditions. 


RPM 

Advance 
Coefficient,  J 

Velocity 
ft/s  (m/s) 

Thrust 
lbs  (kN) 

Torque 
ft-lbs  (N-m) 

Kt 

Kq 

500 

0.972 

22.85  (6.89) 

2560  (11.39) 

1280  (1735) 

0.314 

0.056 

7  X  10« 

The  operating  advance  coefficient  was  selected  to  produce  a  typical  leakage  vortex,  without  unwanted 
flow  conditions  over  the  blade.  The  operating  J  =  0.971  was  increased  relative  to  design  J  =  0.75  to  avoid 
the  formation  of  the  leakage  vortex  at  the  leading  edge.  Additional  increase  in  J  was  avoided  to  prevent  the 
occurrence  of  pressure  side  leading  edge  cavitation  near  the  blade  tip. 

High-Speed  Video  Camera  Photos  of  the  cavitating  vortex  were  taken  with  a  Broadcast  quality,  Sony 
Digital  Betacam  camera.  Examples  of  the  cavitating  vortex  are  shown  in  Figures  12  and  13. 

Laser  Doppler  Velocimetry  An  all  fiber-optic  LDV  system  allowed  for  three  component  velocity  measure¬ 
ments.  The  measurement  volume  was  positioned  at  a  point  in  the  horizontal  plane  containing  the  propeller 
axis.  The  probe  volumes  are  ellipsoidal,  0.7  x  0.7  x  1.3  mm.  The  probe  volumes  for  the  axial  and  tangential 
components  are  coincident,  and  are  oriented  with  the  long  axis  on  a  radial  line.  The  probe  volume  for  the 
radial  component  is  perpendicular  to  the  other  two  volumes,  with  the  long  axis  on  a  tangential  line. 

LDV  measurements  were  taken  upstream  and  downstream  of  the  rotor.  A  window  upstream  of  the  rotor 
permitted  the  measurement  of  axial  and  tangential  velocities  along  a  radial  line  upstream  of  the  rotor.  Figure 
2  shows  the  windows  installed  in  the  shroud  to  permit  measurement  of  the  tip-leakage  flow.  Three  velocity 
components  were  measured.  The  position  of  the  shaft  was  encoded  with  an  8192  counts/revolution  signal, 
which  was  recorded  with  each  velocity  measurement.  The  measurements  are  grouped  into  1024  circumfer¬ 
ential  positions,  each  8  encoder  counts  wdde.  At  each  measurement  point,  300,000  velocity  measurements 
were  made.  This  means  that  each  component  of  each  vector  on  a  plot  respresents  approximately  100  velocity 
measurements.  Doppler  signals  were  analyzed  with  a  TSI  Model  IFA  655  Digital  Burst  Correlator.  The 
processor  performs  a  256-sample,  double-clipped,  autocorrelation  on  each  Doppler  burst,  allowing  the  mea¬ 
surement  of  velocity  even  when  the  signal-to-noise  ratio  is  low.  The  processors  were  operated  in  the  random 
mode.  The  fiber-optic  LDV  system  consisted  of  two  3.25  inch  TSI  model  9832  fiber  optic  probes.  The  two 
probes  were  rigidly  mounted  together  on  a  traverse  that  could  translate  in  the  axial  and  radial  directions. 
The  horizontal  probe  utilized  the  green  (514.5  nm)  and  violet  (476.5  nm)  beams  of  the  argon-ion  laser  to 
measure  the  tangential  and  axial  components  of  velocity,  respectively.  The  vertical  probe  utilized  the  blue 
(488  nm)  beam  of  the  argon-ion  laser  to  measure  the  radial  component  of  velocity. 

Particle  Image  Velocimetry  The  PIV  measurements  were  made  using  two  Quanta-Ray  PRO-Series  Pulsed 
Nd:YAG  lasers  rated  at  800  mJ/pulse  output  at  532  nm.  The  laser  output  was  formed  into  a  light  sheet  of  5 
mm  thickness.  The  light  sheet  passed  through  ^windows  in  the  water  tunnel  and  duct  to  illuminate  the  flow, 
and  the  sheet  was  oriented  parallel  to  the  propeller-shaft  axis  (Figure  2).  The  axial  position  of  the  camera 
and  sheet  were  moved  to  place  the  mean  center  of  the  vortex  in  the  center  of  the  image.  Silcon  Carbide 
particles  of  mean  diameter  1  /xm  were  added  the  flow.  A  LaVision  Flowmaster  3S  PIV/PTV  system  was 
used  to  control  the  firing  of  the  lasers  and  synchronize  image  capture  with  a  digital  imager,  a  1280  x  1024 
pixel  cross-correlating  camera  with  12  bit  resolution.  The  camera  axis  was  perpendicular  to  the  laser  sheet 
and  recorded  the  in-plane  motion  of  the  tracer  particles.  The  light  sheet  was  sufficiently  thick  to  reduced 
the  number  of  particles  that  entered  or  exited  the  light  sheet  due  to  strong  cross-plane  flow  velocities.  The 
camera  was  contained  in  a  water-proof  housing  within  the  test  section  of  the  water  tunnel. 

The  image  field-of-view  was  24  by  30  mm,  and  the  lenses  on  the  camera  were  chosen  to  have  a  depth-of- 
field  much  larger  than  the  thickness  of  the  light  sheet.  The  PIV  images  were  spatially  calibrated  by  taking 
images  of  a  registration  target  in  the  image  plane.  Images  of  the  target  were  taken  in  the  filled  test  section. 

The  double  pulsed  PIV  images  were  analyzed  using  the  LaVision  software  DaVis  version  5.4.4.  An 
adaptive  multi-pass  image  process  was  employed  starting  with  interogation  windows  of  64  x  64  pixels  and 
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repeating  with  a  final  interrogation  windows  of  32  x  32  pixels.  The  spatial  resolution  of  the  velocity  vector 
field  is  0.75  mm.  There  were  typically  fewer  than  10%  bad  vectors  before  postprocessing.  A  three-by-three 
gaussian  smoothing  filter  was  used  on  the  vector  fields.  Averages  were  calculated  from  532  individual  PIV 
images. 

With  the  advantage  of  the  high  spatial  resolution  of  the  PIV  vector  fields,  it  was  possible  to  assess  the 
vortex  chairacteristics  such  as  its  core  radius  by  direct  measurement.  The  core  radius  was  defined  as  the 
radius  of  a  circle  with  an  area  equivalent  to  that  enclosed  by  the  polygon  centered  at  the  point  of  minimum 
speed  and  formed  by  the  points  of  maximum  speed  magnitude  positions  (core  contour).  We  will  refer  to  the 
average  value  of  these  maximum  speeds  as  the  tangential  speeds,  VJ.  The  maximum  tangential  speed  value 
on  the  core  contour  is  denoted  The  strength  of  the  vortex  was  calculated  by  integrating  the  vorticity 

vector  over  the  core  area.  The  vorticity  vector  at  each  cell  (square)  was  calculated  from  the  following  formula, 
Wz  =  (fc  ^  which  is  the  circulation,  or  the  line  integral  of  the  velocity  around  a  closed  path,  c, 

en^closing  the  cell,  divided  by  the  area  enclosed  by  c.  The  path  c  is  a  concentric  square  with  the  cell  of  interest 
and  having  twice  the  side  length.  The  velocities  from  the  eight  neighboring  cells,  four  at  the  sides  and  four 
at  the  corners  of  the  cell  under  consideration  were  involved  in  the  calculation.  The  area  enclosed  by  the  path 
c  was  made  up  of  the  area  of  the  cell  under  consideration  plus  1/2  the  area  of  the  four  side  cells  plus  1/4  the 
area  of  the  four  corner  cells.  Three  methods  of  vector  field  averaging  are  presented  below:  simple-averaging, 
center-averaging  based  on  the  location  of  minimum  speed  near  the  vortex  core,  and  center-averaging  based 
on  the  position  of  maximum  vorticity  near  the  vortex  core. 


3  Flow  Upstream  of  the  Rotor 

Figure  3  shows  the  axial  and  tangential  flow  upstream  of  the  rotor.  The  flow  was  measured  at  x/ Rp  =  —0.362 
where  x  is  measured  along  the  axial  flow  direction,  and  x  =  0  corresponds  to  the  center  of  the  12-inch  (0,3048 
m)  long  propeller  hub.  Rp  is  the  radius  of  the  rotor.  The  LDV  velocity  measurements  were  non-coincident 
and  taken  at  various  radial  positions.  The  advance  coefficient  was  J  =  0.983  and  data  was  taken  at  250  and 
500  rpm.  Ux  is  the  velocity  in  the  axial  direction,  Ut  is  the  tangential  velocity,  and  Uq  is  the  tunnel  flow 
velocity.  The  black  line  shows  the  axial  velocity,  the  red  line  shows  the  tangential  velocity,  the  solid  points 
are  for  250  rpm,  and  the  outlined  points  are  for  500  rpm.  The  boundary  layers  on  the  duct  and  on  the  hub 
can  be  seen.  The  duct  boundary  layer  has  a  thickness  of  about  0.07i2p.  This  corresponds  to  about  5  times 
the  tip-clearance  dimension.  The  tangential  velocities  are  nearly  zero  near  the  tip. 


4  Average  Tip-Leakage  Vortex  Flow 

LDV  and  PIV  were  used  to  measure  the  average  flow  near  the  tip  of  the  blade.  Presented  here  are  planar 
velocity  fields.  The  plane  is  parallel  to  the  axis  and  corresponds  to  an  axial  distance  Ax/i?p  =  0.071  and  to 
the  radial  distance  Ar/iip  =  0,056.  The  axial  position  was  moved  to  center  the  vortex  for  images  taken  at 
different  blade  delays.  As  the  rotating  blade  passes  through  this  plane,  the  tip-leakage  vortex  intersects  the 
image  plane,  and  the  axis  of  the  vortex  is  almost  perpendicular  to  the  plane. 

Figure  4  shows  a  schematic  diagram  of  the  LDV  and  PIV  vector  fields  with  respect  to  the  duct  wall  and 
the  propeller  blade.  The  images  are  taken  such  that  the  blade  is  moving  into  the  paper.  The  leakage  vortex 
is  moving  from  the  pressure  side  of  the  blade  to  the  suction  side  showing  a  counterclockwise  rotation  when 
viewed  from  below  (as  the  PIV  camera  imaged  the  flow.).  The  duct  is  at  the  top  of  the  image  and  the  mean 
tunnel  flow  is  moving  from  left  to  right  through  the  plane.  The  angular  position  of  the  rotor,  0,  is  measured 
with  respect  to  this  plane.  When  0  equals  zero,  the  trailing  edge  of  the  blade  intersects  the  plane,  and  as  0 
increases,  the  blade  moves  farther  away.  This  distance  is  presented  as  5  =  ORjC  where  C  is  the  chord  of 
the  blade. 

It  is  interesting  to  compare  the  average  flow  fields  measured  with  LDV  and  PIV.  Figure  5  presents 
velocity  fields  measured  with  LDV,  on  the  left,  and  PIV,  on  the  right.  The  LDV  is  shown  with  only  2- 
components  of  velocity  for  comparison  with  the  PIV.  The  PIV  has  been  simple-averaged  over  532  images  for 
comparison  with  LDV  measurements.  The  PIV  velocity  field  shows  good  spatial  resolution  and  matches  the 
LDV  measurements. 
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Propiilsor  Inflow  at  250  and  500  rpm 


Figure  3:  Upstream  non-coincident  LDV  velocity  measurements  at  various  radial  positions  at  an  advance 
coefficient  of  J  =  0.983and  at  an  axial  location  of  x/Rp  =  --0.3615.  Ux  is  the  velocity  in  the  axial  direction, 
Ut  is  the  tangential  velocity,  and  Uo  is  the  tunnel  flow  velocity.  These  measurements  were  taken  at  500  and 
250  rpm.  The  black  line  is  the  axial  velocity  and  the  red  line  is  the  tangential  velocity.  The  solid  points  are 
for  250  rpm  and  the  outlined  points  are  for  500  rpm. 


Duct  Wall 


Axial  Row 
Direction 


Figure  4:  Schematic  diagram  of  LDV  and  PIV  velocity  fields.  Note  that  the  images  are  taken  such  that  the 
blade  is  moving  into  the  paper  (see  Figure  2). 
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Figure  5:  LDV  and  PIV  velocity  vectors  for  S-0.0296  at  500  rpm  and  45  psia.  The  contour  color  shows 
the  vector  magnitude.  xJRp  is  the  non-dimensional  distance  in  the  axial  direction  and  v/Rp  is  the  non- 
dimensional  distance  in  the  radial  direction.  The  PIV  is  averaged  over  532  images.  The  LDV  is  shown  with 
only  2-components  of  velocity  for  comparison  with  the  PIV.  The  LDV  has  a  covers  a  larger  domain,  but  the 
PIV  shows  good  spatial  resolution  and  matches  the  LDV  velocity  measurements. 


It  is  important  to  remember  that  the  vortical  flow  near  the  trailing  edge  is  typically  three  dimensional, 
and  that  the  PIV  image  is  only  a  two-dimensional  cut  of  this  flow  averaged  across  the  PIV  imaging  light 
sheet.  This  is  illustrated  in  Figure  6.  Here  a  comparison  is  made  between  2-component  and  3-component 
LDV  measurements.  On  the  left  is  a  2-component  velocity  field  parallel  to  the  PIV  measurement  plane.  On 
the  right  is  the  velocity  field  of  a  plane  that  was  slightly  tilted.  Because  the  vortex  does  not  necessarily  cut 
the  PIV  measurement  plane  at  a  right  angle,  the  vortex  on  the  left  appears  less  axially  symmetric.  The  axial 
symmetry  is  recovered  after  tilting  the  plane.  Consequently,  it  is  important  not  to  interpret  the  PIV  images 
as  perpendicular  cuts  of  a  two-dimensional  vortex.  The  PIV  imaging  plane  remained  parallel  to  the  axial 
flow  direction,  and  the  vortex  may  not  intersect  it  at  a  right  angle. 


5  Variability  of  the  Tip-Leakage  Vortex  Flow 


While  only  a  2-diinensional  realization  of  the  flow,  the  advantage  of  the  PIV  visualization  is  the  possibility 
of  capturing  many  instantaneous  flow  fields.  The  tip  leakage  flow  is  highly  unsteady,  as  evidenced  by  the 
variation  in  single  PIV  images.  Figure  7  shows  four  PIV  images  taken  at  the  same  operating  conditions  and 
downstream  location.  All  vector  plots  axe  for  the  same  blade  angles  or  blade  positions,  the  trailing  edge  is 
1.52®  (5  =  0.0296)  and  27.02®  (S  =  0.5259)  above  the  plane  of  measurement.  Variation  in  the  flow  field  can 
be  readily  observed.  The  position,  strength,  and  core  size  of  the  vortex  varies  between  images.  The  contours 
show  changes  in  the  measured  out-of-plane  vorticity. 


CAV2001:sessionA6.001 


8 


Figure  6:  2'Component  versus  3-component  LDV  measurements.  These  measurements  were  taken  at  500  rpm, 
45  psia,  and  8=0.0296.  The  2-component  LDV  gives  similar  results  to  the  simple-averaged  PIV.  The  contour 
color  corresponds  to  vector  length.  The  3-component  measurements  result  in  a  more  accurate  mean  core 
structure  of  the  vortex.  The  vortex  is  not  necessarily  oriented  with  the  horizontal  plane.  The  3-component 
LDV  measurements  allow  a  non-horizontal  cut  to  be  taken  so  that  the  vortex  is  perpendicular  to  the  in-plane 
velocity  measurements.  r/Rp  and  x/Rp  are  shown. 


T^ble  3.  Comparison  between  simple-averaged  PIV  and  center-averaged  PTV  vortex  quantities. 


5 

Type  of 
average 

Radius, 
a  (mm) 

Circulation, 
r  (m“/s) 

'.2)ra-' 

{misf 

Average  tangential 
velocity,  VJ  (m/s) 

Maximum  tangential 
velocity,  (m/s) 

0.0296 

simple 

average 

8.695 

-0.7146 

-171.1 

13.16 

18.01 

0.0296 

velocity 

centered 

average 

8.511 

-0.7026 

-172.6 

13.14 

17.90 

0.0296 

vorticity 

centered 

average 

8.630 

-0.7034 

-168.3 

13.12 

17.85 

0.5259 

simple 

average 

12.62 

-0.6193 

-61.02 

8.506 

11.32 

0.5259 

velocity 

centered 

average 

12.27 

-0.6024 

-61.06 

8.440 

11.12 

0.5259 

vorticity 

centered 

average 

12.74 

-0.5885 

-54.03 

8.060 

10.56 
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Figure  11  shows  histograms  of  the  vortex-center  location  for  S  =  0.0296  and  S  =  0.5259.  The  location  of 
the  vortex  core  varied  by  ±1  mm  and  ±3  mm,  respectively,  or  11%  and  26%  of  the  vortex  core  radii.  After 
identifying  the  position  of  the  vortex,  it  is  possible  to  shift  the  PIV  images  so  that  a  centered-average  can 
be  computed.  Table  3  shows  three  averages  computed  for  532  PIV  images.  The  simple-average  is  just  the 
non-shifted  average  of  the  PIV  images.  The  velocity-centered  average  used  the  position  of  minimum  velocity 
near  the  vortex  core  as  the  locus  for  the  average,  and  the  vorticity-centered  average  used  the  position  of 
maximum  vorticity  near  the  vortex  core.  Figure  9  shows  the  averaged  PIV  images  for  these  two  blade  values, 
S.  For  the  case  of  S  =  0.0296,  the  effect  of  vortex  wandering  is  small,  with  changes  in  the  core  size  of  the 
centered-averages  less  than  2%  compared  to  the  simple-average.  The  wandering  is  greater  for  the  case  of 
S  =  0.5259,  with  the  differences  approaching  3%.  These  results  are  consistent  with  the  calculations  of 
Straka  and  Farrell  (1992).  They  computed  the  effect  of  vortex  wandering  on  the  resulting  averaged  flow  field 
of  a  two-dimensional  Burgers  vortex  and  showed  that  the  amplitude  of  vortex  must  wander  by  greater  50% 
of  the  core  radius  for  the  simple-average  measurement  to  depart  from  the  center-averaged  measurement. 

Figure  11  presents  histograms  of  the  core  radius,  a,  the  circulation  around  the  viscous  core,  F,  the  average 
tangential  speed,  Vt,  the  maximum  value  for  the  tangential  speed,  and  the  quantity  “(r/27ra)^,  for 

two  blade  angles.  For  S  =  0.0296,  the  core  radius  ranges  between  7.40  and  10.12  mm,  the  average  circulation 
is  between  -0.82  and  -0.66  m^/s,  and  the  average  tangential  velocities  were  between  12.82  and  14.96  m/s 
with  the  maximum  tangential  velocities  between  16.66  and  21.01  m/s. 

As  discussed  above,  variation  in  the  angle  that  the  vortex  axis  makes  with  the  PIV  light-sheet  can  lead  to 
variation  in  the  measured  quantities  of  the  vortex.  Consider  a  line  vortex  with  radius  gq  that  makes  an  angle 
of  a  with  the  PIV  plane.  The  vortex  core  will  be  imaged  as  an  ellipse  with  area  Trag  sec(a).  The  measured 
average  core  diameter  would  then  be  larger  than  the  actual  core  diameter  by  a  factor  of  (sec(Q:))^/^.  A  10° 
variation  in  the  incidence  angle  will  result  in  a  1%  increase  in  the  measured  core  size,  and  a  20°  inclination 
will  result  in  a  3%  variation.  The  actual  angle  of  incidence  is  likely  to  be  much  less  than  ±20°.  Consequently, 
the  observed  variability  in  the  instantaneous  PIV  quantities  is  not  principally  due  to  vortex  wandering  but 
is  the  result  of  true  variations  in  strength  and  size  of  the  vortex. 


6  Development  of  the  Tip-Leakage  Vortex  Flow  Downstream  of 
the  Blade 

The  vortex  core  size  and  strength  vary  as  the  vortex  flow  moves  downstream  of  the  blade  trailing  edge. 
Figure  10  shows  these  quantities  plotted  as  a  function  of  5.  Variation  of  these  curves  are  relatively  smooth, 
except  for  the  region  just  downstream  of  the  trailing  edge  where  the  leakage  vortex  is  rolling  up  and  may 
be  interacting  with  a  trailing  edge  vortex.  For  S  >  0.2312  the  variation  is  quite  smooth.  The  radius  of  the 
core  size  increases  as  the  blade  moves  up  and  away  from  the  plane  of  measurement.  The  circulation  shows  a 
generally  increasing  trend  as  the  trailing  edge  moves  further  out  of  the  plane  of  measurement.  The  average 
tangential  velocity  and  the  maximum  tangential  velocity  show  a  peak  right  when  the  trailing  edge  of  the 
blade  is  in  the  plane  of  measurement  (5  =  0)  and  then  decrease  as  the  blade  moves  further  away.  Lastly, 
the  quantity  (r/27ra)^  is  largest  near  the  trailing  edge  and  is  reduced  farther  downstream.  These  data  can 
also  be  derived  from  the  3-dimensional  LDV  data,  along  with  the  axial  velocity  within  the  vortex  core.  We 
intend  to  compare  the  PIV  and  LDV  results  more  closely. 


7  Developed  Tip-Leakage  Vortex  Cavitation 

The  cavitating  vortex  core  was  visualized  by  reducing  pressure  beyond  inception,  as  shown  in  Figures  12  and 
13.  Visualization  of  the  vortex  shows  that  it  is  angled  relative  to  the  PIV  plane  by  approximately  15°  in  the 
upstream  axial  direction  near  the  trailing  edge  and  about  3°  upstream  of  the  trailing  edge.  Surface  cavitation 
can  be  observed  in  the  gap  region  for  the  developed  cavitation.  High-speed  video  showed  wandering  of  the 
vortex  with  a  displacement  that  ranged  from  very  little  when  the  vortex  was  a  along  side  of  the  blade  (less 
then  one-half  the  core  width)  to  100%  of  the  observed  vapor  core  diameter  about  0.1524  m  (6in)  downstream. 
This  is  consistent  with  the  wandering  amplitude  of  the  non-cavitating  vortex.  Downstream  of  that  the  vortex 
wanders  more,  probably  due  to  the  combining  on  the  trailing  edge  vortex  and  general  wake  roll  up.  At  this 
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S=0.0296 


S=0.5259 


Figure  8:  Histograms  of  the  (X,Y)  position  of  the  vortex  center  showing  the  amount  of  wandering  at  S=0.0296 
and  S=0.5259,  500  rpm  and  45  psia.  The  PIV  measured  the  vortex  core  over  a  range  of  5.36  (mm)  in  the 
x-direction  and  3.06  (mm)  in  the  {/-direction  for  S  =  0.0296.  The  wandering  increases  further  from  the 
blade. 


point  the  vortex  gets  much  more  distorted  and  wanders  up  to  300%  of  the  observed  vapor  core. 

Arndt  and  Keller  (1992)  employed  the  principal  of  angular  momentum  conservation  to  show  that  the 
diameter  of  the  vapor  core  of  a  developed  cavitating  Rankine  vortex  is  smaller  than  the  non-cavitating 
vortex  core  by  a  factor  of  l/>/2.  Here  the  non-cavitating  core  radius  was  measured  to  be  9  mm  downstream 
of  the  trailing  edge,  and  the  measured  vapor  core  had  a  radius  of  approximately  3  mm,  or  0.3  times  the 
original  core  radius.  This  is  much  lower  than  the  expected  value  of  l/\/2  =  0.71.  Similar  measurements  of 
the  cavitating  to  non-cavitating  core  ratio  for  tip  vortex  forming  on  a  stationary  elliptic  hydrofoil  yielded  a 
ratio  of  0.5  (Briangon-Marjollet  and  Merle,  1997). 

Figure  13  shows  a  detail  of  the  leakage  vortex,  where  surface  cavitation  at  the  tip  convects  into  the  vortex. 
The  flow  direction  can  be  observed  along  with  the  increase  in  cavitating  vorticity,  fed  by  the  gap  flow. 


8  Conclusions 

The  tip-leakage  vortex  occurring  on  a  ducted  rotor  was  examined  using  both  LDV  and  PIV.  Preliminary 
analysis  of  the  data  indicates  that  the  there  is  substantial  variability  in  the  instantaneous  vortex  quantities, 
such  as  the  size  of  the  viscous  core  and  the  strength  of  the  vortex.  Inception  of  discreet  cavitation  events 
will  occur  when  nuclei  are  captured  within  the  vortex,  experience  tension,  and  vaporously  grow.  Inception 
can  occur  when  an  unusually  large  nucleus  is  captured  by  the  average  vortex,  and  when  an  average  nucleus 
is  captured  by  an  unusually  strong  vortex.  We  are  presently  conducting  inception  experiments  to  determine 
how  variability  of  the  vortex  influences  inception  of  discreet  cavitation  events.  We  will  examine  when  and 
where  discreet  cavitation  events  occur,  and  attmpt  to  relate  this  to  the  measured  two-  and  three-dimensional 
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S=0.5259 

(a)  (b)  (c) 


Figure  9:  Averaged  and  centered-averaged  velocity  vectors  and  vorticity  contours  from  532  PIV  images  at 
500  rpm,  45  psia,  and  S  =  0.0296  and  0.5259.  An  average  where  the  center  of  each  vortex  is  in  the  same 
location  is  called  ’centered- average’.  The  simple-averaged  image  shows  a  larger  vortex  core  size  and  lower 
tangential  velocities.  Table  3  shows  variations  in  vortex  quantities  between  the  two  averaging  methods. 


flow  quantities. 
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S=0.0296 


S=0.5259 


Figure  11:  Distributions  of  radius  (R),  circulation  (F),  average  tangential  velocity  (Vt),  maximum  tangential 
velocity  (Vt^ax)>  quantity  “(r/27ra)^  for  532  PIV  images  for  S=0.0296  and  S=0.5259  at  500  rpm, 

45  psia.  Note  the  variation  in  vortex  quantities  for  a  single  blade  angle. 
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ABSTRACT 

The  sound  produced  by  the  collapse  of  discreet  cavitation  bubbles  was  examined.  Laser¬ 
generated  cavitation  bubbles  were  produced  in  both  a  quiescent  and  a  vortical  flow.  The  sound 
produced  by  the  collapse  of  the  cavitation  bubbles  was  recorded,  and  its  spectral  content  was 
determined.  It  was  found  that  the  rise  time  of  the  sound  pulse  produced  by  the  collapse  of  single, 
spherical  cavitation  bubbles  in  quiescent  fluid  exceeded  that  of  the  slew-rate  of  the  hydrophone, 
which  is  consistent  with  previously  published  results.  It  was  found  that,  as  collapsing  bubbles 
were  deformed  by  the  vortical  flow,  the  acoustic  impulse  of  the  bubbles  was  reduced. 
Collapsing  non-spherical  bubbles  often  created  a  sound  pulse  with  a  rise  time  that  exceeded  that 
of  the  hydrophone  slew-rate,  although  the  acoustic  impulse  created  by  the  bubbles  was 
influenced  largely  by  the  degree  to  which  the  bubbles  became  non-spherical  before  collapse. 
The  noise  produced  by  the  slow  growth  of  cavitation  bubbles  in  the  vortex  core  was  not 
detectable.  These  results  have  implications  for  the  interpretation  of  hydrodynamic  cavitation 
noise  produced  by  vortex  cavitation. 
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INTRODUCTION 

The  sound  associated  with  the  growth  and  collapse  of  single  cavitation  bubbles  has  been  a 
topic  of  interest  for  some  time  [1-4].  Discreet  cavitation  events  are  often  the  source  of 
significant  cavitation  noise  in  underwater  flows  and  within  turbo-machinery.  Over  time, 
researchers  have  undertaken  to  understand  how  cavitation  bubbles  create  noise  and  how  that 
noise  can  be  scaled  with  changes  in  cavitation  and  Reynolds  number. 

A  widely  used  scaling  method  for  hydrodynamic  cavitation  noise  was  proposed  by 
Fitzpatrick  and  Strasberg  [5]  in  1957.  They  used  spherical  bubble  dynamics  to  analyze  the 
volume  growth  and  collapse  of  a  single  cavitation  bubble  in  an  infinite  fluid.  Fourier  analysis  of 
the  volume  acceleration  was  then  used  to  determine  the  spectral  content  of  the  noise  emitted  by 
the  bubble.  A  time  scale  of  the  problem  is  the  Rayleigh  collapse  time,  , 

tc  =0.915i?^VP^ 

where  R^^is  the  maximum  radius  of  the  bubble  (sometimes  taken  as  the  Equation.3  I AP}  is  the 
pressure  difference  between  the  external  fluid  and  the  bubble  contents,  and  p  is  the  liquid 
density.  With  monopole  sound  generation,  the  acoustic  pressure  generated  by  the  bubble  will 
scale  with  the  volume  acceleration,  an  estimate  of  which  can  also  be  derived  from  Rayleigh’s 
analysis.  The  sound  pressure  produced  by  the  growth  and  collapse  of  the  bubble  is  then  given  by 
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where  t  is  the  distance  between  the  bubble  and  the  sound  pressure  measurement  [6]. 

The  Fitzpatrick  and  Strasberg  model  predicts  that  the  sound  spectrum  generated  by  single 
bubble  collapse  should  scale  as  {(Otci'iox  (otc<\.  The  peak  of  the  spectrum  is  near  =  1 ,  and 
higher  frequencies  will  roll  off  as  {0)1^)'^.  Mellen’s  [7]  measurements  of  bubble  cavitation  noise 
spectra  resulting  from  a  bubble  cloud  formed  at  the  end  of  a  rotating  rod  exhibit  theses  spectral 
features  and  has  been  used  to  support  the  validity  of  this  analysis.  Corrections  to  this  analysis 
have  included  the  effect  of  liquid  compressibility  and  the  introduction  of  non-condensable  gas 
within  the  cavity  [8,9].  These  processes  tend  to  limit  the  predicted  energy  emitted  at  frequencies 
greater  than  10  •  • 

If  the  external  bubble  pressure  greatly  exceeds  the  pressure  of  the  non-condensable 
bubble  contents  at  the  time  of  the  bubble  maximum  volume,  the  bubble  wall  velocity  can 
approach  the  local  liquid  sound  speed,  leading  to  the  formation  of  a  pressure  shock-wave  during 
the  rapid  volume  accelerations  near  the  point  of  minimum  bubble  volume.  Baiter  [10]  has  shown 
that  such  a  shock  wave  can  exhibit  a  pressure  pulse  of  the  form 

=  (3) 

The  peak  pressure  is  Ps{r)  =  and  the  pulse  width  is 

6ir)  =  fz{r,R^,R„,Pc,pc^)-  Rm^  and  R„  and  maximum  and  minimum  bubble  radii,  Pq,  is  the  non¬ 
condensable  gas  pressure  of  the  bubble  at  its  maximum  volume,  and  c  is  the  liquid  sound  speed. 
The  spectrum  associated  with  the  sound  pulse 
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jPsd) 

1  +  (W0)^ 


(4) 


which  is  flat  until  the  cut-off  frequency  l/iojG).  The  cut-off  frequency  is  strongly  affected  by  the 
amount  of  non-condensable  gas  present  in  the  collapsing  cavity.  The  non-condensable  gas  may 
have  existed  within  the  original  gas  nucleus,  and  diffusion  to  or  from  the  liquid  may  either 
increase  or  reduce  the  total  mass  of  non-condensable  gas,  depending  on  the  pressure/volume 
history  of  the  bubble  and  saturation  level  of  the  surrounding  liquid  prior  to  bubble  collapse. 

Measurement  of  noise  spectra  from  discreet  hydrodynamic  cavitation  has  been  reported. 
Typically,  the  noise  produced  by  traveling  bubble  cavitation  has  been  examined  [11-15].  The 
pressure  pulse  created  by  a  collapsing  cavitation  bubble  can  be  characterized  by  the  acoustic 
impulse,  I: 

'-JCV'  (5) 

where  6  is  the  pulse  width.  However,  shearing  and  fission  of  the  bubble  can  significantly  reduce 
the  total  emitted  sound  [16].  The  presence  of  nearby  bubbles  can  also  lead  to  significantly 
reduced  acoustic  emission  [17]. 

Measurements  of  the  spectral  content  of  discreet  hydrodynamic  cavitation  events  have  not 
been  consistent.  Reverberations  in  confined  flow  facilities  and  the  limitations  of  the  fidelity  and 
bandwidth  of  the  measurement  equipment  can  contaminate  the  measured  spectrum.  Much  of  the 
spectral  data  presented  consists  of  the  average  spectra  produced  by  many  cavitation  events. 
Beside  the  data  reported  by  Mellen,  there  is  little  evidence  that  the  spectral  content  of  single 
cavitation  bubble  noise  will  scale  as  predicted  by  incompressible  bubble  dynamics.  Spectra 
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measured  from  traveling  bubble  and  vortex  cavitation  are  not  consistent  [12-15].  Some  spectra 
are  reported  to  exhibit  the  features  predicted  by  Fitzpatrick  and  Strasberg  (i.e.  low  Mach  number 
bubble  collapse),  while  others  are  reported  to  be  broadband.  Moreover,  the  relative  noise 
contribution  due  to  the  growth  and  collapse  phases  of  hydrodynamic  bubble  cavitation  are  often 
considered  comparable,  although  most  researchers  report  experimental  noise  measurements 
principally  from  the  bubble  collapse.  Further  complications  arise  from  the  generally  complex 
and  non-spherical  shape  of  hydrodynamic  cavitation  bubbles,  the  reverberant  environment  where 
experiments  are  conducted,  and  bandwidth  limits  of  measurement  equipment. 

Capture  of  the  true  acoustic  transients  created  by  collapsing  cavitation  bubbles  usually 
requires  the  construction  of  a  high  frequency  response  transducer,  typically  made  from  a  thin 
film  of  polyvinylide  fluoride  (PVDF).  The  pressures  measured  with  these  devices  has  been 
compared  to  the  optically  measured  volume  change  of  the  bubble,  and  it  has  been  shown  that,  by 
far,  the  strongest  noise  is  emitted  when  the  bubble  reaches  a  minima  of  volume  [18,19].  Sound 
is  produced  as  a  shock  wave  that  radiates  away  from  the  bubble,  and  the  rise  time  of  the  shock 
wave  is  on  the  order  of  5  nanoseconds,  the  maximum  peak  pressures  on  the  order  of  hundreds  of 
atmospheres  near  the  bubble,  and  the  pulse  width  is  on  the  order  of  10  to  100  ns.  The  collapse 
of  non-spherical  bubbles  can  produce  less  noise.  Bubbles  collapsing  near  a  solid  boundary  can 
be  deformed,  and  the  lack  of  spherical  focusing  can  significantly  reduce  the  emitted  noise  [18]. 
The  sharp  rise  time  high  pressures  observed  in  these  studies  indicate  that  the  spectrum  of  the 
emitted  sound  is  dominated  by  the  collapse  processes  and  is,  largely,  broadband. 

In  the  present  work,  we  will  measure  the  noise  pulse  and  spectra  associated  with  single 
cavitation  bubble  events  using  techniques  that  are  often  employed  in  the  study  of  hydrodynamic 
cavitation.  Bubbles  will  be  produced  in  the  vicinity  of  a  low-pressure  vortex  core.  We  will 
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show  that  the  noise  emitted  by  the  focused  collapse  of  distorted  bubbles  is  broadband  and  does 
not  exhibit  significant  spectral  content,  even  for  highly  distorted  bubbles.  The  noise  produced 
by  highly  distorted  cavitation  bubbles  is  still  scaled  by  the  dynamics  of  spherically  collapsing 
bubbles.  Lastly,  we  will  show  that  the  magnitude  of  the  noise  produced  during  bubble  growth  is 
significantly  smaller  that  the  noise  produced  upon  collapse. 

I.  EXPERIMENT 

Single  cavitation  bubbles  were  created  in  both  a  quiescent  liquid  and  in  a  vortical  flow. 
Figure  la  shows  the  schematic  diagram  of  the  quiescent  cavitation  bubble  experiment.  A  small 
cubical  Lucite  water  tank  (205  mm  side)  was  created  with  an  optical  glass  window  (100  mm 
diameter)  for  the  introduction  of  laser  light.  The  infrared  light  pulse  (1064  nm  wavelength,  ~10 
nS  pulse  duration,  280  mJ/pulse,  maximum)  was  produced  by  a  Q-switched  Nd:YAG  laser 
(Spectra  Physics  PRO  250).  Through  a  set  of  spherical  lenses,  the  beam  was  expanded  to  a 
diameter  of  100  mm  and  then  focused  to  a  waist  within  the  tank.  The  cone  angle  was  19.4°  in 
air.  The  space  over  the  liquid  in  the  tank  could  be  evacuated,  and  the  gas  pressure  was  recorded 
with  an  Omega  PX203  pressure  transducer.  Distilled  water  was  placed  within  the  tank.  In  order 
to  increase  the  repeatability  of  bubble  production,  a  small  amount  of  drinking  tap  water  was 
added  to  the  tank.  The  tap  water  contained  enough  micron  sized  contaminants  that  led  to  the 
regular  absorption  of  the  laser  light  in  the  location  of  the  beam  waist  and  the  production  of  a 
single  cavitation  bubble,  although  intermittent  production  of  multiple  bubbles  was  also  observed. 
Variation  in  the  static  pressure  above  the  tank  and  the  energy  of  the  light  pulse  resulted  in 
variation  of  the  maximum  bubble  volume  and  collapse  time  of  the  cavitation  bubble. 

Images  of  the  bubble  were  observed  with  a  Cooke  Corporation  “FlashCam”  digital  camera. 
This  camera  has  an  effective  resolution  of  752  by  240  pixels,  and  it  is  capable  of  recording  up  to 
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10  consecutive  exposures  on  the  same  image  frame,  with  varying  integration  times  and  delays 
between  each  exposure.  A  pulse  delay  generator  (Stanford  Research  Systems  model  DG535) 
was  used  to  trigger  the  camera  at  varying  delays  from  the  laser  Q-switch  pulse.  These  images 
were  then  recorded  digitally.  The  acoustic  center  of  a  Briiel  &  Kjasr  type  8103  hydrophone  was 
placed  30  mm  from  the  location  of  the  beam  waist,  and  the  resulting  sound  emission  was 
conditioned  with  a  Bruel  &  Kjaer  type  2635  charge  amplifier  and  digitally  acquired  using  a 
TekTronix  storage  oscilloscope.  The  A/D  was  12-bit  resolution,  and  the  sampling  rate  was  2.5 
million  samples/second. 

Figure  lb  shows  the  setup  for  the  creation  of  single  cavitation  bubbles  in  a  vortical  flow.  The 
water  tunnel  has  a  circular  contraction  downstream  of  a  series  of  flow  management  screens  with 
area  contraction  ratio  6.4:1.  The  test  section  has  a  22.9  cm  (9  inch)  diameter  round  inlet  that  is 
then  faired  into  a  rectangular  test  section  with  widely  rounded  comers.  Four  acrylic  windows 
(93.9  cm  by  10.0  cm  viewing  area)  permit  optical  access  to  the  test  section  flow.  The  flow  in  the 
test  section  can  be  operated  at  pressures  from  vapor  pressure  to  approximately  200  kPa  (30  Psia). 
The  average  velocity  in  the  test  section  is  variable  up  to  18  m/s  (59.05  ft/s).  A  de-aeration 
system  can  be  used  to  vary  the  dissolved  gas  content  of  the  flow,  and  the  inlet  water  is  filtered  to 
1  microns. 

A  vortical  flow  was  created  using  a  cambered  hydrofoil  mounted  to  one  window  of  the  test 
section.  The  hydrofoil  has  a  rectangular  planform  of  9.5  cm  span  and  16.8  cm  chord,  and  the  tip 
of  the  hydrofoil  was  truncated  with  sharp  edges.  The  hydrofoil  mount  allows  continuous 
changes  of  the  incident  flow  angle.  A  series  of  tip  and  trailing  edge  vortices  will  be  shed  near 
the  tip,  and  these  vortices  will  merge  to  form  a  single  vortex  within  one-half  chord  length 
downstream  of  the  hydrofoil  trailing  edge.  The  tip  vortex  produced  by  the  hydrofoil  can  be 
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visualized  with  developed  cavitation,  as  shown  in  Figure  2.  Measurements  of  the  bubble/vortex 
interactions  were  conducted  using  a  free-stream  velocity  of  10  m/s  and  a  variety  of  pressures. 
The  dissolved  oxygen  content  was  measured  with  an  Orion  Model  810  dissolved  oxygen  meter. 
In  order  to  reduce  the  number  of  free-stream  nuclei,  the  free-stream  gas-content  was  reduced  to 
below  1.5  ppm  during  the  measurements. 

Planar  Particle  Imaging  Velocimetry  (PIV)  was  used  to  measure  the  vortical  flow  field  at  a 
station  9.0  cm  downstream  of  the  trailing  edge.  A  double-pulsed  light  sheet  9  mm  thick  was 
created  perpendicular  to  the  mean  flow  direction  using  two  pulsed  Nd:YAG  lasers  (Spectra 
Physics  model  Pro-250  Series).  15-micron  average  diameter  silver  coated  glass  spheres  (from 
Potters  Industries)  were  used  to  seed  the  flow.  An  acrylic  prism  was  optically  mounted  to  a 
window  of  the  test-section  for  viewing  of  the  light  sheet  with  reduced  optical  distortion.  The 
light  sheet  was  imaged  with  a  PIV  image  capture  system  produced  by  LaVision  Inc.  Double- 
pulsed  images  of  the  light  sheet  were  acquired  with  a  digital  camera  with  1280  x  1024  pixels. 
Optical  distortion  of  the  planar  light  sheet  image  was  corrected  through  a  calibration  procedure 
that  employed  the  imaging  of  a  regular  grid  in  the  location  of  the  light  sheet  plane.  Velocity 
vectors  were  produced  from  the  double-pulsed  images  using  the  LaVision  image  analysis 
software  Davis6.0.4.  Multi-pass  processing  with  a  final  window  size  of  32  x  32  pixels  was  used 
with  12%  window  overlap  in  the  final  pass  to  produce  41  x  27  in  plane  velocity  vectors  at  0.62 
mm  spacing.  Since  the  camera-imaging  plane  was  not  parallel  to  the  light  sheet  and  it  had  an 
angular  shift  of  45°  with  the  horizontal  direction,  this  velocity  component  was  corrected  by 
knowledge  of  the  camera  angle  and  the  tunnel  free  stream  velocity. 

Optically  stimulated  vortex-cavitation  bubbles  were  created  at  the  vortex  centerline.  As  the 
cavitation  number  was  reduced,  discreet  cavitation  events  occurred  in  the  vortex  core  due  to  the 
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capture  of  free-stream  nuclei.  This  cavitation  occurred  at  various  locations  along  the  vortex  axis 
downstream  of  the  location  of  vortex  roll-up.  However,  at  such  lower  air  contents,  the  event 
rate  of  this  cavitation  was  quite  low  (<  1  event  per  second).  Vortex  cavitation  could  then  be 
stimulated  by  the  introduction  of  a  nucleus  though  a  focused  pulse  of  IR  laser  light  near  the  axis 
of  the  vortex.  By  varying  the  free-stream  pressure  and  laser  energy,  the  initial  size  of  the  laser- 
induced  bubble  could  be  controlled.  Images  of  the  bubble’s  growth  and  collapse  were  captured 
with  the  pulse-synchronized  digital  camera  after  an  appropriate  delay,  as  described  above.  A 
small  cast  acrylic  tank  of  water  was  placed  in  pool  of  water  on  the  top  window  of  the  test  section, 
and  in  this  tank  the  hydrophone  was  mounted.  This  produced  an  acoustic  path  between  the 
cavitation  event  and  the  hydrophone  of  relatively  matched  impedance,  although  acoustic 
impedance  mismatches  between  the  acrylic  and  the  water  led  to  some  attenuation  and  internal 
reflections. 

II.  RESULTS-  CAVITATION  IN  QUIESCENT  FLUID 

As  discussed  in  the  introduction,  several  researchers  have  examined  in  detail  the  noise 
emitted  by  single  cavitation  events  in  quiescent  fluid.  They  have  shown  that  the  amplitude  and 
frequency  content  of  the  noise  emission  often  greatly  exceeds  the  capabilities  of  standard  high- 
frequency  hydrophones.  In  the  present  work,  we  first  recreated  the  basic  elements  of  these 
experiments  to  characterize  our  process  of  bubble  formation  and  noise  measurement  by 
examining  spherical  and  near-spherical  bubbles  were  created  in  the  small  tank.  The  effect  of 
buoyancy  can  be  characterized  by  the  dimensionless  quantity 


OWEIS,  JASA 


where  g  is  the  gravitational  acceleration  [19].  For  sufficiently  small  values  of  8,  the  effect  of 
buoyancy  can  be  neglected.  Spherical  bubbles  corresponding  to  6  «  0.01  were  created  to 
characterize  the  hydrophone’s  ability  to  record  the  noise  pulse  produced  by  a  single  vaporous 
cavitation  event. 

Figure  3a  shows  multiple  images  of  a  near  spherical  bubble  created  with  a  laser  light  pulse. 
The  hydrophone  acoustic  pressure  signal  produced  by  the  collapse  of  the  spherical  bubble  is 
shown  in  Figure  3b.  Two  distinct  sound  pulses  are  detected.  The  first  occurs  during  the  rapid 
growth  of  the  bubble  volume  when  the  bubble  is  created.  Between  the  pulses,  a  slow  reduction 
and  increase  in  pressure  is  seen  that  corresponds  to  the  relatively  slow  volume  acceleration 
around  the  point  of  maximum  bubble  volume.  The  second  pulse  corresponds  to  the  first  collapse 
of  the  bubble.  Figure  3c  shows  a  time-expanded  pressure  trace  recorded  for  the  collapse  of  the 
bubble.  The  response  is  that  of  a  single-degree-of-freedom  linear  oscillator.  Also  plotted  is  the 
response  of  a  single-degree-of-freedom  oscillator  excited  with  a  narrow  pulse  where  the  natural 
frequency  and  damping  of  the  oscillator  was  chosen  to  best  fit  the  measured  voltage  trace.  The 
natural  frequency  of  the  oscillator,  was  calculated  to  be  105  kHz,  and  the  damping 
coefficient,  C,  was  0.18.  The  natural  frequency  is  in  good  agreement  with  the  specifications  of 
the  hydrophone,  which  reported  a  resonance  frequency  of  nearly  109  kHz.  Variation  in  laser 
input  power  changed  the  maximum  bubble  volume  and  the  corresponding  maximum  pulse 
amplitude,  but  the  pulse  response  of  the  hydrophone  did  not  vary.  This  was  also  true  for  bubbles 
undergoing  non-spherical  collapse  as  well.  Rebounding  bubbles  and  bubbles  collapsing  near  a 
wall  could  sometimes  produce  two  discemable  pulses,  due  to  bubble  fission  or  the  interaction  of 
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jet  impingement  with  the  wall.  But,  the  temporal  response  of  the  hydrophone  was  always  that  of 
a  single-degree-of -freedom  oscillator  responding  to  an  impulsive  pressure  input. 

The  measurement  illustrates  that  the  high-frequency  limit  of  the  measured  sound  spectrum  is 
limited  by  the  hydrophone  response,  even  though  the  hydrophone  would  be  considered 
“broadband”  relative  to  other  commercially  available  hydrophones.  The  pressure  pulse  created 
by  the  bubble  may  be  approximated  as P^(r,t)  =  as  discussed  above.  Pulse  widths 

measured  with  optical  and  piezoelectric  transducers  are  on  the  order  of  10  to  100  ns  [18,19].  An 
oscillator  will  respond  to  such  a  pulse  with  an  under-damped  response  like  that  of  Figure  3c. 
The  peak  amplitude  of  the  response  will  be  proportional  to  Pj,  and  the  initial  pulse  width  will  be 
proportional  tol/tir„,  the  natural  frequency  of  the  hydrophone,  if  l/ur^  »  0.  The  measured  pulse 
widths  were  typically  10  ±  2  ps  =  l/iiT„ ,  which  is  consistent  with  the  reported  hydrophone  natural 
frequency. 

The  actual  acoustic  impulse  will  be  approximately  /(r^) »  Psifn^i’^n)  ihe  location  of  the 
hydrophone,  r„.  The  measured  acoustic  impulse,  /„ ,  will  scale  as 

(7) 

based  on  the  dynamic  impulse  response  function  of  the  hydrophone.  However,  the  measured 
peak  pressure  and  the  pulse  duration  will  be  approximated  by 

respectively.  The  maximum  bubble  volume  accelerations  achieved  and  the  resulting  amplitude 
and  time-scale  of  the  emitted  pressure  wave  is  affected  by  the  amount  of  non-condensable  gas  in 
the  cavity,  the  collapsing  pressure,  and  the  topology  of  the  collapsing  bubble. 
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A  scaling  for  the  acoustic  impulse  can  be  derived  from  on  spherical  bubble  dynamics  and 
energy  conservation  [10].  However,  for  highly  deformed  bubbles,  it  may  be  more  appropriate  to 
scale  the  impulse  based  on  the  initial  energy  of  the  collapsing  bubble,  E^, 

£:^=KAP  (8) 

The  energy  radiated  away  from  the  bubble  after  the  collapse  and  rebound,  is  given  by 

— (9) 

pc 

where  rj^  is  the  proportion  of  mechanical  potential  energy  that  is  converted  into  acoustic  energy. 
Vogel  and  Lauterbom  [17]  have  shown  that  up  to  90%  of  a  bubble’s  mechanical  energy  can  be 
converted  to  acoustic  energy  during  the  first  collapse  of  spherical  laser-produced  cavitation 
bubbles,  but  the  percentage  can  be  reduced  substantially  if  the  collapsing  bubble  is  non- 
spherical. 

The  measured  peak  shock  pressure,  P5,  is  expected  to  scale  with  the  maximum  potential 
energy: 


_L(pcE^\'\_L(lhP£iM.r  (10) 
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The  acoustic  impulse  would  then  scale  as 
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The  functional  relationship  between  6  and  the  other  parameters  of  the  flow  are  not 
straightforward,  especially  in  the  case  of  non-spherical  collapse.  The  amount  of  non¬ 
condensable  gas  within  the  laser-produced  bubble  will  vary  due  to  variation  in  the  process  of 
bubble  creation.  Moreover,  the  amount  of  non-condensable  gas  at  collapse  will  increase  with 
increased  maximum  bubble  volume,  due  to  gas  diffusion,  and  the  minimum  bubble  volume  and 
shape  during  the  last  phases  of  collapse  will  also  vary.  Consequently,  variability  of  the  measured 
impulse  was  observed  given  nominally  constant  laser-energy  pulse  input.  However,  the 
relationship  between  the  impulse  and  the  maximum  bubble  volume  can  be  discerned. 


III.  RESULTS-  VORTEX  CAVITATION 

Next,  the  growth  and  collapse  of  non-spherical  hydrodynamic  cavitation  is  considered. 
Cavitation  bubbles  were  produced  in  the  vicinity  of  a  strong  line  vortex.  The  free-stream 
velocity  was  fixed  at  10  m/s.  Figure  4a  shows  a  vector  map  of  the  average  planar  vortical  flow 
field  measure  perpendicular  to  the  vortex  axis  produced  after  averaging  100  instantaneous  flow 
fields.  Figure  4b  plots  the  average  tangential  velocity,  Ugir),  as  a  function  of  the  distance  from 
the  vortex  center,  r,  along  with  a  fitted  curve  for  a  Gaussian  vortex: 
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Here,  a  « 1.255  which  make  r^.  the  “core  radius”  defined  as  the  radius  where  the  tangential 
velocity  is  maximum; 


Ue(rc)  =  P 


Io_ 

2m-c 


(13) 


where  )8  =  0.715.  For  the  measured  vortex,  the  fitted  core  radius  was  5.6  mm  and  the  strength 
was  0.290  m^/sec.  Examination  of  the  instantaneous  images  indicated  that  the  vortex  core  did 
not  wander  significantly.  The  amplitude  of  wandering  was  consistently  less  than  15%  of  the  core 
radius.  Consequently,  no  correction  for  vortex  wandering  was  needed  [21]. 

The  pressure  depression  due  to  the  vortex  is  given  by 


J  ^ 


f— 


_g-2a(r/rc)^ 

-2a{r  !  Ei{a(r  / 

^2a{r  /  Ei{2a{r  !  ref 


(14) 


where  Ei(x)  is  the  exponential  integral  function.  The  calculated  radial  pressure  distribution  for 
the  vortex  is  also  shown  in  Figure  4b.  Here  we  are  assuming  that  du^ldr^Q,  where  z 
corresponds  to  the  direction  along  the  vortex  axis.  Time  exposed  images  of  these  bubbles  in  the 
vortex  core  indicate  an  axial  velocity  within  95%  of  the  free-stream  speed.  The  pressure  at  the 
axis  of  the  vortex,  r  =  0,  is  given  by 
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(15) 


where  r}  =  0.870. 

Bubbles  captured  by  the  vortex  can  cavitate  if  the  core  pressure  is  at  or  below  vapor  pressure. 
The  free-stream  cavitation  number  is  defined  as  a„  =  (P„  -  Py)ljpUl ,  where  U„  is  the  free- 
stream  velocity.  Captured  bubbles  can  grow  if  P^.  s  Py,  making  the  cavitation  number  at 
inception,  cr„,, 


<7, 


(16) 


For  the  measured  vortex,  this  corresponds  to  a  free-stream  cavitation  number  of  1.18. 

Before  creating  bubbles  within  the  vortex,  near-spherical  bubbles  were  created  in  the  test 
section  but  without  any  flow.  Figure  5a  shows  ten  superimposed  images  of  the  laser  produced 
bubble  taken  with  a  time  interval  At  =  40  ps,  and  Figure  5b  shows  the  resulting  acoustic  pulse. 
The  test-section  static  pressure  was  31.1  kPa.  Again,  when  the  laser  pulse  initiates  the  bubble,  an 
acoustic  pulse  is  created.  The  bubble  passes  through  its  maximum  volume  and  collapses, 
creating  a  second  pulse  of  sound.  The  sound  pulse  created  upon  collapse  is  still  limited  by  the 
response  of  the  hydrophone.  Numerous  bubbles  were  created  with  fixed  static  pressure.  The 
maximum  size  of  the  bubbles  varied  with  the  amount  of  laser-light  energy  deposited  in  the 
liquid.  The  laser  pulse  energy  could  be  varied  manually,  but  the  actual  amount  of  light  absorbed 
by  the  liquid  varied  with  the  size  and  density  of  small  particulates  that  were  present  in  the  focal 
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volume  of  the  laser  during  the  pulse.  This  is  the  principal  cause  for  the  variability  of  the  laser- 
induced  bubbles. 

The  acoustic  pulse  created  during  the  collapse  of  the  bubble  was  recorded,  and  the 
acoustic  impulse  was  calculated  by  time  integration  of  the  pressure  pulse  using  the  trapezoidal 
rule.  The  maximum  bubble  diameter  was  recorded  with  the  digital  camera.  Figure  5c  shows  the 
variation  of  the  acoustic  impulse  with  the  maximum  bubble  radius.  While  there  is  significant 
scatter,  the  data  indicate  that  the  impulse  increases  with  the  bubble  volume,  as  expected.  The 
spectrum  analyzer  was  triggered  to  capture  the  noise  signal  produced  by  the  collapse,  although 
the  residual  echo  from  the  initiation  pulse  was  also  present.  To  calculate  the  impulse,  the 
average  reverberant  pattern  was  subtracted  from  the  individual  noise  traces,  and  then  the  collapse 
pulse  was  identified  and  measured.  This  reverberant  pattern  from  the  bubble  initiation  pulse  was 
highly  repeatable  when  the  bubble  and  hydrophone  locations  were  fixed. 

Figures  6, 7,  and  8  present  similar  series  of  data  for  bubbles  induced  in  the  center  of  the 
vortex.  Bubbles  were  created  with  three  free-stream  cavitation  numbers:  =  2.15,  1.73,  and 

1.49.  Wither^  =  2.15,  the  bubbles  collapse  immediately  after  they  are  created,  with  =  1.73, 
the  bubbles  slowly  collapse  after  growing  along  the  vortex  axis,  and  with  =  1.49,  the  bubbles 
do  not  collapse  but  continue  to  grow  along  the  axis. 

The  normalized  acoustic  impulse,  /,  is  defined  as: 


17 


OWEIS,  JASA 


While  the  values  of  r]^  and  6  are  unknown,  the  pulse  width,  6,  is  expected  to  be  on  the  order  of 
10  ns,  for  the  case  of  spherical  collapse  [18,19].  The  energy  conversion  efficiency,  rj^,  can  vary 
with  the  collapse  geometry  of  the  bubble  [18],  but  is  expected  to  be  >  60%  for  the  bubbles  under 

consideration.  Figure  9  shows  a  plot  of  the  normalized  /  against  the  bubble  maximum  radius 
with  r]^  =  1  and  0  =  10  ns  for  all  the  cases  studied.  The  normalized  acoustic  impulse  of  bubbles 
imploding  in  quiescent  fluid  as  well  as  bubbles  that  are  not  greatly  deformed  by  the  vortex  field 
(  a„  =  2.82  and  o„  =  2.15)  are  independent  of  the  bubble  maximum  radius,  as  expected  from  a 
proper  scaling.  The  deformed  bubbles  (  cr„  =  1.73)  have  a  normalized  impulse  that  is  clearly 
lower  than  the  spherical  and  less  deformed  bubbles.  If  it  is  assumed  here  that  the  collapse 
pressure  pulse  width  does  not  change  with  the  bubble  collapse  geometry.  The  results  then  imply 
that  highly  deformed  bubbles  are  less  efficient  (have  lower  rj^)  in  converting  their  mechanical 
energy  into  radiated  acoustic  pressure  in  comparison  to  the  slightly  deformed  and  spherically 
imploding  bubbles.  Still  in  all  of  these  highly  deformed  as  well  as  spherical  bubbles,  resolution 
of  the  collapse  pressure  pulse  is  limited  by  the  hydrophone  response.  The  measured  bubble 
collapse  time  (time  from  maximum  volume  to  minimum  volume)  relative  to  the  Rayleigh 
collapse  is  plotted  in  figure  10.  This  figure  shows  that  highly  deformed  bubbles  ((  =  1.73) 

have  a  longer  collapse  time  than  the  more  acoustic-efficient  spherical  bubbles,  but  the  collapse 
time  is  on  the  same  order. 

A  second  case  was  considered  where  the  bubbles  did  not  collapse  but  continued  to  grow 
at  C7„  =  1.19,  and  this  is  shown  in  figure  11.  Here,  the  laser  energy  was  greatly  reduced,  created 
only  a  small  bubble  that  grew  slowly  under  the  influence  of  the  low-pressure  in  the  vortex  core. 
In  this  case,  the  nucleus  grew  slowly  and  no  noise  was  detectable.  The  magnitude  of  the  volume 
acceleration  during  the  natural  bubble  growth  is  much  smaller  that  the  acceleration  during  the 
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last  stages  of  bubble  collapse.  As  the  bubble  grows  under  the  influence  of  the  low  vortex  core 
pressure,  R  ~  -^J-AP/p,  and  R~0.  Consequently,  the  minimal  amount  of  noise  produced  during 
the  slow  process  of  bubble  growth  was  not  detectable  over  the  background  noise  created  by  the 
water  tunnel  flow. 

Lastly,  the  pressure  spectra  recorded  by  the  hydrophone  for  many  bubble  collapses  was 
collected  for  the  cases  described  above.  The  sound  traces  for  1000  events  were  collected.  These 
are  shown  in  figure  12.  These  spectra  do  not  include  the  noise  pulse  produced  by  the  initial 
bubble  growth,  but  do  include  some  of  the  reverberation  resulting  from  bubble  production.  The 
case  of  the  small  bubble  production  (sb)  corresponds  to  the  background  noise  of  the  tunnel 
operating  at  10  m/s.  The  spectral  features  below  30  kHz  correspond  to  reverberant  modes  of  the 
facility.  As  the  impulse  of  the  bubbles  increases,  the  energy  in  the  reverberations  increases  as 
well.  Also,  the  energy  in  frequencies  greater  than  30  kHz  increases,  but  the  spectra  is  flat  until 
100  kHz.  This  is  consistent  with  the  input  of  a  sharp  impulse  into  the  bandwidth-limited 
hydrophone. 


IV.  DISCUSSION  AND  CONCLUSIONS 

Acoustic  transients  by  discreet  hydrodynamic  cavitation  events  were  examined  in  both  a 
quiescent  and  a  vortical  flow,  and  the  following  conclusions  were  found: 

The  hydrophone  used  was  not  capable  of  resolving  the  dynamics  of  the  pressure  pulse 
produced  upon  the  collapse  of  a  spherical  cavitation  bubble. 

The  acoustic  impulse  produced  by  a  collapsing  bubble  scales  with  the  collapse  pressure 
and  the  maximum  volume,  even  for  non-spherical  collapse.  The  noise  pulse  produced  by 
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deformed  collapsing  bubbles  is  often  sufficiently  narrow  to  be  unresolved  by  the 
hydrophone. 

The  normalized  impulse  created  by  the  collapse  of  spherical  bubbles  represents  the  most 
efficient  collapse,  converting  potential  mechanical  energy  into  acoustic  energy. 
Deformed  bubbles  collapse  with  less  transfer  efficiency. 

The  noise  produced  by  the  natural  bubble  growth  in  the  vortex  is  not  measurable. 

Vortex  cavitation  occurs  when  naturally  occurring  nuclei  are  captured  by  a  concentrated 
vortex,  and  they  grow  and  collapse.  The  results  of  this  study  suggest  that  significant  noise  is 
created  only  when  the  bubbles  collapse,  even  if  they  are  highly  distorted.  The  noise  produced  is 
broadband,  with  a  high-frequency  cut-off  beyond  the  bandwidth  of  most  conventional 
hydrophones.  However,  the  impulse  created  by  the  collapsing  bubbles  can  be  scaled  with 
spherical  bubble  dynamics.  Such  scaling  will  provide  the  upper  bound  on  the  noise  that  a  bubble 
of  a  given  maximum  volume  can  create  when  collapsing  under  a  given  pressure  difference. 

Finally,  an  unsteady  pressure  in  the  vortex  core  is  required  for  a  single  nucleus  to  experience 
both  rapid  growth  and  collapse.  The  vortical  flow  employed  in  the  present  study  is  steady,  and 
natural  cavitation  occurred  only  when  the  core  pressure  dropped  below  vapor  pressure.  At  this 
point,  captured  nuclei  continued  to  grow  and  ultimately  filled  the  vortex  core  as  the  bubble 
convected  downstream.  The  case  of  growth  without  collapse  can  be  contrasted  with  cavitation 
accompanying  more  complicated  vortical  flows.  In  these  flows  variation  in  the  local  mean 
pressure,  vortex-vortex  interactions,  and  other  flow  unsteadiness  can  result  in  rapid  growth  and 
collapse  of  bubbles  in  the  vortex  core.  Such  cavitation  occurs  in  turbulent  shear  flows  and  in  the 
wake  vortex  system  of  complex  turbo-machinery.  However,  the  results  of  this  study  suggest 
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that  spherical  bubble  dynamics  can  be  used  to  scale  the  upper  bound  of  the  noise  impulse 
produced  by  discreet  cavitation  events,  even  when  the  bubbles  are  quite  deformed. 
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FIGURE  CAPTIONS 

Figure  1:  (a)  Schematic  diagram  of  the  quiescent  cavitation  bubble  experiment;  (b)  the  bubble- 
vortex  interaction  experiment.  A  single  laser  pulse  is  used  to  create  a  cavitation  bubble  in  the 
bulk  of  the  fluid.  Images  of  the  bubble  are  captured  with  a  multi-exposure  digital  camera,  and 
the  acoustic  emission  of  the  bubble  is  captured  with  a  hydrophone. 

Figure  2:  Photograph  of  the  hydrofoil  in  the  water  tunnel  test  section  with  developed  tip-vortex 
cavitation  with  =  10  m/s,  and  =  1.4.  Flow  is  from  left  to  right 

Figure  3:  (a)  Five  superposed  exposures  of  a  laser  produced  spherical  cavitation  bubble  taken 
every  6  ps  in  a  quiescent  water  container;  (b)  the  resulting  acoustic  pressure  trace  recorded  from 
the  hydrophone;  (c)  the  expanded  hydrophone  pressure  trace  produced  upon  bubble  collapse 
(solid  line)  along  with  the  second-order  response  of  a  SDOF  oscillator  to  a  Dirac-delta  input  with 
system  tUn  =  105  kHz  and  damping  ratio,  ^=0. 18  (dashed  line). 

Figure  4:  (a)  The  average  vector  map  of  the  flow  perpendicular  to  the  vortex  axis  90  mm 
downstream  of  the  tailing  edge  for  10  m/s;  (b)  the  measured  tangential  velocity,  Ug(r), 
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along  with  the  velocity  distribution  for  the  Gaussian  vortex  with  Vf-  =  5.6  mm  and  =  0.290 
mVs  and  the  resulting  calculated  pressure  depression,  P{r)  - 

Figure  5:  (a)  Ten  images  of  a  laser  produced  bubble  taken  40  ps  apart  in  the  water  tunnel  test 
section  without  flow  with  =  31  kPa;  image  dimensions  (V  23.1  x  H  31.9  mm);  (b)  the 
pressure-time  trace  recorded  by  the  hydrophone;  (c)  the  acoustic  impulse  created  during  bubble 
collapse  versus  the  maximum  bubble  radius. 

Figure  6:  (a)  Ten  exposures  of  a  laser  produced  bubble  created  in  the  vortex  taken  54  ps  apart 
with  cr„  =  2.15,  f/„=10m/s;  the  bright  region  on  the  right  results  from  pixel  saturation  from  the 
original  laser  pulse;  image  dimensions  (V  15.9  x  H  21.0  mm);  (b)  the  pressure-time  trace 
recorded  by  the  hydrophone;  (c)  the  acoustic  impulse  created  during  bubble  collapse  versus  the 
normalized  maximum  bubble  radius  (normalized  by  the  vortex  core  radius). 

Figure  7:  (a)Ten  exposures  of  a  laser  produced  bubble  created  in  the  vortex  taken  94  ps  apart 
with  a„  =  1.73,  (/„=10m/s;  the  bright  region  on  the  right  results  from  pixel  saturation  from  the 
original  laser  pulse;  image  dimensions  (V  15.9  x  H  21.0  mm);  (b)  the  pressure -time  trace 
recorded  by  the  hydrophone;  (c)  the  acoustic  impulse  created  during  bubble  collapse  versus  the 
normalized  maximum  bubble  radius  (normalized  by  the  vortex  core  radius). 

Figure  8:  (a)  Ten  exposures  (last  two  exposures  outside  the  frame)  of  a  laser  produced  bubble 
created  in  the  vortex  taken  154  ps  apart  with  =  1.49,  (7„=10m/s;  the  bright  region  on  the 
right  results  from  pixel  saturation  from  the  original  laser  pulse;  image  dimensions  (V  15.9  x  H 
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21.0  mm);  (b)  the  pressure-time  trace  recorded  by  the  hydrophone.  At  this  cavitation  number, 
the  bubble  continues  to  grow  along  the  vortex  axis  after  initiation  with  the  laser. 

Figure  9:  The  dimensionless  acoustic  impulse,  I  =  I.rff  l(^r]^pc  AP  ,  versus  the 

dimensionless  maximum  bubble  radius,  =  R^,  tr^.,  for  P„=  10  kPa  without  flow  {•),?„  =  31 

kPa  without  flow  (*),  cr„  =  2.82  (+),  =  2.15  (o),  and  a„  =  1.73  (<).  For  this  plot  the 

acoustic  conversion  efficiency  was  chosen  to  be  100%,  and  the  collapse  pulse  width  Q  10 
nanoseconds. 

Figure  10:  The  dimensionless  collapse  time  tc/TR^yieigh  versus  the  dimensionless  maximum  radius, 
R^  =Ri^  I r^,  for  P„=  10  kPa  without  flow  (•),  Poo=  3 1  kPa  without  flow  ( * ),  a„  =  2.82  (+),  a„  = 
2.15(0),  and  a„  =  1.73(<). 

Figure  11:  (a)  Ten  exposures  of  a  bubble  created  by  a  low  energy  laser  pulse  in  the  vortex  taken 
300  ps  apart  with  cr„  =  1.19;  the  initial  nucleus  is  small  and  not  visualized  in  the  image;  flow  is 
from  right  to  left;  (b)  the  pressure-time  series  recorded  by  the  hydrophone,  which  is  equivalent 
to  the  baseline  hydrophone  signal  before  inception. 

Figure  12:  Average  power  spectral  density  of  the  hydrophone  signal  for  the  5  different  cases 
corresponding  to  macroscopic  bubbles  created  without  test-section  flow,  and  in  the  vortex  for 
=  2.8,  1.70.  1.28  The  spectrum  of  the  noise  recorded  when  the  microscopic  bubbles  (sb) 
were  created  at  =  1.19  is  also  included,  and  this  is  equivalent  to  the  baseline  spectrum 
without  cavitation.  The  spectrum  analyzer  was  triggered  to  analyze  the  noise  signal  produced 
due  to  the  bubble  collapse  only,  and  the  spectra  from  1000  acoustic  events  were  averaged. 
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Figure  1:  (a)  Schematic  diagram  of  the  quiescent  cavitation  bubble  experiment;  (b)  the  bubble-vortex  interaction  experiment.  A  single 
laser  pulse  is  used  to  create  a  cavitation  bubble  in  the  bulk  of  the  fluid.  Images  of  the  bubble  are  captured  with  a  multi-exposure 
digital  camera,  and  the  acoustic  emission  of  the  bubble  is  captured  with  a  hydrophone. 
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Figure  2;  Photograph  of  the  hydrofoil  in  the  water  tunnel  test  section  with  developed  tip-vortex  cavitation  with  =  10  m/s,  and 
=  1.4.  Flow  is  from  left  to  right 
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(a) 


Figure  3:  (a)  Five  superposed  exposures  of  a  laser  produced  spherical  cavitation  bubble  taken  every  6  p.s  in  a  quiescent  water 
container;  (b)  the  resulting  acoustic  pressure  trace  recorded  from  the  hydrophone;  (c)  the  expanded  hydrophone  pressure  trace 
produced  upon  bubble  collapse  (solid  line)  along  with  the  second-order  response  of  a  SDOF  oscillator  to  a  Dirac-delta  input  with 
system  xun  =  105  kHz  and  damping  ratio,  ^=0,18  (dashed  line). 
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Figure  4:  (a)  The  average  vector  map  of  the  flow  perpendicular  to  the  vortex  axis  90  mm  downstream  of  the  tailing  edge  for  f/„=  10 
m/s;  (b)  the  measured  tangential  velocity,  Ma(r),  along  with  the  velocity  distribution  for  the  Gaussian  vortex  with  =  5.6  mm  and 
Tq  =  0.290  m'/s  and  the  resulting  calculated  pressure  depression,  P(r)  -  P„. 
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Figure  5:  (a)  Ten  images  of  a  laser  produced  bubble  taken  40  ps  apart  in  the  water  tunnel  test  section  without  flow  with  =  31  kPa; 
image  dimensions  (V  23.1  x  H  31.9  mm);  (b)  the  pressure-time  trace  recorded  by  the  hydrophone;  (c)  the  acoustic  impulse  created 
during  bubble  collapse  versus  the  maximum  bubble  radius. 


OWEIS,  JASA 


Tlmo,  |i* 


(b) 


(c) 

Figure  6:  (a)  Ten  exposures  of  a  laser  produced  bubble  created  in  the  vortex  taken  54  \is  apart  with  =  2.15,  {7„=10m/s;  the 
bright  region  on  the  right  results  from  pixel  saturation  from  the  original  laser  pulse;  image  dimensions  (V  15.9  x  H  21.0  mm);  (b)  the 
pressure-time  trace  recorded  by  the  hydrophone;  (c)  the  acoustic  impulse  created  during  bubble  collapse  versus  the  normalized 
maximum  bubble  radius  (normalized  by  the  vortex  core  radius). 
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Figure  7;  (a)Ten  exposures  of  a  laser  produced  bubble  created  in  the  vortex  taken  94  ^s  apart  with  -  1.73,  =10m/s;  the  bright 

region  on  the  right  results  from  pixel  saturation  from  the  original  laser  pulse;  image  dimensions  (V  15.9  x  H  21.0  mm);  (b)  the 
pressure-time  trace  recorded  by  the  hydrophone;  (c)  the  acoustic  impulse  created  during  bubble  collapse  versus  the  normalized 
maximum  bubble  radius  (normalized  by  the  vortex  core  radius). 
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Figure  8:  (a)  Ten  exposures  (last  two  exposures  outside  the  frame)  of  a  laser  produced  bubble  created  in  the  vortex  taken  154  ps  apart 
with  =  1.49,  f/„=10m/s;  the  bright  region  on  the  right  results  from  pixel  saturation  from  the  original  laser  pulse;  image 
dimensions  (V  15.9  x  H  21.0  mm);  (b)  the  pressure-time  trace  recorded  by  the  hydrophone.  At  this  cavitation  number,  the  bubble 
continues  to  grow  along  the  vortex  axis  after  initiation  with  the  laser. 
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Figure  11:  (a)  Ten  exposures  of  a  bubble  created  by  a  low  energy  laser  pulse  in  the  vortex  taken  300  fxs  apart  with  =  1.19;  the 
initial  nucleus  is  small  and  not  visualized  in  the  image;  flow  is  from  right  to  left;  (b)  the  pressure-time  series  recorded  by  the 
hydrophone,  which  is  equivalent  to  the  baseline  hydrophone  signal  before  inception. 
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Figure  12;  Average  power  spectral  density  of  the  hydrophone  signal  for  the  5  different  cases  corresponding  to  macroscopic  bubbles 
created  without  test-section  flow,  and  in  the  vortex  for  C7^  =  2.8,  1.70.  1.28  The  spectrum  of  the  noise  recorded  when  the 
microscopic  bubbles  (sb)  were  created  at  (7^^  =  1.19  is  also  included,  and  this  is  equivalent  to  the  baseline  spectrum  without 
cavitation.  The  spectrum  analyzer  was  triggered  to  analyze  the  noise  signal  produced  due  to  the  bubble  collapse  only,  and  the  spectra 
from  1000  acoustic  events  were  averaged. 
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Abstract 

The  capture  of  a  bubble  by  a  line  vortex  for  both  non-cavitating  and  cavitating  conditions  using  a 
one-way  Particle-Tracking  Model  (PTM)  and  Direct  Numerical  Simulation  (DNS)  is  examined. 
These  results  are  compared  to  experimental  observations.  The  particle-tracking  model  can 
successfully  predict  the  capture  time  for  small  bubbles  by  a  line  vortex.  For  the  small  bubbles 
starting  far  from  the  vortex,  variation  in  capture  time  with  differing  choices  of  the  lift  coefficient 
are  relatively  small.  The  accuracy  of  the  PTM  is  reduced  as  the  bubble  moves  to  within  the  core 
radius,  and  the  DNS  solutions  suggest  that  the  PTM  will  under-predict  the  capture  time  in  the 
final  stages  of  collapse.  The  effect  of  bubble  deformation  and  bubble  vortex  interactions  will  be 
important  only  at  the  final  stages  of  bubble  capture.  The  non-cavitating  PTM  can  be  used  when 
considering  the  capture  of  cavitation  nuclei.  The  bubble  volume  changes  very  slowly  until  the 
final  stages  of  bubble  capture.  Consequently,  it  is  useful  to  separate  the  process  into  the  non- 
cavitating  capture  of  the  nucleus  and  then  the  growth  of  the  nucleus  in  the  low-pressure  core 
region.  The  growth  of  the  nucleus  prior  to  cavitation  inception  is  relatively  small,  but  such 
bubble  growth  will  lead  to  a  somewhat  faster  capture  as  the  bubble  comes  near  the  core. 
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1.  Background 

Hydrodynamic  lifting  surfaces  associated  with  propulsion  and  control  systems  often 
develop  concentrated  trailing  vortices.  Small  bubbles  near  the  concentrated  vortex  will  be 
drown  toward  to  axis  of  flow  rotation,  and  if  the  core  pressure  falls  below  vapor  pressure,  these 
small  bubbles  can  act  as  nucleation  sites  for  vortex  cavitation,  and  Arndt  (2002)  has  recently  this 
phenomenon.  The  probability  that  a  cavitation  bubble  will  occur  is  related  the  concentration  and 
size  distribution  of  the  free-stream  nuclei  {i.e.  the  nuclei  number  distribution)  and  the  likelihood 
that  the  bubble  is  found  in  the  possibly  transient  region  of  low  pressure.  The  bubble  may  already 
be  present  in  the  core  region  when  an  unsteady  flow  phenomenon,  such  as  vortex  stretching, 
reduced  the  core  pressure.  Or,  the  nucleus  could  be  captured  by  the  low-pressure  vortex.  The 
likelihood  of  bubble  capture  is  related  to  the  probability  of  the  nuclei’s  presence  within  or  near 
the  vortex  core  and  to  the  probability  that  nuclei  may  be  captured  by  the  vortex  while  the  core 
pressure  in  the  vortex  is  reduced.  This  second  process  is  governed  by  the  time  a  nucleus  takes  to 
move  close  to  the  vortex  axis,  the  capture  time.  With  an  estimate  of  the  capture  time,  it  is  then 
possible  to  identify  which  nuclei  surrounding  an  unsteady  vortex  will  lead  to  discreet  vortex 
cavitation  bubbles. 

The  vortex-induced  capture,  interaction,  and  growth  of  small  bubbles  have  been  studied 
by  several  researchers.  Sridhar  and  Katz  (1999)  studied  the  effect  of  an  entrained  bubble  in  a 
vortex  ring  using  a  particle  image  velocimetry  and  showed  that  in  certain  cases  the  entrained 
bubble  can  substantially  change  the  structure  of  the  vortex.  Hsiao  et  al.  (1999)  examined  the 
capture  of  a  cavitation  nucleus  by  a  tip  vortex  employing  a  Reynolds  Averaged  Navier  Stokes 
(RANS)  model  of  the  non-cavitating  flow  along  with  a  particle-tracking  model  to  describe  the 
dynamics  of  the  nucleus.  This  approach  assumes  that  the  dynamics  of  the  bubbles  can  be 


predicted  with  spherical  bubbles  dynamics  and  a  one-way  coupled  momentum  balance  on  the 
bubble.  Chahine  and  co-workers  have  used  using  both  RANS  and  boundary  integral  methods  to 
study  vortex  bubble  interaction  (Chahine,  1990;  Chahine  and  Hsaio,  1994;  Hsiao  and  Chahine 
2001;  Chahine  and  Hsaio,  2002).  They  have  shown  how  captured  bubbles  can  significantly 
deformed  as  they  are  captured  by  the  vortex,  grow,  and  collapse. 

Models  of  bubble-flow  interactions  often  employ  one-way  coupling  between  the  flow  and 
the  motion  of  the  bubble  whereby  the  presence  of  the  bubble  does  not  alter  the  underlying  flow 
field.  An  equation  of  motion  for  the  trajectory  of  the  bubble  is  solved  using  a  force  balance  on 
the  bubble.  For  cavitating  flows,  the  Rayleigh-Plesset  equation  is  employed  to  describe  the 
dynamics  of  the  cavitating  bubble.  Such  Particle  Tracking  Models  (PTM)  often  assumes  that  the 
bubbles  remain  spherical  during  the  bubble  capture,  growth,  and  collapse.  However,  as  the 
volume  of  a  bubble  increases,  the  pressure  and  flow  gradients  near  the  vortex  can  cause  the 
bubble  to  deform  and  can  lead  to  significant  modification  of  the  forces  on  the  bubble.  Higher 
fidelity  simulations  of  the  bubble-flow  interactions  are  possible,  but  can  be  computationally 
intensive.  It  would  be  preferable  if  the  one-way  particle  tracking  models  can  be  used  to  yield  an 
accurate  estimate  of  the  bubble  capture  times. 

In  the  present  work,  we  will  examine  the  capture  of  a  bubble  by  a  steady  line  vortex  for 
both  non-cavitating  and  cavitating  conditions  using  a  one-way  particle-tracking  model  and  Direct 
Numerical  Simulation  (DNS).  Experimental  observations  of  bubble  capture  are  then  compared 
with  the  predicted  results  of  the  PTM.  These  data  are  used  to  explore  the  accuracy  and 
limitations  of  the  one-way  particle- tracking  models. 


2.  One-Way  Coupled  Particle  Tracking  Model  of  Bubble  Capture 

The  dynamics  of  small  bubbles  in  the  flow  can  be  calculated  assuming  that  their  presence 
does  not  significantly  alter  the  non-cavitating  flow.  The  effect  of  the  flow  on  a  bubble  can  be 
determined  through  the  equation  of  motion  for  the  particle  (Maxey  and  Riley  (1983);  Magnaudet 
and  Eames  (2000)): 

P.V,^.V,(p.-p,)g-V.Vp  +  lp,V',(^-i^)  +  F„+r,  +  F.*F,  (1) 

where  Pg  is  the  bubble  density,  is  the  fluid  density,  is  the  volume  of  the  bubble,  is  the 
bubble  velocity,  u  is  the  fluid  velocity  in  the  absence  of  the  bubble  at  the  center  of  the  bubble 
(the  unperturbed  velocity).  is  the  drag  force,  is  the  lift  force,  F^  is  the  Basset  history 
force,  and  Fj,  is  the  force  due  to  volume  changes.  In  the  present  analysis,  the  drag  force  is  given 

by 

FD=C'oip^AB|u-UB|(u-U3)  (2) 

where 

Co  =—(l  +  0.197Re/-'' +  2.6(10-' )Re/'*)  (3) 

and  Re^  =|u-Uo|2/?/v^  is  the  bubble  Reynolds  number  based  on  the  bubble  radius,  R,  and 
projected  area,  Ag  (Haberman  and  Morton  (1953)).  We  will  neglect  the  Basset  force  term,  as  the 
bubble  relative  acceleration  is  expected  to  be  small  (a  detailed  discussion  of  the  relative  size  of 
the  Basset  term  is  found  in  Maxey  and  Riley  (1983)).  The  force  due  to  the  bubble  volume 
variation  is  given  by  Johnson  and  Hsieh  (1966)  as: 


The  buoyancy,  pressure  gradient,  and  added  mass  terms  are  also  included  in  equation,  with  an 
added  mass  coefficient  of  1/2  for  spherical  bubbles.  We  will  ignore  the  effects  of  buoyancy  on 
the  small  bubbles  considered  here. 

The  lift  forces  on  the  bubble  is  related  to  the  slip  velocity  and  the  circulation  of  the  un¬ 
perturbed  fluid  at  the  center  of  the  bubble: 

°  Q  T  |(^  -  )  =  f  ~  ^ 

where  a  =  |£u|i?/|u-Ug|  is  the  dimensionless  shear  rate,  and  O)  is  the  fluid  vorticity  of  the 
undisturbed  flow  at  the  center  of  the  bubble.  The  lift  coefficient  is  a  function  of  both  the  shear 
and  vorticity.  Auton  (1987)  and  Auton  et  al.  (1988)  showed  that  for  weak  shear  (a  « 1)  the  lift 
coefficient  is  given  by 

“s'® 

Saffman  (1965)  showed  that,  for  low  Reynolds  numbers  (Re^  <  20), 

Q^=5.82Re3‘°-^a‘’^  (7) 

For  higher  Reynolds  numbers.  Dandy  and  Dwyer  (1990)  used  numerical  analysis  to  show  that 
for  0.0005  <a<  0.4  and  Re^  >  20 , 

(8) 

Lastly,  Sridar  and  Katz  (1995)  experimentally  determined  the  lift  coefficient  on  bubbled  with 
diameters  between  500  and  800  p.m  and  20  <  Re^  <  80,  and  found  relatively  high  lift  coefficients 
Qj;,  =0.59a°^^  (9) 

To  determine  the  volume  variation  for  a  spherical  bubble,  the  Rayleigh-Plesset  equation 
given  is  used  to  determine  the  time- varying  bubble  radius,  R{t) : 


(10) 


dt  2\dt 
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Pv  +  Pgi\  Pb  ^  ^ 


This  relation  is  derived  for  a  spherical  bubble  in  an  unbounded  fluid,  and  pg  is  the  pressure  far 
away  from  the  bubble,  p^i  is  the  pressure  of  the  non-condensable  gas  initially  in  the  bubble,  k  is 
the  polytropic  exponent,  and  a  is  the  surface  tension.  For  the  particle  tracking  model,  pg  is 
often  taken  as  the  local  fluid  static  pressure  at  the  location  of  the  bubble  center.  If  the  bubble  is 
small  and  the  pressure  gradients  near  the  bubble  are  mild,  then  Vp  -  R!  p  «\  at  the  bubble 

center,  and  we  may  model  Pg  ^  p- 


3.  Gaussian  Vortex  Model 

We  will  consider  the  capture  of  bubbles  by  a  Gaussian  vortex.  Also  known  as  Lamb- 
Oseen  vortex,  the  vorticity  distribution  at  a  radius,  r  and  at  time,  t  is  given  by: 
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The  tangential  velocity  around  the  axis  is  given  by: 
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In  many  cases  of  high  Reynolds  number  flows,  the  time-scale  associated  with  vorticity  diffusion, 
»  rj'  lAVf,  is  long  compared  to  the  time  of  bubble  capture.  We  can  then  assume  that  the 
vorticity  distribution  is  constant,  with  a  resulting  velocity  distribution 


(13) 


where  =1.255,  and  the  maximum  tangential  velocity,  u^,  occurs  at  the  core  radius  and  is 
given  by 
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where  172  =  0.715.  The  pressure  depression  at  the  vortex  centerline,  r  =  0,  is  given  by 
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At  the  vortex  core, 
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where  773  =  /(O)  =  0.870.  Here  we  are  neglecting  the  possibility  that  there  may  be  non-uniform 


flow  along  the  vortex  axis. 

The  vortex  Reynolds  number  is 

Rec=-^  (17) 

The  Weber  number  of  a  bubble  with  radius  Rq  is 

(18) 
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The  vortex  cavitation  number  is  defined  as: 
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4.  Direct  Numerical  Simulations  of  Bubble  Dynamics 

The  numerical  method  used  in  this  study  to  compute  the  Navier-Stokes  equations  is  a 
Direct  Numerical  Simulation  (DNS)  front-  tracking  method  for  multiphase  flows  developed  by 
Unverdi  and  Tryggvason  (1992)  and  improved  by  Esmaeeli  and  Tryggvason  (1998,  1999).  The 
fluids  inside  and  outside  the  bubbles  are  taken  to  be  Newtonian  and  the  flow  is  incompressible 
and  isothermal,  so  that  the  densities  and  viscosities  are  constant  within  each  phase.  The 
unsteady,  viscous,  incompressible  Navier-Stokes  equations  are  solved  by  a  finite-difference 
approximation  on  a  fixed  staggered  grid,  while  the  interface  between  the  bubbles  and  fluid  were 
explicitly  represented  by  a  separate,  unstructured  grid  that  moves  through  the  stationary  one. 
The  front  (interface)  keeps  the  density  and  viscosity  discontinuities  sharp  and  also  provides 
means  for  calculation  of  the  surface  tension  forces.  A  numerical  method  based  on  the  Marker- 
And-Cell  method  is  used  to  solve  the  unsteady  Navier-Stokes  equations  over  the  entire 
computational  domain.  This  “one-field”  formulation,  which  computes  both  the  properties  within 
the  bubbles  and  in  the  ambient  fluid,  is  used  in  conjunction  with  a  linked  list  of  markers  that  is 
used  to  explicitly  track  the  position  of  the  fluid  interfaces.  A  single  Navier-  Stokes  equation  with 
variable  density  and  viscosity  is  solved  for  the  entire  domain,  and  the  surface  tension  is  added  as 
a  body  force  concentrated  at  the  fluid  interfaces.  The  numerical  scheme  used  is  a  conservative 
second-order  accurate,  centered  difference  scheme  for  spatial  variables,  and  explicit  second- 
order  time  integration.  Yu  et  al.  (1995)  used  this  front  tracking  method  to  examine  collapse  of 
bubbles  is  shear-flow,  and  the  model  developed  in  this  study  is  used  in  the  present  work. 
Resolution  and  accuracy  studies  for  the  cavitation  models  were  reported  in  that  study. 


5.  Experimental  Setup 

Figure  la  shows  the  setup  for  the  creation  of  single  cavitation  bubbles  in  a  vortical  flow. 
Single  cavitation  bubbles  were  created  near  a  steady  line  vortex  using  the  University  of  Michigan 
9-Inch  Cavitation  Tunnel.  The  water  tunnel  has  a  circular  contraction  downstream  of  a  series  of 
flow  management  screens  with  contraction  ratio  6.4:1.  The  test  section  has  a  22.9  cm  diameter 
round  inlet  that  is  then  faired  into  a  rectangular  test  section  with  widely  rotmded  comers.  Four 
acrylic  windows  (93.9  cm  by  10.0  cm  viewing  area)  permit  optical  access  to  the  test  section  flow. 
The  flow  in  the  test  section  can  be  operated  at  pressures  from  vapor  pressure  to  approximately 
200  kPa.  The  average  velocity  in  the  test  section  is  variable  up  to  18  m/s.  A  de-aeration  system 
can  be  used  to  vary  the  dissolved  gas  content  of  the  flow,  and  the  inlet  water  is  filtered  to  1 
microns. 

A  vortical  flow  was  created  using  a  cambered  hydrofoil  mounted  to  one  window  of  the 
test  section.  The  hydrofoil  has  a  rectangular  planform  of  9.5  cm  span  and  16.8  cm  chord.  The 
tip  of  the  hydrofoil  was  truncated  with  sharp  edges.  The  hydrofoil  mount  allows  continuous 
changes  of  the  incident  flow  angle.  A  series  of  tip  and  trailing  edge  vortices  will  be  shed  near  the 
tip,  and  these  vortices  will  merge  to  form  a  single  vortex  within  one-half  chord  length 
downstream  of  the  trailing  edge.  Measurements  of  the  bubble/vortex  interactions  were  taken  for 
a  free-stream  velocity  of  10  m/s  and  a  variety  of  pressures.  The  dissolved  oxygen  content  was 
measured  with  an  Orion  Model  810  dissolved  oxygen  meter.  In  order  to  reduce  the  number  of 
free-stream  nuclei,  the  free-stream  gas-content  was  reduced  to  below  1.5  ppm  during  the 
measurements.  A  photograph  of  the  hydrofoil  with  developed  tip  vortex  cavitation  is  shown  in 
Figure  lb. 


Planar  Particle  Imaging  Velocimetry  (PIV)  was  used  to  measure  the  vortical  flow  field  at 
a  station  9.0  cm  downstream  of  the  trailing  edge.  A  double-pulsed  light  sheet  9  mm  thick  was 
created  perpendicular  to  the  mean  flow  direction  using  two  pulsed  Nd:YAG  lasers  (Spectra 
Physics  model  Pro-250  Series).  15.3  micron  average  diameter  silver  coated  glass  spheres  were 
used  to  seed  the  flow.  An  acrylic  prism  was  optically  mounted  to  a  window  of  the  test-section 
for  viewing  of  the  light  sheet  with  reduced  optical  distortion.  The  flow  was  captured  with  a  PIV 
image  capture  system  produced  by  LaVision  Inc.  Double-pulsed  images  of  the  light  sheet  were 
acquired  with  a  digital  camera  with  1280  x  1024  pixels.  Optical  distortion  of  the  planar  light 
sheet  image  was  corrected  through  a  calibration  procedure  that  employed  the  imaging  of  a 
regular  grid  in  the  location  of  the  light  sheet  plane.  Velocity  vectors  were  produced  from  the 
double-pulsed  images  using  the  LaVision  image  analysis  software  Davis6.0.4.  Multipass 
processing  with  a  final  window  size  of  32  x  32  pixels  was  used  with  12%  window  overlap  in  the 
final  pass  to  produce  41  x  27  in  plane  velocity  vectors  at  0.62  mm  spacing.  Since  the  camera¬ 
imaging  plane  was  not  parallel  to  the  light  sheet  and  it  had  an  angular  shift  of  45°  with  the 
horizontal  direction,  this  velocity  component  was  corrected  by  knowledge  of  the  camera  angle 
and  the  tunnel  free  stream  velocity. 

For  these  experiments,  the  free-stream  velocity  was  fixed  at  10  m/s.  Figure  2a  shows  a 
vector  map  of  the  average  planar  vortical  flow  field  measure  perpendicular  to  the  vortex  axis 
produced  after  averaging  300  instantaneous  flow  fields.  Figure  2b  plots  the  average  tangential 
velocity,  Ug{r),  as  a  function  of  the  distance  from  the  vortex  center,  r,  along  with  a  fitted  curve 
for  a  Gaussian  vortex  and  the  computed  pressure  depression  as  a  result  of  the  vortex.  For  the 
measured  vortex,  the  fitted  core  radius  was  =5.6  mm  and  the  strength  was  Fq  =  0.290  m^/sec. 
Examination  of  the  instantaneous  images  indicated  that  the  vortex  core  did  not  wander 


significantly.  The  amplitude  of  wandering  was  consistently  less  than  10%  of  Consequently, 
no  correction  for  vortex  wandering  was  needed,  as  discussed  by  Boulon  et  al.  (1999). 

Optically  stimulated  vortex-cavitation  bubbles  were  created  at  the  vortex  centerline.  As 
the  cavitation  number  was  reduced,  discreet  cavitation  events  occurred  near  the  vortex  core  due 
to  the  creation  of  a  nucleus  after  the  collapse  of  the  initial  laser-induced  bubble.  Natural 
cavitation  also  occurred  are  various  locations  along  the  vortex  axis  downstream  of  the  location  of 
vortex  roll-up.  However,  at  such  low  air  content,  the  event  rate  of  this  cavitation  was  relatively 
low  (<  0.1  event  per  second).  Bubble-capture  and  any  resulting  vortex  cavitation  could  then  be 
examined  after  the  introduction  of  a  nucleus  though  a  focused  pulse  of  laser  light  near  the  axis  of 
the  vortex.  By  varying  the  fi'ee-stream  pressure  and  laser  energy,  the  initial  size  of  the  laser- 
induced  bubble  could  be  controlled.  Images  of  the  bubble  were  observed  with  a  Cooke  Corp. 
“FlashCam”  digital  camera.  This  camera  has  an  effective  resolution  of  750  x  240  pixels,  and  it  is 
capable  of  recording  up  to  consecutive  10  exposmes  on  the  same  image,  with  varying  integration 
times  and  delays  between  each  exposure.  A  pulse  delay  generator  (SRS  model  DG535)  was  used 
to  trigger  the  camera  at  varying  delays  from  the  laser  light  pulse. 

6.  Non-Cavitating  Bubble  Capture  Predicted  With  the  Particle  Tracking  Model 

The  one-way  coupled  particle-tracking  model  (PTM)  was  used  to  predict  the  trajectory 
and  capture  time  for  spherical  bubble  capture  by  a  Gaussian  vortex,  I An  implicit 

second  order  numerical  scheme  was  used  to  integrate  the  radial  and  tangential  components  of 
equation  1  with  a  time-stepping  procedure.  The  case  of  non-cavitating  bubbles  was  first 
considered.  Bubbles  of  varying  sizes  were  released  at  varying  radial  positions  away  from  the 
vortex  axis,  ,  with  an  initial  velocity  equal  to  that  of  the  fluid  at  the  bubble  center.  A  bubble 


was  considered  captured  when  it  intersected  the  annulus  defined  by  1/4  of  the  core  radius,  which 
approximately  corresponds  to  the  region  of  lowest  core  pressure.  In  these  simulations,  the 
bubble  remained  spherical,  which  is  equivalent  to  We  «  1 . 

For  bubbles  very  far  away  from  the  vortex,  we  can  assume  that  the  bubble  accelerations 
and  lift  forces  are  negligible  and  that  there  exists  a  balance  between  the  pressure  gradient 
induced  force  and  the  drag  force  on  the  bubble,  VgVP  “  •  If  we  assume  that  the  largest 

component  of  the  bubble  velocity  is  in  the  radial  direction,  Ug^,  the  force  balance  far  from  the 


The  predicted  capture  times  are  shown  in  Figure  3  for  both  equation  22  and  the  results  of  the 
PTM.  The  approximate  capture  time  under-predicts  the  capture  time  compared  to  that  of  the 
PTM,  and  the  extent  of  the  under-prediction  increases  with  increasing  bubble  size.  However,  the 
relationship  reveals  the  trends  in  capture  time  with  r^lr^,  rc/Ro,  and  Re^,.  The  pressure 
gradient  induced  force  increases  with  R\  while  the  drag  force  increases  with  R^.  Consequently, 
in  regions  of  slow  bubble  acceleration,  smaller  bubble  will  take  much  longer  to  move  toward  the 


vortex  center. 


Once  the  bubbles  come  closer  to  the  vortex  axis,  the  lift  and  added  mass  forces  contribute 
significantly  to  the  bubble  motion.  Figure  4a  and  4b  shows  the  capture  time  for  the  case  with 
Re^  =  Mcrc/v  =  3xlO'’,  which  is  equivalent  to  the  experimentally  measured  vortex.  The 
solution  is  shown  for  the  four  proposed  lift  coefficients.  The  force  on  the  bubble  resulting  fi'om 
the  local  pressure  gradient  drives  the  bubble  toward  the  vortex  axis,  while  the  drag,  lift,  and 
added  mass  forces  act  to  retard  the  bubble  motion. 

Variation  in  lift  coefficient  leads  to  some  change  in  the  capture  time,  although  the  effect 
is  secondary.  We  can  scale  the  average  bubble  Reynolds  number  as 

Re^  ~  =  Re^  (23) 

tc  V/  tc  Tc 

Here  r^/rc~l,  Rec~10'‘.  Then  the  bubble  Reynolds  number  will  be 

Re^  ~  10(Ro  /rc)"*,  and  the  bubble  Reynolds  number  will  decrease  with  increasing  capture  time 
and  decreasing  bubble  size.  The  lift  coefficients  proposed  by  Saffman  (equation  7)  is  valid  for 
Rcg  <  20 ,  and  the  largest  bubbles  approach  this  limit.  The  relationships  proposed  by  Dandy 
and  Dwyer  (equation  8)  and  Sridar  and  Katz  (equation  9)  are  valid  for  somewhat  higher  bubble 
Reynolds  numbers,  up  to  ~100.  The  highest  Re^  will  occur  near  the  point  of  capture. 

The  highest  levels  of  vorticity  occur  for  r/r^  <2,  along  with  the  highest  bubbles  slip 
velocities.  The  capture  time  when  the  bubble  is  within  rlr^  <2  is  on  the  order  of  unity,  making 
the  relative  velocity  ~  .  The  non-dimensional  shear  in  the  core,  a^,  can  then  be  scaled  as 

^  _  |<^c|^o  2U(~Rq  1  ^  2Rq 

^  Iu-IbI  ''c  “c  'b 


(24) 


The  highest  local  values  of  occur  within  the  vortex  core,  however.  For  the  smallest  bubbles 
considered  here,  «1.  For  the  largest  bubbles  considered  here,  ~0.1  to  1.  The  lift 
coefficients  proposed  by  Auton  (Equation  6)  and  Dandy  and  Dwyer  (Equation  8)  are  for  the  case 
of  a  <1,  making  them  less  appropriate  for  the  final  portion  of  the  capture  process  of  larger 
bubbles.  Figure  5  shows  the  variation  of  the  capture  time  with  the  vortex  Reynolds  number  for 
the  case  of  /Rq  =  50,  using  the  lift  coefficient  of  Auton  (equation  6).  As  expected,  the  time  for 
bubble  capture  increases  with  increasing  fluid  viscosity  and  decreases  as  ~  l/Re^, . 

6.  Direct  Numerical  Simulations  of  Non-Cavitating  Bubble  Capture 

The  particle-tracking  model  does  not  account  for  bubble  deformation  and  bubble  flow 
interaction  that  may  occur  during  the  final  stages  of  collapse.  It  is  expected  that  these  effects 
would  be  the  most  important  for  relatively  large  bubbles  as  they  near  the  vortex  core.  To 
examine  this  process,  bubbles  with  /Rq  ~  10  were  computed  using  DNS  for  various  Weber 
numbers  and  cavitation  numbers.  In  these  simulations,  the  effective  vortex  maximum  tangential 
velocity  and  core  size  are  u^=5  m/s  and  rc=2.5  mm,  respectively.  The  length  of  the 
computational  domain  is  4  core  radii,  or  10  mm.  The  vortex  Reynolds  number  =  250, 
which  is  about  50  times  lower  than  the  experimentally  examined  vortex.  However,  even  at  the 
lower  Reynolds  number,  the  amount  of  vortex  diffusion  during  the  computation  was  relatively 
small,  with  the  core  growing  less  than  5%  over  the  duration  of  the  computation.  The  kinematic 
viscosity  and  density  of  the  gas  phase  were  1/10  that  of  the  fluid.  The  initial  bubble  radius  was 
varied  as  a  proportion  of  the  core  radius,  and  the  release  position  of  the  bubble  was  varied  along 
the  radius.  Again,  the  bubble  was  given  an  initial  velocity  equal  to  the  mean  of  the  surrounding 
fluid.  Bubbles  were  released  at  y/L  =  0.5  and  varying  x/L,  with  the  vortex  axis  at  x/L  =  ylL  =  0.5. 


A  grid  resolution  study  was  performed  for  the  capture  of  a  bubble  with  IRg  =  8  and  We  =  0.13 
for  a  non-cavitating  bubble.  The  bubble  was  released  at  r^lr^  =0.5,  a  region  of  large  flow 
gradients.  Simulations  with  resolutions  of  96^  128^,  160^  were  compared,  and  bubble  center 
position  at  non-dimensional  time  t*  =  1.77  were  within  1%  of  the  domain  length.  It  was 
concluded  that  a  resolution  of  128^  would  be  sufficient. 

Figure  6a  shows  the  trajectory  of  a  non-cavitating  bubble  with  IRq  =  8  and  We  =  0.13 
for  hf.  =1.0,  and  Figure  6b  shows  the  vorticity  in  the  x-y  plane.  It  is  clear  that  the  bubble 
locally  modifies  the  flow,  and  a  wake  develops  behind  the  bubble.  Figure  7  shows  the  trajectory 
of  three  bubbles  released  at  r^lr^=  1.5,  1.0,  and  0.5  along  with  the  prediction  of  the  PTM. 
Again,  the  bubble  is  considered  captured  when  the  bubble  interface  crosses  the  /4 .  The 
capture  times  predicted  with  the  DNS  for  the  three  bubbles  are  t*  =  4.61,  2.22,  and  0.74.  These 
are  compared  with  t*  =  4.51,  2.74,  and  0.85  for  the  PTM.  Here,  the  density  ratio  and  vortex 
properties  of  the  DNS  calculation  were  used  in  the  PTM,  and  the  lift  coefficient  of  Auton  was 
used.  The  low  Weber  number  leads  to  relatively  little  bubble  deformation,  and  the  capture  time 
are  within  20%  of  those  from  the  PTM.  This  variation  is  on  the  order  of  the  changes  in  the  PMT 
prediction  that  result  from  the  use  of  the  different  lift  coefficients. 

Figure  8a  and  8b  show  the  effect  of  varying  Weber  number  of  the  trajectory  and  capture 
time  for  two  different  release  positions,  =  1.5  and  1.0.  The  time  decreases  with  increasing 
Weber  number,  implying  that  increased  bubble  deformation  leads  to  faster  bubble  capture.  The 
capture  times  computed  with  the  PTM  and  the  DNS  for  Ir^  =  0.5  are  similar.  But  the  PTM 
under-predicts  the  capture  times  for  the  bubble  released  at  1.5,  suggesting  the 

accumulated  effect  of  varying  lift  coefficient  on  the  final  capture  time.  Figure  9a  and  9b  show 
the  effect  of  varying  bubble  size  on  capture  time,  with  the  larger  bubble  with  IRq  =  8  being 


captured  faster  that  the  bubble  with  =12  for  We  =  0.1.  The  PTM  under-predicts  the 

capture  time  for  both  cases,  suggesting  that  the  lift  and  drag  coefficients  of  the  PTM  are  too 
small  when  the  bubble  is  close  to  cature,  even  though  the  lift  coefficient  used  is  the  largest  of  the 
four  considered.  A  third  trajectory  was  calculated  with  an  arbitrarily  increased  drag  coefficient 
chosen  to  match  the  DNS  trajectory.  The  increase  in  the  drag  coefficient  was  1.5  times  for 
'c  1^0  ^  ^  ^  times  l^o 


7.  The  Capture  of  Cavitation  Nuclei 

The  capture  of  potentially  cavitating  nuclei  are  now  considered.  As  the  gas  bubbles 
experience  reduced  pressure,  volume  of  the  bubbles  will  increase  isothermally.  Bubbles  far  from 
the  low-pressure  region  will  grow  quasi-statically  and  isothermally  {k=  1),  with  the  radius  given 


by 


Pv  -  Pir)  +  Pa, 


2o  Rq 

~R^W) 


=  0 


(25) 


The  radius  of  the  bubble  will  grow  quasi-statically  until  p{r)  ^py.  As  the  pressure  is  decreased 
further,  the  bubble  can  begin  to  grow  rapidly.  A  stability  analysis  of  equation  10  yields  the 
critical  radius  R.(r)  above  which  bubbles  will  grow  explosively  for  a  given  liquid  tension  (see 
Brennen  (1995)  for  a  complete  discussion): 
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Vaporous  cavitation  can  occur  when  the  pressure  at  the  vortex  axis  is  less  than  or  equal  to  the 
vapor  pressure,  Pc^Py  From  the  relationship  for  the  core  pressure  (equation  16)  and  our 


/g,(r)  g,We  ( _ 1 _ 

^  \f{rlrc)lr}^-G^IOi 

Small  bubbles  are  expected  to  move  toward  the  vortex  axis,  quasi-statically  changing  volume 
until  they  reach  a  radial  position,  r^,  where  /?,(?;)  =  R(r^)  when  unstable  bubble  growth  will  occur. 

Consequently,  the  parameters  for  nuclei  capture  are-^.,  u^.,  Re^,  r^,  Rq,  We,  and  tT„. 
We  will  consider  the  bubble  capture  for  the  experimentally  examined  vortex,  thus  setting  the  first 
three  parameters.  The  nuclei  to  be  considered  are  r^lR^^  50,  200,  500,  and  1000 
corresponding  to  bubbles  of  110,  28,  11,  and  6  micron  radii.  This  corresponds  to  We  =  26,  6.5, 
2.6,  and  1.3  for  air  bubbles  in  clean  water.  We  will  consider  freestream  cavitation  numbers 
ranging  from  o^lo^  =  1,  7/8,  3/4,  1/2,  and  1/4.  It  should  be  noted  that,  in  practice,  it  is  difficult 
to  sustain  large  tensions  in  the  vortex  core  for  steady  line  vortices.  When  the  tension  is 
sufficiently  high,  even  the  smallest  nuclei  will  cavitate,  and  the  vortex  core  will  fill  with  an 


annulus  of  vapor. 


Solutions  for  /?,.(r/r^)and  Rirlr^)  were  computed  for  the  four  conditions  described 
above.  The  vertical  dashed  line  is  the  radial  position  when  the  fluid  pressure  reaches  the  vapor 
pressure.  Figure  10  shows  the  case  of  IRq  =  50,  the  largest  bubbles.  As  the  bubbles  come 
closer  to  the  core,  the  equilibrium  radius  changes  by  a  factor  of  up  to  3  times  before  inception. 
All  of  these  bubbles  are  initially  large  enough  to  cavitate  when  they  pass  into  the  fluid  below 
vapor  pressure,  even  before  the  added  volume  increase  as  they  approach  the  core.  Figure  11 
shows  the  case  v^IRq  =  200,  and  Figure  12  shows  the  case  r^^IRg  =  500.  The  quasi-static 
radius  change  of  these  bubbles  is  even  smaller,  about  2  times.  But,  even  with  this  increase,  the 
nuclei  are  still  barely  big  enough  to  cavitate.  Once  they  pass  into  the  region  of  tension,  they  will 
continue  to  grow  and  will  then  cavitate  after  reaching  the  core  for  the  lowest  cavitation  numbers. 
Finally,  for  the  case  of  IRq  =  1000  in  Figure  13,  the  bubbles  change  their  radii  the  least  as 
they  approach  the  core,  and  these  small  bubbles  will  only  cavitate  when  the  core  pressure  is  in 
strong  tension  and  after  the  bubbles  have  grown  quasi-statically. 

The  capture  time  of  the  cavitating  nuclei  will  not  differ  significantly  from  that  of  the  non- 
cavitating  case  for  bubbles  that  originate  far  from  the  vortex.  First,  the  capture  criteria  will 
differ,  since  the  bubble  can  be  considered  captured  once  it  begins  to  grow  rapidly,  and  this  can 
occur  at  a  radius  that  is  farther  than  rc/4,  the  criteria  chosen  for  the  non-cavitating  capture. 
Second,  as  the  bubbles  begin  to  grow  quasi-statically,  the  rate  of  capture  will  increase.  Third,  as 
the  bubbles  begin  to  grow  dynamically,  the  added  mass  term  due  to  volume  changes  (equation  4) 
will  retard  the  bubble’s  inward  motion,  but  the  relatively  slow  rate  of  the  bubble  volume  change 
during  the  quasi-static  growth  of  the  bubble  makes  this  term  negligible  up  until  the  rapid  bubble 
growth.  All  of  these  effects  are  relatively  small,  leading  to  a  decrease  in  the  capture  time  of 
At*  ~  1  for  the  larger  bubbles.  Otherwise,  the  capture  time  for  the  smaller  bubbles  is  not 


substantially  different  between  the  cavitating  and  non-cavitating  cases,  especially  for  bubbles 
starting  far  from  the  vortex  core. 


Figure  14  shows  the  DNS  solution  for  cavitating  bubble  volume  changes  during  the  last 
part  of  the  bubble  capture.  The  simulation  is  for  a  bubble  with  an  initial  r^.  IRq  =  8,  We  =24, 
released  at  =  0.6.  The  cavitation  number  is  cr„  =  1.0  based  on  the  pressure  imposed  at  the 
computational  boundary.  Once  the  bubble  begins  to  grow,  it  is  pulled  in  toward  the  vortex  axis 
and  begins  to  deform.  Figure  15  shows  the  volume  history  of  the  bubble  for  a’„=1.0  and  2.0. 
Once  a  spherical  bubble  begins  to  grow  rapidly,  the  asymptotic  rate  of  radius  growth  is 
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Consequently,  the  volume  should  increase  as 
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Also  plotted  in  figure  15  is  the  volume  growth  predicted  by  equation  35,  where  t*  =  1*^- with 
the  offset  =0.18  and  0.45  for  a„=1.0  and  2.0,  respectively.  These  offsets  are  necessary  to 
account  for  the  time  when  the  computed  bubble  begins  to  grow  rapidly.  Once  the  computed 
bubble  begins  to  grow,  the  growth  rate  is  approximately  that  of  the  asymptotic  value,  although 
the  growth  rate  increases  once  the  bubble  is  large  and  captured  by  the  vortex. 


8.  Experimental  Observation  of  Bubble  Capture  and  Growth 

The  capture  of  bubbles  was  examined  experimentally  for  the  vortex  described  in  Section 
4.  The  size  of  the  bubble  created  varied  with  the  laser  power  and  the  freestream  cavitation 
number.  However,  it  was  not  possible  to  precisely  control  the  size  of  the  bubble  created  with  the 


laser.  Figure  16a  and  16b  shows  8  images  of  a  laser-induced  bubble  as  it  is  captured  by  the 
vortex  for  =  3.0.  The  vortex  centerline  and  core  radius  are  marked,  and  the  time  duration 
between  the  images  is  350  microseconds.  The  capture  time  as  measured  for  multiple  bubbles, 
and  the  results  are  shown  in  Figure  17.  The  estimated  uncertainty  in  the  capture  time  is  +/-  6%. 
The  experimentally  determined  capture  times  fall  within  the  computed  capture  times  given  the 
variability  of  the  initial  bubble  size.  The  bubbles  under  consideration  can  be  considered 
cavitating  nuclei,  as  their  volume  does  grow  as  they  enter  the  low-pressure  region  near  the  core. 
However,  as  discussed  above,  the  non-cavitating  PTM  predicts  the  capture  time  with  reasonable 
accuracy. 

9.  Conclusions 

The  particle-tracking  model  can  successfully  predict  the  capture  time  for  small  bubbles 
by  a  line  vortex.  For  the  small  bubbles  starting  far  from  the  vortex,  variation  in  capture  time  with 
differing  choices  of  the  lift  coefficient  are  relatively  small.  The  accuracy  of  the  PTM  is  reduced 
as  the  bubble  moves  to  within  the  core  radius,  and  the  DNS  solutions  suggest  that  the  PTM  will 
under-predict  the  capture  time  in  the  final  stages  of  collapse.  The  effect  of  bubble  deformation 
and  bubble  vortex  interactions  will  be  important  only  at  the  final  stages  of  bubble  capture.  The 
PTM  will  provide  a  usefully  accurate  estimate  of  the  capture  time  using  standard  relationships 
for  lift  and  drag  on  the  spherical  bubble,  especially  for  the  cases  when  the  bubble  starts  far  from 
the  vortex  axis.  The  non-cavitating  PTM  can  be  used  when  considering  the  capture  of  cavitation 
nuclei.  The  bubble  volume  changes  very  slowly  until  the  final  stages  of  bubble  capture. 
Consequently,  it  is  useful  to  separate  the  process  into  the  non-cavitating  capture  of  the  nucleus 
and  then  the  growth  of  the  nucleus  in  the  low-pressure  core  region.  The  growth  of  the  nucleus 


prior  to  cavitation  inception  is  relatively  small,  but  such  bubble  growth  will  lead  to  a  somewhat 
faster  capture  as  the  bubble  comes  near  the  core. 
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Figure  Captions 


Figure  1:  (a)  Schematic  diagram  of  the  hubble-vortex  interaction  experiment.  A  single  laser 
pulse  is  used  to  create  a  cavitation  bubble  in  the  bulk  of  the  fluid.  Images  of  the  bubble  are 
captured  with  a  multi-exposure  digital  camera,  and  the  acoustic  emission  of  the  bubble  is 
captured  with  a  hydrophone;  (b)  photograph  of  the  hydrofoil  in  the  water  tunnel  test  section  with 
developed  tip-vortex  cavitation  with  U„-  10  m/s,  and  cr„  =  1.4,  with  flow  from  left  to  right. 


Figure  2:  (a)  The  average  vector  map  of  the  flow  perpendicular  to  the  vortex  axis  90  mm 
downstream  of  the  tailing  edge  for  (/„=  10  m/s;  (b)  the  measured  tangential  velocity,  u^ir), 
along  with  the  velocity  distribution  for  the  Gaussian  vortex  with  =  5.6  mm  and  =  0.290 
m^/s  and  the  resulting  calculated  pressure  depression,  p(r)  -  p^. 


Figure  3:  Bubble  capture  time,  t* ,  versus  the  release  position,  tr^.,  for  varying  bubble  sizes, 
Rq  computed  with  the  PTM.  Also  shown  is  the  solution  for  the  approximate  capture  time 
from  equation  22. 

Figure  4a:  Bubble  capture  time,  /*,  versus  the  release  position,  r^lr^,  for  varying  bubble  sizes, 
r^-l  Rq,  and  the  lift  coefficient,  cl ,  computed  with  the  PTM. 

Figure  4b:  Same  as  figure  4a. 

Figure  5:  Bubble  capture  time,  t  * ,  versus  the  release  position,  r^lr^,  for  varying  Renolds 
number.  Re ,  computed  with  the  PTM. 

Figure  6:  (a)  The  DNS  simulation  of  a  non-cavitating  bubble  being  captures  by  a  line  vortex. 
//?o=  8,  We  =  0.13  and  released  at  rjr^=\.  (b)  The  DNS  simulation  of  the  vorticity  contours 

and  the  streamlines  at  t*=0.85  with  the  same  conditions  as  in  figure  6a. 

Figure  7:  The  trajectory  of  a  non-cavitating  bubble  for  varying  the  release  position, 
computed  with  the  PTM  and  DNS.  The  capture  time  criterion  is  0.25  *  . 

Figure  8a:  The  trajectory  of  a  non-cavitating  bubble  for  varying  the  Weber  number.  We , 
computed  with  the  PTM  and  DNS.  The  capture  time  criterion  is  0.25  *  . 

Figure  8b:  Same  as  figure  8a 

Figure  9a:  The  trajectory  of  a  non-cavitating  bubble  computed  with  the  PTM  and  DNS.  Also  is 
shown  the  trajectory  of  a  non-cavitating  bubble  with  an  increase  drag  coefficient  by  150%, 

PMT  * .  The  capture  time  criterion  is  0.25  *  . 


Figure  9b;  The  trajectory  of  a  non-cavitating  bubble  computed  with  the  PTM  and  DNS.  Also  is 
shown  the  trajectory  of  a  non-cavitating  bubble  with  an  increase  drag  coefficient  by  200%, 

PMT  * .  The  capture  time  criterion  is  0.25  *  . 

Figure  10:  The  quasi-static  growth  of  the  bubble  radius  is  plotted  as  a  function  of  distance  from 
the  vortex  axis  for  a  bubble  with  an  initial  size  /i?Q=  50,  We  =  26  for  varying  vortex  cavitation 

numbers  (solid  lines).  Also  plotted  is  the  solution  of  equation  33,  which  yields  the  critical  radius 
for  rapid  bubble  growth  as  a  function  of  I for  varying  cavitation  numbers.  The  vertical 

dotted  lines  occur  when  p{r)  =  for  a  given  .  The  pressure  field  is  that  of  the  experimental 
vortex. 

Figure  11:  Same  as  Figure  10  but  with  r^l R^=2QQ,  We  =  6.5. 

Figure  12:  Same  as  Figure  10  but  with  rJR^=  500,  We  =  2.6. 

Figure  13:  Same  as  Figure  10  but  with  I R^-  1000,  We  =  1.3. 

Figure  14:  The  DNS  simulation  of  a  cavitating  bubble  being  captures  by  a  line  vortex,  “ 

8,  We  =  24,  a„  =  1.0.  The  bubble  was  initially  in  equilibrium  and  released  at  =  1 . 

Figure  15:  The  rate  of  volume  growth  for  cavitating  bubbles  captured  by  a  line  vortex  from  the 
DNS  simulation  for  varying  We  and  a„ .  /i?o=  4,  and  the  bubble  is  initially  in  equilibrium 

and  released  at  /  r^  =  1.  Also  plotted  is  the  offset  asymptotic  volume  growth  of  a  cavitating 
spherical  bubble  given  by  equation  35. 

Figure  16a:  8  images  of  a  laser-induced  bubble  as  it  is  captured  by  the  vortex  for  cavitation 
number,  o„  =3.0.  The  release  position  of  the  bubble,  /  a;  =1 .  The  delay  between  each  image  is 

350  microseconds.  The  capture  time,  r  *=1.25.  The  bubble  size,  /i?o=6. 

Figure  16a:  Same  as  Figure  16a  but  with.  The  capture  time,  r  *=0.9.  The  bubble  size, 

r,//?o=10. 

Figure  17;  Bubble  capture  time,  t  * ,  versus  the  release  position,  for  varying  bubble  sizes, 

r^/Rp,  and  the  lift  coefficient,  cl ,  computed  with  the  PTM.  Also  shown  is  the  experimental  data 
for  varying  the  release  position,  r^tr^. 
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Figure  1:  (a)  Schematic  diagram  of  the  bubble-vortex  interaction  experiment.  A  single  laser 
pulse  is  used  to  create  a  cavitation  bubble  in  the  bulk  of  the  fluid.  Images  of  the  bubble  are 
captured  with  a  multi-exposure  digital  camera,  and  the  acoustic  emission  of  the  bubble  is 
captured  with  a  hydrophone;  (b)  photograph  of  the  hydrofoil  in  the  water  tunnel  test  section  with 
developed  tip-vortex  cavitation  with  =  10  m/s,  and  C7„  =  1.4,  with  flow  from  left  to  right. 
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Figure  2:  (a)  The  average  vector  map  of  the  flow  perpendicular  to  the  vortex  axis  90  mm 
downstream  of  the  tailing  edge  for  t/„=  10  m/s;  (b)  the  measured  tangential  velocity,  u^ir), 
along  with  the  velocity  distribution  for  the  Gaussian  vortex  with  =  5.6  mm  and  F^  =  0.290 
m^/s  and  the  resulting  calculated  pressure  depression,  p{r)  -  p„. 
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Figure  4a:  Bubble  capture  time,  t*,  versus  the  release  position,  for  varying  bubble  sizes, 

r^l  Rq,  and  the  lift  coefficient,  cl ,  computed  with  the  PTM. 
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Figure  5:  Bubble  capture  time,  t  * ,  versus  the  release  position,  Ir^,,  for  varying  Renolds 

number,  Re ,  computed  with  the  PTM. 
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Figure  6:  (a)  The  DNS  simulation  of  a  non-cavitating  bubble  being  captures  by  a  line  vortex. 

=  8^  We  =  0. 1 3  and  released  at  r,  /  =  1 .  (b)  The  DNS  simulation  of  the  vorticity  contours 

and  the  streamlines  at  t*=0.85  with  the  same  conditions  as  in  figure  6a. 


Figure  8a:  The  trajectory  of  a  non-cavitating  bubble  for  varying  the  Weber  number,  We, 
computed  with  the  PTM  and  DNS.  The  capture  time  criterion  is  0.25  *  . 


Figure  9a:  The  trajectory  of  a  non-cavitating  bubble  computed  with  the  PTM  and  DNS.  Also  is 
shown  the  trajectory  of  a  non-cavitating  bubble  with  an  increase  drag  coefficient  by  150%, 
PMT  * .  The  capture  time  criterion  is  0.25  *  . 


Figure  9b:  The  trajectory  of  a  non-cavitating  bubble  computed  with  the  PTM  and  DNS.  Also  is 
shown  the  trajectory  of  a  non-cavitating  bubble  with  an  increase  drag  coefficient  by  200%, 
PMT  * .  The  capture  time  criterion  is  0.25  *  . 
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Figure  10;  The  quasi-static  growth  of  the  bubble  radius  is  plotted  as  a  function  of  distance  from 
the  vortex  axis  for  a  bubble  with  an  initial  sizer^  /  Rq  =  50,  We  =  26  for  varying  vortex  cavitation 
numbers  (solid  lines).  Also  plotted  is  the  solution  of  equation  33,  which  yields  the  critical  radius 
for  rapid  bubble  growth  as  a  function  of  for  varying  cavitation  numbers.  The  vertical 
dotted  lines  occur  when  p{r)  =  for  a  given  cr„ .  The  pressure  field  is  that  of  the  experimental 


vortex. 
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Figure  11:  Same  as  Figure  10  but  with  /;  //?o”  200,  We  =  6.5. 
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Figure  12:  Same  as  Figure  10  but  with  /i?o=  500,  We  -  2.6. 


Figure  14:  The  DNS  simulation  of  a  cavitating  bubble  being  captures  by  a  line  vortex, 
8,  We  =  24,  CT„=  1.0.  The  bubble  was  initially  in  equilibrium  and  released  dX  rjr^=\. 


Figure  15:  The  rate  of  volume  growth  for  cavitating  bubbles  captured  by  a  line  vortex  from  the 
DNS  simulation  for  varying  We  and  cr„ .  r^l R^-  4,  and  the  bubble  is  initially  in  equilibrium 

and  released  at  /r^  =  1 .  Also  plotted  is  the  offset  asymptotic  volume  growth  of  a  cavitating 
spherical  bubble  given  by  equation  35. 
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Figure  16a:  8  images  of  a  laser-induced  bubble  as  it  is  captured  by  the  vortex  for  cavitation 
number,  a,  =3.0.  The  release  position  of  the  bubble,  rjr=\.  The  delay  between  each  image  is 

350  microseconds.  The  capture  time,  t  *=1.25.  The  bubble  size,  rJR^=^6. 
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Figure  16b:  Same  as  Figure  16a  but  with.  The  capture  time,  r*=0.9.  The  bubble  size, 
rJRo-lO. 
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Figure  17:  Bubble  capture  time,  t*,  versus  the  release  position,  r^lr^,  for  varying  bubble  sizes, 
r^l  Rq,  and  the  lift  coefficient,  cl ,  computed  with  the  PTM.  Also  shown  is  the  experimental  data 
for  varying  the  release  position,  //c . 


